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ABSTRACT 


A continuous flow well-stirred tank photochemical 
reactor of known surface area was used to obtain measure- 
ments of the rate constants for the ozonation of free 
cyanide at various UV intensities and 22°C. 

The ozonation of free cyanide, with concentrations 


Ne Vile acy ake ine was carried out 


ranging from 1.14x 10 
at pH 12 with and without UV. The reaction scheme was 
regarded as consisting of the main cyanide oxidation reac- 
tion, a parallel ozone decomposition reaction, and a conse- 
cutive cyanate oxidation reaction. The rate constant and 
stoichiometry of the cyanate reaction were assumed from 
Balyanskii's work, while the ozone decomposition rate cons- 
tant was predicted by Stumm's and Czapski's correlations 
interchangeably. A fast reaction regime model was used to 
decouple the mass transfer from the overall reaction. The 
rate constant for cyanide was found to be independent of UV 
intensity and had an average value of 19086 Masa with an 
uncertainty of 18%. 

A transition from fast to instantaneous reaction 
regime model was also applied to the free cyanide ozonation 
reaction, and the value of the rate constant obtained was 
higher than that of the fast reaction model by 5%. The gas 


phase resistance only accounted for 2% of the total resis- 


tance. 
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The ozonation of ferricyanide, with concentrations 


PM conan os 107°, was studied at four 


ranging from 1.84x10 
different UV intensities and pH of 7 and 12. The rate con- 
stant was obtained by characterizing the reaction to be in 
a transition from slow to fast reaction regime using either 
a flat or quadratic cyanate concentration profile. A reac- 
tion scheme similar to that of the free cyanide ozonation 

was assumed. The rate constant for the ferricyanide ozona- 
tion reaction depended on the UV intensity when the cyanide 


content was lower than 1.5x 107° M or 40 ppm above which a 


PeivemceMiCO MGs sa> was obtained independent of UV inten- 
sity. The rate constant at pH 12, where the ozone decompo- 
Sition accounted for 90% of the total ozone consumed, agreed 
with that at pH 7 where the ozone decomposition was negli- 
gible. The agreement at different pH's supports the model 
proposed. 

The dependence of the rate constant, K, on UV intensity 
and ferricyanide concentration was interpreted by proposing 
a molecular mechanism involving the formation of an ozonide 
intermediate. Two other reported mechanisms, i.e. free radi- 
cal and ferricyanide dissociation mechanisms, were also con- 
sidered. Even though these two mechanisms do not comply with 
the behavior of K observed, the possibility of the presence 
of these mechanisms cannot be excluded. The ozonation reac- 
tion proceeds most likely through a combination of the above 
three mechanisms. Further work to elucidate the mechanism 


and reaction kinetics of ozonation reactions is required. 
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CHAPTER 1 
INTRODUCTION 


Cyanide is one of the most poisonous chemicals known 


(1) 


by man. Wuhrmann discovered in 1948 that the toxicity 
of cyanide comes chiefly from the molecular species HCN, 
but not the ion CN. Generally speaking, the toxicity 
increases with temperature and the age of the victim and 
decreases with the oxygen content in the air and the pH. 


(2) oy 


As reported by Doudorff and Mckee , cyanide causes 
osmoregulatory and excretory failures, internal hemorrhages 
as well as the inhibition of the functioning of oxidase 
responsible for transferring oxygen from the blood to the 
tissues. Exposure to air containing a high concentration 

of HCN (10 ppm by volume) causes the symptoms of acute 
Cyanide poisoning to appear almost instantaneously, namely, 
giddiness, headache, unconsciousness, convulsion and finally 
cessation of respiration as a result of the paralysis of the 


(4) 


respiratory center located in the brain An exposure 
time longer than one hour may cause serious health damage 
and even death. At concentrations higher than 200 ppm, 
fatality happens almost immediately. With rainbow trout 

as a test species, the 48 hour LC, 4 for cyanide (the concen- 
tration of cyanide that kills half of the test organisms in 


48 hours) is 0.07 mg/L. This serves as a guide-line for 


cyanide discharge. The Ohio River Valley Sanitation 
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Commission, U.S.A. recommended the permissible level of 
cyanide in receiving waters (1960) to be 0.025 mg/L, and 
a discharge level of 0.1 mg/L; the same standard had been 
adopted by Canada'!>), 

The use of cyanide for industrial processing is not 
extensive in Alberta?) However, on a national scale, it 
is by no means negligible. Cyanide waste streams result 
from industries such as gold and silver mining, electro- 
plating, photographic processing, coke furnaces and synthe- 
tics manufacturing. The major source of waste cyanide is 


the electroplating industry ‘®"7) 


in which a cyanide bath is 
: +2 RA Uy 

used to hold ions such as Zn and Cd in solution as well 
as to give a uniform and lustrous metallic coating on the 
surface of the object to be plated. "Drag over" of the 
plating solution containing cyanide and metal cyanide com- 
plexes (e.g. zn (CN), + cd (CN) 7°) contaminates the rinsing 
baths. The cyanide in the rinse water ranges from 10 mg/L 


coe 7 00 mg/t '®) , 


For a typical medium size plant, the waste 
contains an average cyanide content of 33 mg/L at a 
discharge rate of 200,000 gallons per month. Another 

major source for waste cyanide is the gold mining industry 
in which cyanide is used to extract the precious metal from 
the ore. The gold mine operated in Yellow Knife, North West 
Territories, discharges 200,000 ~ 300,000 gallons of barren 
liquid per day containing 30~ 50 mg/L of cyanide. This 


cyanide has a major impact on the environment and, in order 
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to comply with the government's environmental protection 
regulations, the cyanide level has to be lowered to 0.1 
mg/j, before the effluent can be discharged. 

The major industrial sources of ferricyanide are 
blast furnaces in steel plants, photographic processes and 
the inadvertent mixing of streams containing ferric ion and 
cyanide. Fe(II) and Fe(III) can both react with cyanide to 


form extremely stable complexes ferrocyanide and ferri- 


ahs) 42 (8) 


cyanide, with the instability constants of 10 arid LO 


respectively. The complexes are so stable that the tests 
for cyanide and iron are both negative in aqueous solution. 
The complexes are not only thermodynamically stable but 


also inert ina kinetic sense as reflected by the slow 


ligand isotope exchange rate (t,=52 hrs at Mee gence 
2 


Nevertheless, ferricyanide and ferrocyanide may undergo 
photo-decomposition under sunlight to liberate the poi- 
sonous compound HCN. In order to keep the free cyanide 
level below the 0.1 mg/L limit, the allowable ferricyanide 
discharge concentration has to be lower than 1.5 TaG)/ Ls 


The current popular cyanide treatment methods include 


Ba eat ton 2 ion Spychpenersion eee biological BeseaeRt tetra 


(13,14) (15) 


electrolytic oxidation carbon adsorption alkali 


Ler and ozonation. Among these methods, 


Cie Loe Os) 


chlorination 
only the last two receive widespread attention 
As a comparison of alkali chlorination and ozona- 

(20) 


tion , chemical oxidation by chlorination requires a 


high concentration of chlorine to drive the reaction towards 
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completeness owing to the lower potential for chlorine 

(E= 1.36 v for Cl., and 2.07 v for ozone). The residual 
chlorine built up in the effluent has to be dechlorinated 
while the residual ozone at 2 ppm level can quickly decom- 


pose by itself into a harmless product, O and also there 


2! 
is no formation of the undesirable final products like the 
chlorinated compounds in chlorination. Furthermore, ozone 
may be generated only as needed so that there is no trans- 
portation or storage problem. On the other hand, the basic 
chlorinator is less expensive than the ozonation equipment, 
but requires extensive support equipment. With the advent 
of modern effective ozone generators, the total equipment 
cost tends to be equalized, and the more economical opera- 
tion of ozone equipment will more than compensate for any 
equipment cost advantage remaining to chlorine. The trend 
of ozonation replacing chlorination can be foreseen. 

While the reaction between ozone and free cyanide as 
well as other metal cyanide complexes is quite spontaneous, 
the reaction between ozone and ferricyanide is slow parti- 
cularly at low ferricyanide concentration. For example, 
ozonating for 75 minutes causes negligible decomposition of 


(21) 


ferricyanide at the 57 ppm level Ferricyanide is 


Classified as a "highly refractory compound" in waste water 
treatment technology, and has a Refractory Index RFI of 
(22) 


270 (cf. 0.41 for free cyanide, a slightly refractory 


compound). RFI is a measure of the difficulty of oxidation 
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of a species, and is defined as 


B° (mg/L) +, (hrs) 
RFI = - 


A° (mg/L) 


m,, (mg /min) t, (min) 


V (L) 


where B° is the cumulative ozone pumped into the solution 
from t=0 to eee per liter of the solution, ty is the 
time required for 50% conversion of the reactant component, 
A, is the initial amount of the reactant component, Mp is 
the ozone mass feed rate to the reactor, and V is the solu- 
tion volume in the reactor. 

While the rate of oxidation of ferricyanide can be 
enhanced by elevated temperature and high ozone concentra- 
tion, neither of these appears to be as effective as UV 


(23,24,20) The UV intensity is an important 


irradiation 
factor and should be optimized. The ozone treatment of 
ferricyanide has been proven to be technically and econo- 
mically feasible for industrial purposes. 

A literature review of the 03/CN- reaction 
shows that the reaction rate is limited by the ozone mass 
transfer from the gas to the liquid phase. Presently, the 
great excess of data lumps the mass transfer and the inter- 
nal reaction processes into one overall effective "rate 
constant". This makes the interpretation of data almost 


impossible. These "rate constants" are valid for a parti- 


cular contactor and at certain narrow gas flow rate and 


( 


(er tay, J (ms Om) om 
———- - € 7 7 
Lv 1 
| one 
“oitulor eft ocdni beqmua chdsn ovixealuayvo sce ‘i | slotiad 
"e aK 
sis #: ,¢ .solyufoe sd2 to seth! we 3s = 7 oo Oe? RSS . 
; ‘Ww * 7 
: . | ae ' 
og : ae 3 
,fcenheamo> tnaicts: etd to. nereyevneg £fe 10? seeiapes om: 
at 


ai «Sf mernanon Jrstores ent te teaver inkates ans “A oh 

egio2a ofa er VY BS , YooCas= Ssat.o7, Sved best ateam SNOSO | 

»todceen e439 gk smulov 

ad ooo shinsyvirrs? to notstsbino 20 exten sdd oliiW 

-s7I0a0ACO arose Adptd Eas stutsteqm?s Sedevets, yd “<a , 

VU as evijoetts as od of en58Sags' S2'sH3 | to tedtiea a0 

imetucans. As 25 ut hematni YWoHif . (05,88 .€E) etiatakbd ea 

90 Fnheanicoa ts Bey veneer, . DBS £m GO su Bisoce fei init 

~onsos bre Yil asset ad a3 ee RAPE at asd shia’ ° 198 ' 

 Readg Ty fst laa nh w52 paeeeeny “er 

HOLINESS % “WAN, 0 erty to weives osdenerth a 

agen stoxo sda yd nesting ai ecez nodtoaas oat 
ody ,yltneserd el ants a die 

~re70i ete Bis olan " 


concentration ranges. They are obviously not suitable for 
design purposes. The mass transfer characteristics, reac- 
tion kinetics and mechanism of the UV-Ozone-Ferricyanide 
reaction still remain untouched. It is the main purpose 

of this research project to further explore the CN. and 
Fe(CN), /0,/0V reaction systems, and hopefully, from the 
experimental measurements, to determine kinetic rate cons- 
tants that are entirely separated from hydrodynamic 
influences. The ozonation reactions are conducted ina 
continuous flow, well stirred reactor with known interfacial 
area and mass transfer characteristics, and the theory of 
simultaneous absorption and chemical reaction is applied to 
determine the parameters that characterize the chemical 


kinetic processes. 
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CHAPTER 2 


LITERATURE REVIEW 


2.1 Chemistry of Ozone, Cyanide and Ferricyanide 
2.1.1 Ozone Chemistry 


Ozone is usually prepared by the action of a static 
electric discharge upon 0, or air and concentrations up to 


(25) 


10% oF are obtainable Different sizes of ozone gener- 


ators are available commercially from companies (such as 


The Welsbach Corp. and PCI Ozone Corp. 2° 


) ranging from 
laboratory scale (3-48 gO,/day) to industrial scale 
(50~1000 1b0,/day) . 

The chemical reactivity of ozone as a hostile com- 
ponent of the environment is widely recognized. Ozone is 
characterized as a very toxic compound and has a charac- 
teristic odor readily detected even at low concentration. 
Ozone can be regarded as an allotropic form of oxygen. 
Because of its oxidizing power, it tends to be reduced to 
the more stable 0. form either by reaction with other 


reducing agents or by self-decomposition. One proposed 


mechanism for ozone oxidation involves the formation of 


the active oxygen a tome.) (singlet oxygen) as an inter- 
mediate, 
ae + O 
0, Sai 0-5 


and this singlet oxygen is believed to be the key component 
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responsible for smog formation in the atmosphere. The 
rate of singlet oxygen production is increased by UV 
irradiation as indicated by the strong absorption band 


in the UV region (2000~3000°A, Hartley band 28) 


). 


Being a powerful oxidant with a standard half-cell 


potential Of. 2,07 .v, 


+ 
+ — 
0, ZHAY+ eZee 0, = H,O 207 Rv. 


ozone has an oxidizing power similar to that of the peroxo- 
disulfate ion and larger than those of permanganate, 
chromate and hydrogen peroxide. Because of its high 
potential, ozone will oxidize most of the metals and non- 
metals to their highest oxidation state by the 'transferr- 
ing' of an oxygen atom through a free radical mechanism. 

The active oxygen atom in ozone can also be transferred to 
the substance to be oxidized through an ozonide intermediate. 


(29) 


Ozone reacts with alkene to give an ozonide which con- 


sists of a 5 member ring, and this ozonide can rearrange 


rapidly to give another different 5 member ring compound 


tele Pomel — ag Nes 
O Da \ / 


which is quite stable at low temperature but becomes 


7N 


C=C 


explosive at high temperature. This ozonide can decompose 


into aldehyde, ketone or carboxylic acid by hydrolysis in 
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the presence of a reducing agent,e.g.Zn, or it may undergo 
hydrogenation by H, with Pd as catalyst. 


The decomposition reaction 


(30) 


is very exothermic (AH=-142 Kcal/mole) and is cata- 


lyzed by catalysts such as platinum, nickel and oxides. 
In the absence of catalyst, the decomposition reaction was 
found to be first order with respect to ozone. The reaction 


rate, r mole/f£-s, can be expressed as (31,32) 


r= = Ky [03] 


where Ky) is the rate constant at a certain temperature and 
PH. The effect of pH on the decomposition rate is still 


controversial. While some workers found the decomposition 


(300733) 
a 


rate to be lower in an alkaline medium Shambaugh 


(31) 


and Melnyk held an opposite opinion. They claim that 


Stumm's correlation can be used to estimate Kp in the pH 
range from 7.6 to 10.4 


Kp = KtoH]°*’? sec? at 20°C (1) 


and k = 242 if [OH ] is expressed in the unit of molar con- 


(323) 


centration. At the same time, Czapski correlated Kp 


and pH in a somewhat different manner 


Ky = 700{0H ] Pace at 250 (2) 
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The dependence of ozone decomposition rate constant 
on PH can be explained by the mechanism proposed by Alder 


ene nett (34) 


K 
n: + o 
+ 
03 HO —> HO, + OH 
K 
+ - 2 
HO. + OH = 2HO 


3 
3 2 


K 


4 
03 + HO, eo HO mar 20, 


K 


= 
HO> + HO. ——— >: H,0 3 0. ° 


The change in ozone concentration can be expressed by the 
following reaction rate equation 
Ex 
nee 
dt 


+, -_b 1. 
= - 3 Kk, [HO,] *[0H Apel e eels : 

Both correlations indicate the positive effect of 
hydroxyl ion concentration on the decomposition rate, and 
yield similar Kp values at the same pH. The fact that the 
decomposition of ozone becomes significant at higher pH 
imposes a limit on the use of ozone as an oxidant in the 
chemical treatment of waste water at high pH. For instance, 
Stumm and Czapski predict a Ky of 7.65 and 7 respectively 
at pH 12, which means that half of the ozone will be 
destroyed by decomposition in 0.1 sec (half life); thus a 
large percentage of ozone will be wasted by decomposition 
especially when the reaction between ozone and the pollu- 


tants is slow, as in the case of ferricyanide-0, reaction. 
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Owing to its high reactivity, the popularity of ozone 
is only second to that of chlorine in water treatment. 
Numerous articles have been published in that Fea eo 
In summary, the typical applications of ozone include the 
following: 

a) Treatment of industrial waste water containing 


cyanide, phenol, sulfur compounds, surfactants and 


polymers; 

b) Disinfection of process water and deionized water; 

-) Algaecide for cooling towers; 

d) Bleaching agent for textile and paper industries; 

e) Odor control in water treatment plants and industrial 
plants; 

-} Treatment of domestic waste water (tertiary treatment); 

g) Disinfection of biologically treated waste water. 


2.1.2 Cyanide and Ferricyanide Chemistry 


Cyanide is one of the most poisonous chemicals known 
by man. Being a strong complexing agent, cyanide has the 
ability to form stable complexes with virtually every heavy 
metal in the transition metal block. In aqueous solution, 


it can undergo hydrolysis easily to give a volatile weak 


acid, HCN, which has a dissociation constant of 4.93*x 10719 


ara 52oC (32d. 


Fe(II) and Fe(III) can both react with cyanide to 


form extremely stable complexes ferrocyanide and ferri- 


: - —-42 (8 
cyanide, with the instability constants of 10 35 and 10 (8) 
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respectively. The complexes are so stable that the tests 
for cyanide and iron are both negative in aqueous solution. 
The complexes are not only thermodynamically stable, but 
also inert in a kinetic sense as reflected by the slow 
isotope exchange rate (ty = 52 nrsvat BH oes ae Both 
complexes may undergo photo-decomposition under sunlight 
to liberate the poisonous compound HCN. Devonshire and 


(38) 


Weiss found that UV irradiation of ferrocyanide and 
other anions leads to an unstable species with an absorp- 
tion spectrum identical with the spectrum of the hydrated 
electron (~700 nm). The mechanism of the photo-chemical 
primary process in aqueous ferrocyanide solution is a 
predissociation of the ferrocyanide ion leading to the 


trapping of an electron at some distance from its parent 


center. 


= 2505 

Fe (CN), + hv — Fe (CN) excited ion 
—4* =3 

Fe (CN) ¢ —> Fe(CNn) ¢ +e 


en} orap «—>.ce trapped electron 


t 


e, may react with any oxidants in aqueous solution, for 


ee 
: + (39) 
instance, H30 , dissolved oxygen etc. Ohno 


investiga- 
ted the other aspect of photo-chemical dissociation of 
ferrocyanide in aqueous solution by using a high pressure 
Mercury lamp (366 nm) to irradiate TOmePett Oe M potassium 
ferrocyanide solution in a quartz vessel. The following 


reactions are observed: 
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An hv a = 
Fe (CN) ¢ + HO aoe Fe (CN), (OHA) + CN 


CN” ee HO <— HCN + oun 


The species Fe(CN),OH, is identified by its UV absorption 
peak which ranges from 410 nm at pH 1.14 to 450 nm at 
pH=6.5. Ferricyanide reacts with UV in a similar 


fashion» 


, a dissociated product Fe (CN), (OH) 5” is pro- 
duced upon UV irradiation at the wavelengths of 254, 313, 


405 nm. 


2.2 Physical Properties of Ozone 
2.2.1 Solubility of Ozone in Aqueous Media 


The values of Henry's law constant are not consistent 


in the literature. A value of 67680 atm/mole/cc at 20°C 


(41) whereas 79400 atm/ 


mole/cc at 20°C is from Matrozov's pCbIeet ona The 


is quoted from Perry's Handbook 


temperature coefficient of Henry's law constant AH/AT is 
provided by Matrozov's work. The coefficient in the range 
from 20 to 25°C is found to be 2870 atm/mole/cc per °C by 


interpolation. 


Pee DITLUSIVItY 
The values of diffusivity of oF in water can be 
estimated from the Wilke and Chang correlation or obtained 


from experimental works. The values from various sources 


are as follows: 
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2.3 Ozone/Cyanide Reaction 


The ozone and free cyanide reaction has been studied 


in a batch reactor with continuous 0. flow by Sondak and 


3 
(44,45) (46) | palyanskii et a1. (47), 


(49) 


, Khandelwal et al. 
(48) 


Dodge 


Kandzas et al. as well as Matsuda et al. For the 


most part these studies are unsatisfactory because they 
fail to decouple the ozone mass transfer step from the 
internal reaction step. For the ozonation reaction 
expressed as CN” +0,(g) > 9 (2) + CNO, the rate law can 


be written as 


TCN 


—s. ' n 
at = Ke LONy 


when the partial pressure of 03 is constant. Sondak and 
Dodge found that the overall oxidation rate constant K' 

is almost independent of temperature, which implies a 
diffusion control mechanism. They reported that the reaction 
is vee order with respect to cyanide for cyanide concentra- 
tions greater than 4 ppm, and pseudo first order for con- 
centrations between 2 to 4 ppm. In addition to the 
variation of reaction order, K' was found to depend on the 


gas flow rates, and the ozone and cyanide inlet concentra- 


tions in a very complicated way. In summary, Dodge and 


Sondak found 
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n= 0 [CN] . > 4 ppm 
n= 1 4 ppm > [CN] . > 2 ppm 
K'= ) 


f£(Fg, [CN] ;,P, 


However, Khandelwal et al. found that n is equal to 
1/3 over a large concentration range. 

Balyanskii et al. investigated the CN /0, reaction 
in more detail at pH 11-12. The oxidation of cyanide ions 
by ozone was described by two consecutive reactions 
(unbalanced ) 


Gr (LPernrss 0, (g) BR MEONCS (a) 


CNO™ (2) + 4.5 03(g) SOT NO MON LS NRE HCO DEAE) 


3 2 3 


Both processes are pseudo first order with respect to CN. 
and CNO~ respectively. They report that the rate of cyan- 
ide oxidation is about 5.1 times higher than that for 
cyanate (CN ,0.634 mint, CNO7,0.124 min? at Fg = 0.34 
L/min). From reaction (a) and (b), the ozone consumption 
is found to be far above the value predicted by simple 
stoichiometry. This may be due to the fact that the possi- 
bility of ozone decomposition has been overlooked. 


(36) also indicated that the ozonation of cyanide 


Selm 
ion is very fast with no ozone existing in the solution 
until the cyanide is completely oxidized. The overall 


conversion of cyanide was determined by how fast the ozone 


can be transferred from the gas phase to the liquid phase. 
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(50) 


Heist suggested that ozone oxidation of cyanide 
involves an electrophilic reaction either initiated by a 
singlet oxygen or proceeded through an ozonide intermediate 
as Mentioned in Section 2.2.1 Ozone Chemistry. It is 
Still not clear whether the ozonation of cyanide proceeds 
through a free radical, the formation of ozonide or a com- 
bination of the two. 

It has been mentioned in Section 2.1.2 Cyanide and 
Ferricyanide Chemistry that the reaction between ozone and 


(21) 


ferricyanide is slow , and it can be enhanced by UV 


1. (20723724) ctuagied the uv- 


irradiation. Prengle et a 
Ozone-Ferricyanide reaction in both laboratory and pilot 
scale reactors. A batch reactor with a continuous flow of 
ozone was used to study the reaction characteristics, and 
the gas phase was well dispersed to enhance the ozone mass 
transfer. The concentration of cyanide was recorded at 
different time intervals and the effect of UV, temperature 
and concentration of ozone studied qualitatively by com- 
paring the cyanide concentration vs time curves running at 
various conditions. However, basic kinetic parameters and 
the reaction mechanism cannot be obtained by this approach 
because no attempt is made to decouple the mass transfer 
from the reaction step. Generally speaking, Prengle et al. 


found that the reaction is favored by high UV intensity, 


high temperature and high ozone partial pressure. 
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A six stage continuous flow reactor system has been 
developed by the Houston Research eroup voce in which the 
concentration of ferricyanide can be reduced from 4000 ppm 
at the first stage down to 0.3 ppm at the sixth stage. 

UV irradiation (75 min) is only applied to the last stage 
where the ferricyanide concentration is below 20 ppm. 


There seems to be no apparent advantage to use UV at ferri- 


cyanide concentrations higher than 20 ppm. 
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CHAPTER 3 


THEORY 


3.1 Background for Mass Transfer Modelling '°1/52+93) 


Each gas/liquid two phase reaction involves the 
transfer of a component, A, from the gas phase to the 
liquid in which the reaction takes place and simultaneous 
transfer and reaction in the liquid phase. The liquid 
phase reaction may be a very complicated one consisting 
of a series of parallel and sequential reactions involv- 
ing A and B, the liquid reactant, as well as a number of 
intermediates. The technique of studying a complicated 
gas-liquid system is to classify the system into a reaction 
regime and then reduce the system into a number of simpler 
elementary steps, from which useful kinetic information 
can be obtained by the application of models. The film 
model and the penetration model serve equally well for 
this purpose. 

Consider a case where the gas component A trans- 
fers into the liquid and reacts with B irreversibly 


according to the following reactions 
AA) ee Ed 
ns 
Gh: REN eas Jee 


If the reaction is first order with respect to both A and 


B, the reaction rate can be expressed as 
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By material balance on A and B, the following differential 


equations can be written 


oH aC, 
Da ay oar Ways ese : (4) 
| | | 
Diffusion accumulation reaction 
a°c, 9C, 
Da ay 7 neem oe LS Aes : (5) 


At steady state, the equations can be simplified as 


3 Cy 
D Ni Gace Ed & (6) 
A 5X2 AB-A~B 
and 
CP 
D ail Okra Od © 2 (7) 
B 5x2 AB-A~B 


In the film model, it is assumed that a fictitious 
stagnant film exists near the surface, outside which there 
is no concentration gradient. The thickness of the film 
6 is not directly measurable, it is even questionable 
whether such a film exists at all. Nevertheless, the 


thickness is related to a well defined physical quantity 


Kor the mass transfer coefficient, by 
D 
me A 
where 
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The penetration model takes as its basis the replace- 
ment at intervals of elements of liquid at the surface by 
the liquid from the interior which has the local mean 
bulk composition. While the element of liquid is at the 
surface and is exposed to the gas, it absorbs as though it 
were quiescent and infinitely deep. The rate of absorp- 
tion is a function of the time of exposure of the element 
(contact time 6), and can be obtained by solving the 
unsteady state partial differential equations for material 
balance (Equations (4), (5)). The physical mass transfer 
coefficient K, is directly related to the contact time 6, 


0 


which depends solely on the hydrodynamic conditions: 
D 
A 
ens e2 8 R 10 


In case of mass transfer with chemical reaction, the trans- 
fer of A across the surface will be augmented owing to the 
additional driving force on the gas molecule. Thus for 


simultaneous mass transfer and chemical reaction 


= = ; Be 
ie. SK, (Cag Cro) (Ee 


K, is always greater than Ko: The extent of the increase 


in mass transfer because of chemical reaction is expressed 
quantitatively by the ratio of Ky to Kor known as the 


enhancement factor 


K 
age Sr ae (12) 
n Ko 
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borrowed conceptually from the penetration model have to 


Two more important parameters t 


be introduced in order to define the reaction regimes: 


A 
0 
uy 
t Sh seme : (14) 
B Kapp 


The diffusion time, t is the average life of the surface 


Db” 
element; in other words, the interval of time among 


successive mixing processes which makes the concentration 


within the liquid element uniform. The parameter, ty: 


depends solely on the hydrodynamic condition. The reaction 


LImMme,) c is the length of time required in order that 


R’ 
the reaction may proceed to an appreciable extent. tp 
is only affected by the kinetics of the reaction but not 
the hydrodynamic condition. The solution of the simul- 
taneous ordinary differential equations, Equations (6) and 
(7), for the film model depends on the boundary conditions 
at the surface (x=0) and the thickness of the film (x= 6); 
the boundary conditions are, in turn, governed by the 
relative magnitude of tp and th: Several regimes are 
defined for the sake of convenience and simplifying the 


(51) provides both the cri- 


solution procedure. Astarita 
teria for defining the regimes and the methods for deter- 
mining the effect of the liquid phase reaction on K: 


A slow reaction regime refers to the case where the 


reaction is slower than the mass transfer, i.e. 
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The solution would probably (but not always) be saturated 
with A; the system is regarded as reaction controlled. 
Either using a mathematical approach by setting the appro- 
priate boundary conditions or solely by physical argument, 


it can be shown that K. is equal to K If the reaction is 


L 0° 
in the kineti - 
n c sub-regime where SK Cas Sig VK a peas Ro 
9Cy 
1 = = We =-=- —_ 
4 A Er eMMADN GAS CR io) a ht oe 


Since the reaction is unaffected by mass transfer and Cas 
is well defined (P,/H,) + the reaction can be studied by 
the conventional time-dependent batch reactor, thus the 
value of Kap can be calculated from the slope of the log on 


vs t plot. For a CSTR, the mass balance will be 


0) ns _ 
J += FLAC, = Kapoas’ BoY . Cds) 


In the case of slow mass transfer i.e. baa ty where 
the reaction is fast yet not fast enough to affect the con- 
centration profile of B significantly, and if the reaction 
can be regarded as pseudo first order, the system is said 


to be in the fast reaction regime. The solution of 


Equations (6) and (7) is 
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A system sitting somewhere between the above two 
extremes is said to be in a transition regime, where 
te = th: An analytical solution can be obtained from 


Equations (6) and (7) by the boundary conditions 
=O = = 
x CG 6; P/HA 
+ aa CH= 0 : (19) 
The solution is 


/DxKeac 


K,. - A AB Bo i (20) 
Dake: C 
ental A AB Bo | 
ND 


This general solution can be reduced to the two extreme 


cases by taking the limits 


lim K. = K slow reaction regime (21) 
L 0 
K,+© 
0 
. — Nhe foe . 7 . 
Ay Kk, DAKa RoR, fast reaction regime (18) 
(4) +0) 
D 


which is in perfect conformity with the solutions of the 
fast and slow reactions. 

If the reaction is faster than that described by the 
fast reaction regime, the concentration of B will not be 
constant any more, and the assumption of pseudo first 
order reaction does not hold; the reaction is said to be 


in a transition from fast to instantaneous reaction regime. 
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The boundary conditions for the differential equations are 


Cy = Cag at x = 0 (22) 
ape. 

ax 

Cn = 0 at x = 6 

Ce ': CB 


An analytical solution of this set of equations is not 


(52) (1948) 


available, but van Krevelen and Hoftizer 
computed an approximate solution and expressed the enhance- 


ment factor E as follows: 


E:. -E 
1 
ae Bony 
2 el : (23) 
i E 
tanh M(=—_7) 
af 
where 
DaKap@Bo D 
wer Pe tp 
0 
pac 
ren eee ee Eo 
1 qDaCns 


q is the stoichiometric factor for the reaction A+qB>P; 
Kap can be solved by trial and error. 
Danckwerts suggested that if 


/M >10 E. 


t D.C 
oe. kK ee alii et) (24) 
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the reaction is in the instantaneous reaction regime, and 


the solution can be reduced to 


E = E. 
1 
K ne 
aes = Pe ALE ss : (25) 
0 ole Not Me 


Similar approach can be applied to the penetration 


model. In the transition regime, the second order reaction 


is again approximated by a pseudo first order one. The 


initial and boundary conditions are 


sete U Fe i Bs lng Cn = Cr0 Cp = a (26) 
tae x= 0 Cy = Cas Ca = Sie 
x > © Cr # +0 C. = C3, 


then Equation (4) can be solved readily. The flux J' is 


found to be (21152), 


oCy exp (-K,,C, t) 
'= = res = Sf 
J Bees ) ee Cag’ ante Pater f¥Fagte, t —— } 
ay AB Bo 
2%) 


where 


X 
| exp (-x7) ax 5 
0 


Taking the average J' over a period of time from t = 0 to 


t = 6, the life time of an element, then it can be written 


as: 
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Despite the minor difference, the result is entirely 
compatible with that of the film model, especially in the 
limiting cases of slow and fast reaction regimes, where they 
are identical. Either one of the models can be applied satis- 


factorily to the O,-ferricyanide-UV reaction under study. 
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The approximate solution for the transition from fast 


to instantaneous reaction regime is provided by Brian 


et oo. 
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Using a two resistances in series model, the overall 


mass transfer coefficient in liquid phase, K can be 


OL’ 


expressed as a series combination of the two resistances 


in the gas phase, KS and the liquid phase, K 


ie 
o-=ptee (31) 
OL L g 


where H has units of atm/mole/cc and ia has units of cm/s 


H'='.67680 + 2870(T, — 20) . 


1 1 RT 
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If the gas phase transfer resistance is very small compared 
with that of the liquid phase, the above equation can be 


reduced to 
K = K F (33) 
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J. Rowley measured the gas phase resistance for 
the same reactor under a similar operating condition, a 
value of 1.53 cm/sec was obtained for es hence the term 
SH oes would be equal to'2.15 cm/sec at 22°C. 

If the reaction regime is well characterized, the 
basic parameters can be obtained by applying the appropriate 


solution for that regime. 


3.2 Modelling of Ozonation Reaction 
3.2.1 Stoichiometry of the Ozonation Reactions 

It has been revealed in Chapter 2 Literature Review 
that the ozone consumption in Balyanskii's work is far above 
the value allowed by simple mass balance. This may be due 
to the fact that ozone decomposition is significant at high 
pH; the additional consumption of ozone by decomposition 
causes a higher apparent ozone/cyanide or cyanate ratio. 
Attributing the excess consumption of ozone to decomposi- 
tion and completing the material balance in the following 


unbalanced overall reactions provided by Balyanskii'*’? 


CN + 1.30, ar ONO S a) 


CNO + 4.503 —— id, NO, + O.15N, + HCO, b) 


the second reaction can be expanded into 2 balanced 
parallel reactions c) and d), one of which is the cyanate 


oxidation reaction while the other is the decomposition 


reaction. 
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By summing up these two reactions, a balanced reaction b) 


(reaction e) can be obtained 


Ont CNOM+t A 0.5.0. o> 0 Si INOM mh iL5<N ot. HCO 


0. H 
Brel; 3 3 2 3 


aa ieee Coe 0.7H* e) 


The cyanate oxidation reaction itself is composed of two 
parallel reactions. Cyanate can be oxidized directly to 


No, or hydrolyzed in water to give NH., which is then 


3 
The oxidation products depend 


3° 


on the residence time of cyanide in the reactor. 


oxidized by ozone to NO 


Based on this analysis it can be assumed that 1 mole 
of cyanate reacts with 3.25 moles of ozone after the cor- 
rection for ozone decomposition has been made. Stoichio- 
Metric ratios of 3.25 for cyanate, 1 for cyanide and 1 for 


ozone self-decomposition are used for all subsequent calcu- 


lation: 
dono 3.25 
Ton = 1 
Het . 
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at 2.2 Reaction Scheme - 0,/CN_ Reaction 

The purpose of this research project is to model the 
CN and Fe (CNY /0,/0V reaction systems and hopefully, from 
the experimental measurements, to be able to produce a more 
meaningful and interpretable value of K which is entirely 
separated from the hydrodynamic influences. In the modell- 
ing the whole 2 phase system has been broken into a number of 
simple consecutive and parallel elementary processes which 
can be handled proficiently by the models mentioned in the 
last section after the proper regime is defined. 

Schematically, the system can be described as a mass 


transfer process in series with a family of elementary re- 


actions 
(g) 
K — 
ec) peels fy aay Sys + CNOe 
3 (f£ ) 3 Ke 2 
in| Ky 5 | onogas ag 
1 0. NO3,N,,HCO, Pa Pi 


where (f) is the mass transfer process; reaction (g) is 
the free cyanide oxidation reaction, and Ke the rate con- 
stant; reaction (c) is the cyanate oxidation reaction, and 
Kp the rate constant; reaction (i) is the ozone decomposi- 
tion reaction, and Ky the rate constant. 


According to Balyanskii's work, reactions 9, C, 1 canbe 


written in a more detailed form 
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The reaction rates are assumed to obey simple rate laws. 


rate of ozone decomposed K, [03] 
rate of cyanide oxidized = K- [03] [CN ] 


rate of cyanate oxidized = K, [0] [CNO ] 


From Balyanskii's work, the ratio of Ke to K, can be appro- 


ximated by 


= eT (34) 
B 
From Stumm's ROCrola tone a 
K. = 242 foa} °c at..20°C : (1) 
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3.2.3 Reaction Scheme - UV/0,/Ferricyanide Reaction 

In contrast to the chemistry of the 0,/CN /UV reac- 
tion, the reaction system of 0,/Fe (CN) 2 /uv is more compli- 
cated. There is little work done on the properties of 
cyano group in the ferricyanide molecule; only the cyanide 


(87,88) oA 


ligand exchange rate has been studied so far 
the destruction of cyanide in the ferricyanide complex has 


not been attempted at all. The work shown below only 
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deals with the mass transfer characteristics and the rate 
of consumption of ozone. No attempt is made to elucidate 
the chemistry of the whole system; further research is 
warranted in this area. 

Similar to the free cyanide case, the reaction 


scheme for ferricyanide can be postulated as 
i] 


K 
0, (9g) — 0, (2) —> 0, + CNO- + Fe(CN) ou”? 


a 5 
Fe (CN)¢ 
Ky Ky [3.25 0 0, | Ky 


150 Final products Fe (CN) 


3 (OH) 2+ CNO +0 


4 2 


etc 


The only difference between the free cyanide and ferri- 
cyanide cases is the assumption of the additional consecu- 
tive oxidation of the ferricyanide molecule and its 
intermediate until all the cyanide groups in the complex 
are exhausted. Fe (CN) ,OH ~, Fe (CN) , BS 


are the postulated reaction intermediates containing 5 to l 


(OH) 5°... Fe (CN) (OH) 


cyanide groups, and Ker K,+--Kj the corresponding rate 
constants. The rate of cyanide destruction can be written 


ass: 


3 


' -3 ' - ' -3 
R= [0] {Kk [Fe (CN) ¢ ] + Ke [Fe (CN) , (OH) Jato. -K, [Fe (CN) (OH). }} 


(35) 
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The subscripts from 1 to 6 are to indicate that these 
cyanides are different chemically, and thus would have 
different reaction rates. In order to simplify the model, 
Ky-++Ke are lumped into a collective term K, which can be 
regarded as the effective rate constant, and then R can be 


rewritten as 
R= K[03] [CN] (36) 
where 


[CN] = i[Fe(CN),(OH)77,], 


|ma38S 


the total cyanide content. The above expression holds for 
one of the following special cases: 

1) All cyanides are identical so that they can react 
with ozone at the same rate. By taking into account the 
stoichiometric ratio, the following relationships between 


individual Ke can be established: 


Ae ae “eg H —— : aC a ts Orie 
Ke: Ko: Ky t e+. ¢ Kj = 6: 5 1 
because the number of cyanides in Fe (CN) 7° is) 6, 
Fe(CN) OH >, 5 and Fe (CN) , (OH) 5”, 42 etc.o Then 
eK! BE cra oe be CN) , (OH) 7. ] 
)_K [Fe (CN) ; (OH). 3] = ()_K; ilFe(CN), aes 
ey i=l 
= K, [CN] (37) 
A Rusk (38) 
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ii) The subsequent oxidations of ferricyanide inter- 
mediates are fast comparing with the oxidation of ferri- 
cyanide, i.e. 


[| i] t 
3°° °K, >> Ke 


Ki +K5,K 
then 


kK = & : (39) 


In this case, the intermediates do not exist. Applying 
the approximation in Equation (36), the complex reaction 
scheme can be reduced readily into a simple one identical 
to that assumed for the ozonation of free cyanide 

except the different rate constant. Because of the 
approximation made in Equation (36) K would be a weak 


function of cyanide concentration. 


woe Moge Liang of “the CN /0,/UV Reaction 
3.3.1 Transition from Slow to Fast Reaction Regime - 
Film Model 
The derivation here is to develop the formulation 
which holds for general cases from slow to fast reaction 
regimes with bulk ozone concentration not equal to zero. 
Referring to the reaction scheme shown in Section 3.2.2 
Reaction Scheme - 0,/CN_ Reaction, the general material 
balances for CNO , CN “and 0, at steady state can be 


expressed by a set of ordinary differential equations: 
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A a2 KC,C, Tate 5 KCala + KRCy (40) 
wor 

D — Pe <p Gait Bs 

B 9x7 f-AB 

a2 Ce 

Da mie: = “K,CAC, + KROnSo 


where Cyr e Co are the concentrations of ozone, cyanide 


B’ 
and cyanate respectively. The equations might be solved 
by a tedious numerical analysis method. However, the 
equations can be simplified tremendously by making several 
bold yet valid assumptions. Providing the reaction to 

be in a fast or slow reaction regime, the cyanide 
concentration can be readily assumed to be constant from 
the interface to the bulk solution (i.e. a flat concen- 
tration profile), so the second order reaction of cyanide 
and ozone can be regarded as a pseudo first order one 
(Section 3.1 Background for Mass Transfer Modelling, 
Equation (5)). Since cyanide is in a much larger amount 
than ozone in solution, the flat cyanide concentration 
profile assumption should be valid as long as the reaction 
is not in the instantaneous reaction regime. The material 
balance of cyanide need not to be considered. To assume 
a flat profile for cyanate is not so easily justifiable 

as that for cyanide because the concentration of cyanate 


in the film has to be higher than that in the bulk solution. 
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However, cyanate is present only in a minute amount compar- 
ed with cyanide for the residence time used in the experi- 
ments (cf. Mixing Pattern section) and since the reaction 
rate of cyanate is only 1/5 of that of cyanide, the error 
associated with the over-simplification will not affect 

the result to any significant extent. The validity of this 
assumption will be discussed in Section 6.1 CN/O ,/UV 
Reaction. The above 3 differential equations can then be 


reduced into one 


—— = = K-CiC.. teecm paket © + KC (41) 


x = 0 Cy = Crs 
xi = 6 Cy = Cro 
where Crol Ceol Cao are the bulk concentrations of ozone, 


cyanide and cyanate respectively. 


‘ ope {Sed 
The solution of the differential equation ee ) 


C. 2 sinh (a,x) + Cas sinh[a, (6-x) ] 


= 42 

Cn a sinh (a, 6) (2) 
where 
eae 
1 Da 

ie is 
Ky = K Ca, + 3.25 KBCo, Ky 
eqraenes 
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a) Mass Balance of Ozone in the Film 


The amount of ozone reacted in the film ry is 


I 


al 


0, coming in - 03 going out 
dC dC 
A 
el ABP eH ieee) ; 
X= 6 aly: AB x=0 


II 
o) 
WN 


From Equation (42): 


aC Cc cosh(a_.6) - Cc. 
weal os Ao 1 As 
( ) 2, | sinh (a, 6) | (44) 


1 


aC C -C cosh(a,. 6) 
AGM AS | (45) 


sinh (a, 6) 


= A 1 
oe r,= sinh (a, 6) [cosh(a,6) - 1] [C,,+C, 1] . (46) 


The amount of ozone reacted in the bulk phase ry ss 


The total amount of ozone reacted, RT, is: 


D,S a, 
RT = ri + ro = Sta Mee - 1] [Cag + Cro! + KC, OY 


b) Mass Balance of Ozone at the Interface 


ac 
The flux, J'*=-- sp, (= *) = total ozone reacted + 


x=0 


ozone flowing out in the liquid phase. 
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c) Bulk Concentration of Cyanate: G 
Sn wyanate? CA, 


Working on the cyanide material balance in a CSTR orl c 


can be shown that 
pe Ce a [oi Cerca] y 49 
ayaa bs 1 Bo fe) ( ) 


where Cai is the inlet cyanide concentration and Cc’ the 


A 
effective ozone concentration averaging over the film and 


the bulk. At steady state 9C, fot = 0, 


vad ' 
Parc... (2 ele KeCn Cav j (50) 


Similarly, the balance on cyanate gives 


ss ' oss ' 
Fil Cel = Koa Cao¥ ek of V ° (51) 


Co f°-Bo A 


Solving Co5 from Equations (50) and (51) 


ee ge Le (52) 


d) Evaluation of RT, the Total Ozone Reacted 
The ozone consumption can be calculated from several 
possible ways: 


VORA is related to R by the rate constants 


' i} ' 
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pp Kya t KpCRoCa + 3-25 KRCQCa ey 
ean Oh oe KC Chan hn. 
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where CA is the effective ozone concentration in the solu- 


tion 
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ii) RT is directly measured from the gas phase ozone 


mass balance, 


RR =sR0:.. (56) 


3.3.2 Transition Regime with Cro Equal to Zero -Film Model 
Equations (47), (48) can be tremendously simplified 


by setting C equal to zero; then a standard solution of 


Ao 
the form 
rane sae Ness es. (57) 
L tanh (a, 6) SCa 5 
can be obtained from the fecetire oe Equations (47) and 


(48) in Section 3.3.1 are replaced by Equation (57) with the 
definition of aye Kir 6 remaining unchanged; then the only 

unknown is Ky which can be solved by Mueller's method a aden 
(Appendix 1). The concentration of ozone drops fairly 
rapidly as the reaction proceeds from the slow reaction 


regime to transition regime. 


3.3.3 Fast Reaction Regime 
The system is simplified a step further by the assump- 
tion of fast reaction regime. The reaction is ina fast 


reaction regime if ss tyr and for this case 
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(58) 


Similarly, replacing Equations (47) and (48) by (58), Ke 


can be solved. 


3.3.4 Transition from Fast to Instantaneous Reaction 
Regime 
For the reaction to be in the fast reaction regime, 


the following criterion should be fulfilled: 


i *D aes Cae ESEINBO, t 
Se aera 


If not, the reaction would be in the transition or even 
the instantaneous reaction regime. In the transition from 


fast to instantaneous reaction regime, an approximate 


52) 


solution is available from Danckwerts for the following 


differential equations: 


d¥c. 
Dy 5 - K,C,C, = 0 (6) 
ax 
d2c 
D, FT ~ Anpk Cap = 0 (7) 
ax 
dc, 
x (O-, Beate te ene (22) 
oe aaa) Cn = (0), Ca = Cao : 


(cf. Section 3.1 Background for Mass Transfer Modelling, 


Equation (23)). 
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In case of free cyanide, the differential equations for 


the reaction scheme developed previously can be written 


as 


saa 

Da oar - K CAC, ~ KyCn - KC,Co = ( (59) 
Ge 

D, eazy _ Gapk Canon = 0 . (60) 


A transformation technique is used to convert these equa- 


tions into a simpler form 


* = 
Cr = Con (61) 
Ge Maaci*+ Cc 
B > CUB 2 (62) 
where Cy and Cc. are arbitrary constants, and 
oak 
ag 
B * A* 
es = 3 
Dy 5 K -C,C, 0 (63) 
dx 
Pugs BLCR GR fon CB on i 
ax 
K K 
B 
se ef os —— + 0 ; (65) 
2° Kp Kp, C 
From Equation (60) 
= (66) 
C1 Cy (Cp ss C.) BxCe 
C 
2 <i ae ; (67) 
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Strictly speaking, Ci is not a constant because it is a 
function of C. which varies in the film. However, Ci can 
be assumed to be constant if the reaction is in the transi- 
tion regime where the concentration of B (cyanide) at the 
film is nearly the same as that at the surface. Also, a 


flat cyanate concentration profile is assumed. The solution 


for the transformed equations is 


ES - Ey ES = Ey 
Ey = EEN EKA tanh Sr (68) 


where 
ih K,/Ko 
pie: ehe 
O, 
EB, = (1+ 3) = 1+ 
TaB-A~As As 
* 
DrKeCay 
M = 5 
K 
0 
K K 
*k D B 
Bo Bo £ Ke CoO 
ere CE0°as 
As Ky Ka 
Cots +t ex C 
Bo P P Co 


ils (85) - ‘ 
It has been indicated by J. Rowley that Da = Dai and 
the stoichiometric ratio of ozone to cyanide, dap’ is known 


to be 1 for free cyanide. 
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3.4 Modelling of the Perricyanide/0,/UV Reaction 
3.4.1 Flat CNO Concentration Profile 


The approach is essentially the same as that of the 
free cyanide, so the derivation developed earlier (3.3 
Modelling of the CN_/0,/UV Reaction) for free cyanide 
can also be applied to ferricyanide with only minor 
modifications. Ke in free cyanide case is replaced by 
its counterpart K, the rate constant of the ozone/ferri- 


cyanide/UV reaction, and K the rate constant of ozone/ 


B! 
cyanate reaction,is assumed to be hy SMk fae K, determined 
from cyanide decomposition measurements. From Equation (1) 
Ky has the value of 8.5~x 10° aes at pH = 7 and 700[0OH] 
Q@iePHl2 as, calculated from the. Czapski, correlation, 

C. in the equations stands for the total cyanide content 

in the solution, and the bulk concentration of ozone, Cro! 
is equal to zero. Since the reaction between ozone and 
ferricyanide is rather slow, the transition from slow to 
fast reaction regime is tried (cf. 3.3.2 Transition 


Regime with Cn Equal to Zero - Film Model). 


3.4.2 Quadratic CNO Concentration Profile - Film Model 
While the assumption of flat CN concentration profile 

holds for both cases, the approximation of constant CNO 

concentration is questionable in the case of ferricyanide 


because the reaction between ozone and ferricyanide is 


considerably so much slower such that cyanate plays a more 
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Sioniticant role in determining the overall ozone mass 
transfer; even a slight difference in cyanate concentra- 
tion profiles may change the value of K completely. In 
order to improve the proposed model, the possibility of 
non-flat CNO profiles has to be explored. 

Considering the mass balances of ozone and cyanate, 


the following ordinary differential equations can be 


written: 
bac 
D = e 
A we KC,CRo + 3.25 KnC,Ca + KCa (69) 
9°, 
Da ay =— KCC, + KnCaCc 2 (70) 


The mass balance of cyanide need not to be considered 
because it iS in excess compared with ozone. There are 
a number of routes to solve these equations, for instance, 


(56/57/58) and collocation technique. 


numerical analysis 
The latter does not require as much computer effort and 
can give a result as good as the first one. In addition to 
its simplicity, the collocation technique can also furnish 
a better physical picture. 

In the collocation technique, a quadratic CNO con- 
centration profile is assumed. The three parameters 


involved i.e. by, bos b, in Equation (71) can be solved 


easily by imposing the boundary conditions onto the 


equation 
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i +b, 5 +b, 5 (71) 


where bir bo, b are the dimensionless parameters, and C 


Co 
the bulk concentration of CNO . The boundary conditions 
are: 

7 = = 
) x 6 Co = Cay C72) 
aC, 
a) x 37,0 “S Do aaa ee = FilCao : tyes) 


The amount of cyanate diffusing out of the film into the 
bulk phase per unit time is equal to the flow rate of 


cyanate going out of the reactor. 


dc, 
yet) x = 0 Asa = 0 Co = Cas (74) 
x=0 
where C is the surface cyanate concentration. The 


Cs 
boundary conditions and the concentration profile can be 


represented schematically in the following diagram. 


45 


f 
bd 


+ 
. 
4 
C201 
ry P . se7SO > oe f » «2 
" Jd 2M oe 4 we gap = @ 
oO.) - 
—. i 
(foitibnoo yasbriod sat... 
4 
: ~ =. - 
~~ " 
oS wi Qf 
’ 


7 : re. To 3! onienmtzrop os 6n &YD io ane 
jJes weld. she cs {6099 os ane ting _ _— 


wenden Sit Bo s0e enked 9 


ry) = xy ‘| = Pas - ae P Te 
at en! *P a 


f J 
/ 


ve As 
sy jolteitasesds efeasys soatane edd ea 


aa 7 , 
wy 2 ‘ey ¢ ie 
ad asp Sslitoug soptewrassnes. edt Ske enobtebowe ' 46 Batt 


oe S 
ms ais: 


b 


eapsrd oniwoalkhoeit sade go. viles it eMetioe, boaae 


From the above boundary conditions, 


can be obtained: 
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Solving Equations (38), (39), (40), 
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the following equations 
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(76) 


(77) 


(78) 


(79) 


This is obvious because the cyanate concentration in 


the film from which cyanate is generated has to be higher 


than that in the bulk solution. 
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S in the experiments, C 


Cs 


can be related to C 


Using the typical values 
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by 


Ca. = 2.34 Coo : (80) 


Cao can be obtained from the cyanide concentration by 


Equation (52) and it is only a function of K. 

A single ordinary differential equation containing 
Car Ca and x aS variables can be got by combining two 
differential equations in Equation (69) 


d’c d7c 


A — 
mn 5 3.25 D, 7 = 4.25 KC, OC 


+ . 
A Span 


é =— im 
(AM25eKC. ot RK pC eas 254 cae Ca, (81) 


The only variable is Car so Equation (81) is simply a 
second order homogeneous ordinary differential equation 


of which solution can be obtained readily. 


Cr = C, (complimentary) + C, (particular) 
= C,exp (/b, x) + C,exp (-vb, x) + b,/b, 
= Cc, sinh (vb, x) + C,cosh (vb, x) + b,/b, (82) 
Che C, or Cyr C. are the integration constants to be speci- 


fied by boundary conditions, be and bi are constants defined 


as 
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Substituting these conditions into Equation (82) 


Cas = Cc. + b./b, 
0 = C, sinh (vb, 6) + C,cosh(vb, 6) + b./b, 


Solving Ci, Co 


values into Equation (82), 


a. Sha b./b 
c, = [48 —2_*) sinh [/b, (6-x) ] + ———+*— 
i sinh (vB, 6) sinh (/b, 6) 
[sinh (vb, 6) - sinh (vb, x)] (84) 
dc Day bs b b aw 
> jeedoyy yeaa) patel uae EEE {4 [cme - 521cosh (vB;6) | 
L Cag 9% x=0 c.sinh(¥b,8) \°4 4 
As 4 
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The assumption of zero bulk ozone concentration is 
always valid as long as the reaction is in the transition 
regime. In fact, ozone bulk concentration drops fairly 
rapidly as the reaction proceeds from the slow reaction 
to the transition regime as indicated by the value of Cro 
calculated in 3.3.1 Transition from Slow to Fast Reaction 


Regime (not listed in the thesis). Equation (87) is a 


combination of Equations (85) and (86) 


SD,Vb, b. bd, 
RT = ————___ es [Cg - 2cosh (v5, 6) } (87) 
sinh (vb, 6) 4 4 


where by and by are defined by Equation (83) and 6=D,/Ko- 
3.4.3 Penetration Model 
Assuming flat concentration profiles for both 


CN and CNO , the unsteady state ozone mass balance can be 


expressed as (cf. Equation (4)) 


dG dC 


£ A 
D = pe + KC,Cao fi 3 yi 2 KCC, + Ky 2 (88) 


The solution to this partial differential equation is 


shown in Equation (28) 


RT 1 exp (-K, @) 
Ky = <+ /Kp,{ertvE 511+ yds] + | (89) 
Cas 1 TK, 9 
where 
K = KC aie is As) 8 tell © + K 
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The penetration model differs from the film model only by 


the expression for Ko. 


50 


* ~ 
x Sig-yqias | 5: ma test iy 


a 
ah 7 
Wist oft moct eretTib febom agitentense ett 
i ~*% 
A tot aobesd 1a 2 al 4 


i 


- 


“a 


CHAPTER 4 
EXPERIMENTAL 


4.1 Equipment Description 
4.1.1 Reactor and UV Lamp Assembly 

The reactor is shown schematically in Figure 1. The 
reactor is a QVF glass section, 18.6 cm in height and 
10.84 cm in diameter. It is sealed by silicone rubber 
gaskets between a stainless steel plate at the top and a 
quartz plate at the bottom. A stainless steel ring of 
inner diameter 10.3 cm (4 in) is used to hold the quartz 
plate (12.7 cm diameter, 7 mm thickness) in place by 6 
bolted rods from the stainless steel plate to the ring at 
the bottom. Another silicone rubber gasket is inserted 
between the steel ring and the quartz plate. The bolts 
have to be tightened gently and evenly with extreme care. 
Excessive torque on any one of the bolts would create 
strain and probably crack the quartz plate. The quartz 
plate is used because it is transparent to UV. 

The stirrer is also constructed of stainless steel. 
The end of the shaft (1 cm diameter) is slotted into a 
Teflon bearing (1.6 cmx 2.54 cm diameter) glued at the 
center of the quartz plate and three sets of stirring 
blades are attached to the shaft. An upper propeller is 
used to provide uniform mixing in the gas phase; a stirrer 


located slightly above the liquid surface is used to 
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minimize the gas phase resistance, and the bottom stirring 
blades (7.3 cmx 0.95 cm) are used to agitate the liquid. 
The stirring speed can be set precisely by changing the 
gear ratio. A stainless steel baffle of dimension 

8.9 cmx 4.13 cm is used to eliminate ripples and vortex 
formation in the liquid phase. It is possible to have a 
well mixed liquid phase and yet maintain a quiescent inter- 
facial area in the experiment, thus the interfacial area 
can be approximated by the geometrical surface area. The 
internal diameter at the liquid surface is 10.392 cm as 
measured by a Vernier scale, so the surface area is equal 
to 84 om? after accounting for the stirring shaft. 

Two pieces of stainless steel tubing of 4.76 mm i.d. 
are glued into the two holes drilled on the quartz 
plate as the inlet and outlet for the solution. Likewise, 
two tubes are screwed into the top plate for ozone. It 
is crucial to maintain a reproducible liquid level and 
surface area in the experiments, thus a mark is placed 
on the side of the reactor to indicate the proper liquid 
level. 

A 200 watt Hanovia medium pressure Mercury lamp 
(model 65410) is mounted on an aluminum reflector of 
dimension 44.5 cm lengthx 10.16 cm height x 5.08 cm width. 
The reactor is seated 2 cm above the reflector. To avoid 
stray UV light, the whole reflector is enclosed by metal 


except for a 11.43 cmx 5.08 cm space via which UV light 
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is directed through the quartz bottom plate of the reactor. 
One end of the reflector is open while the other end is 
slightly widened to house an air fan used to dissipate the 
heat generated by the UV.lamp. The UV lamp is hooked to a 
Hanovia UV ballast transformer. The whole reflector assem- 
bly is mounted firmly on a pole, and the position and orienta- 
tion are carefully marked so that it can be placed in the 
Same position for each run. Screens, placed between the 
reactor and the reflector, are used to vary the intensity 
of the UV radiation. 

In addition to the forced turbulent mixing, the dis- 
persion of ozone is further augmented by the upper propel- 
ler and the middle interface stirring blade. The latter 


can also reduce the gas phase resistance in the reactor. 


4.1.2 General Description of the Set-up and Flow Diagram 
The reactor is a CSTR with continuous flows of 0, 

and cyanide solution. The UV lamp is directly beneath the 

reactor. The set-up is shown schematically in Figure 2, 


and the following instruments and symbols are used: 


Ozonator: PCI Ozone Corp., model C2P-9C4 
Stanford, Conn., U.S.A. 
Capacity 9 gO,/hr from dry air. 
Gas-dryer: Drierite CaSO, 
W.A. Hammond Drierite Co. 


Xenia, Ohio, U.S.A. 
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Whitey, 3 way ball valves, 316 stainless 


satu Weferehin 

Whitey, on-off valves, 316 stainless steel 
body. 

NUPRO, ‘S' series, 216 stainless steel. 


Manometer. 


of various components of the set-up are des- 


cribed as follows: 


Drierite : 


Cooling water: 


M1 : 


M2 : 


Needle valve: 


Back pressure 


Dampers: 


absorbs the moisture in the air that otherwise 
May cause damage in the ozonator. 

removes the heat generated in the bearing. 
gives the pressure at the gas inlet. 

gives the reactor pressure. 

is located at the gas outlet to control the 
reactor pressure in such a way that the 
pressures in the reactor should be the same 

in the run and the outlet gas sampling, so that 
the steady state will not be disturbed by 
sampling. 

valve: holds the liquid in the reactor so 

that the liquid will not drain out by gravity 
or by pressure. 

serve to smooth the pulse of circulating 
liquid resulting from the positive displace- 


ment pumps. 


or 
: tatate @I1£ .asvisv iLbad 
ny 
— - ~ © eo ate f * . , q 
suintese Bit ,asvisv +2o=n0 
ay Soc 
a 
roe 
[uss trase GIS ,asize® $2"  ,OneDM = 
2 
Tas enionsM. -: : 
: 7 ws 7 7 
«+ eye Uu-cee sft to astvforsures avotsayv To Boks te 
* 
tawollied 2s 
fan? sits ails ot svuvdgerom sis @dioeas + ost 
Ospeeso Of af Spsmeb eauasd yer Je 
; = a Vp, 
ecciee es En Fe end ae Be bade the : >t Cee 
ae i ar ry ‘stip XJ pon Sis BSYOMNG SL » wong, pak Looe 
a+ | 
luezsig oft asvip : Ue 
(esetgq. 1otonet gid eeviie |: ’ 
7 ealis ' 
oll » 
aq> Lottiios ‘si3vo asd sax 26 Keteool ef :ovigv « { he 
~ > - 
i Seth “ate oa a2 - a : “eX ot 
ofi< 7 e6nS Y , 5B OOS o Wess ag Tos VSe52 - 
62 edit ad Sitede wos9sst offs nk SHuttacewy 
360.7 


, - tas ; i 
Qa ioto8ox sft mi bewpkl add shiod vewiw omens 
« pas aaa 


Si Winiga 


Ytivexy yd joo ciatb Oniétw biypht edt geek 


7 one ar V 
‘ ame S é ae yO 
A , are Pi rf 


Bubblers with fritted disc gas disperser: contains 4% KI 
solution for ozone sampling. 

Wet test meter: records the inert air flow rate after ozone 
is stripped by the bubblers. 

CO stops the ozone supply by connecting to 
"VENT' position. 

lg samples inlet 0, by turning to 'BUBB' 


position from 'REAC'. 


tbs Be samples outlet 03 by turning to 'BUBB' 
from 'NEED'. 
T4: samples outlet liquid by turning to 


"SAMP' from ‘EXIT’. 
Dis closes to stop 0, to the bubblers. 
D2: opens to purge the 0, trapped in the tub- 


ings by air. 


Polyflo tubing deteriorates slowly in the presence of 0, 
and must be replaced periodically. A detailed operating 


procedure is given in Appendix 2. 


a2c4.5 Scatety Measures 


Ozone is a poisonous gas, SO care should be taken to 
avoid ozone escaping into the working area. The ozonation 
experiment has to be done inside a well-ventilated fume 
hood. If the characteristic odor of ozone is sensed, the 
ozonator should be shut down immediately until the faulty 


part is fixed. Cyanide should also be handled with 
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extreme care; gloves should be worn to protect hands from 
@irect contact with cyanide or its solution. Cyanide 
solutions must be kept at a pH higher than 12 to avoid 
formation of the volatile compound HCN. After the 

cyanide solution is drained, sufficient water should be 
used to rinse the sink as well as to dilute the cyanide 
solution. It is a good practice to immerse the ozone 
outlet into the waste cyanide solution so that the cyanide 
level is lowered by ozone oxidization before it is dis- 


charged. Hands should be washed after experiment. 


4.2 Reactor Characterization 
4.2.1 Mixing Pattern 

Since the analyses of results depend on the assump- 
tion of ideal mixing in the continuous flow stirred tank 
reactor, this assumption has to be justified before 
further work can be done. Should any non-ideal mixing 
pattern be detected, the reactor has to be redesigned to 
facilitate more efficient mixing until the requirement is 
met. Ideal mixing refers to the case of complete mixing 
so that the properties of the reaction mixture are uniform 
in all parts of the vessel and are the same as those in 
the exit stream. In other words, there would be no con- 
centration or temperature gradient in the reactor. 

The mixing behavior of a CSTR can be studied by a 
residence time distribution RTD experiment $771 59) 


which a 'step function' input is imposed on the reactor, 
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and the response in the exit stream is recorded. The resi- 
dence time ea is the average time the molecules spend in 
the reactor. In this experiment, diluted HCl was used 
to provide the step change, and a pH meter and glass 
electrode were used to monitor the exit pH. Details are 
given in Appendix 3. 

A typical recorder output for a stirring rate of 
85 rpm is shown in Figure Al. The time delay ty is indi- 
cated in the time axis by a sudden drop in pH, and equal 
to 2.48 min as measured directly from the recorder output. 


At time larger than t the pH drops down exponentially 


L’ 
and@ethen Vlevels¥of f° afteres minutes. mhLogiCe -i@)7(C-.j}- Co) 
at different stirring rates are plotted against t in 
Figures A2-A5. For the first three graphs (Figures A2-A5) 
forwhich the stirring rates are greater than 0 rpm, 
straight lines can be drawn through the points. th can 

be measured readily from the graphs as the time at which 
the lines bend sharply from a zero slope to some negative 
values. The values of th for, all*fours cases: including sce 
one without stirring are 2.48 minutes. The mean residence 
times t. are calculated correspondingly from the slopes of 


R 
the lines and found to be 7.84, 7.87 and 7.46 minutes for 


Miepeos and 59.5 rpm. 
The results are tabulated in Table l. The values of 
t, and cn from graphical methods (Figures A2, A3, A4, A5 in 


Appendix 3) are compared with those by calculation 


(Equations (A22), (A23) in Appendix 3), and good agreement 
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RPM Graphical method CSTR 


th min tp: min assumption 
102 Zed 7.84 good 
85 Lao Firsts e good 
i ee ig 7.46 a 
0 4.82 6.07 no 
Liquid flow rate, E. 71 ml/min 
Reactor volume, V ee tai | 
Inlet void volume (from the feed tank to the reactor) 66 ml 


Outlet void volume (from the reactor to the glass electrode) 110 ml 


Calculated time delay, ty =(66 + 110)/71 = 2.48 min 


Calculated average residence time, = = 533/71 = 7.51min 


Table 1. Comparison of t. 


calculation at different stirring rates. 
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can be seen. Finally, it can be concluded that the mixing 
is perfectly ideal with stirring rate down to 70 rpm; a 
slight deviation from ideal mixing begins to show up at 
59.5 rpm which causes a decline of an from 7.87 to 7.46 
minutes. The assumption of ideal CSTR that has been made 


earlier is valid at a stirring rate higher than 70 rpm. 


4.2.2 Physical Mass Transfer 

The physical mass transfer coefficient for 0. in water, 
Kor can be measured either directly or indirectly. 
The first method involves measurements of the partial 
pressure and the concentration of the gas component in the 
gas and solution phase, while the latter the correlation 


of K, of the gas component to that of CO, by the assumption 


0 2 
of equal contact time 6 in the penetration model. Being 


able to yield K, without making use of any model, the direct 


0 
method is deemed to be better than the indirect one. How- 
ever, the indirect (H,0-CO,) method should be as good as the 
direct one if the contact time @ only depends on the hydro- 
dynamic conditions such as the detailed design of the 
reactor and the way the solution is stirred, but not on the 
nature of the gas molecules itself. CO.-H,0 method has 
received wide-spread acceptance, and finds use in 

cases where the gas does not dissolve enough in water to 
give a detectable concentration in water, or no suitable 


analytical method exists, or the gas reacts instantaneously 


with water so that Ko cannot be obtained. co. rts) 
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particularly used because there is no gas phase resistance 
nor any known fast chemical reactions which may enhance 
the transfer. 

The 0,-H.0 method is quite straightforward. Ores 
in contact with distilled water in the continuous flow 
well stirred reactor shown in Figures 1 and 2 and the ozone 
concentration in the outlet water was measured. O.was mea- 


es 


sured by iodometric method while CO, wasby acid base 


2 
Litracion. The detailed operating procedure, theory and 
calculation are shown in Appendix 4. 

The Ky from indirect (CO. 5 


measurements are plotted in Figure 3 as a function of stirr- 


-H,0) and direct (O,-H O) 
ing speed. Two straight lines can be seen,with the indirect 
value 3.8% lower. The close agreement between these two 


values validates the CO, correlation approach, which is 


2 
indeed a powerful tool whenever the direct measurement 
method is not feasible. The slopes of these two lines are 


(ol) The 


both 0.744 as compared with 0.8 from Sata's work 
physical mass transfer coefficient of ozone is correlated 


with the stirring rate by the following expression: 


log Ky meal) 2 4 0.) 441.00) rpm oi. (90) 


The values of Ko by the two methods are also compared 
in Table 2. It can be observed that Ko depends strongly on 
the degree of agitation. The increase of Kp with stirring 
rate is comprehensible, because the surface renewal rate is 


enhanced in the penetration model point of view, or the film 
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Figure 3. Plots of log K,S by both 03-H,0 and CO5-H.,O 
methods vs log rpm at 22°C. 
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Table 2. Comparison of K 37H, 
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thickness decreases in terms of the film model. Besides the 
stirring rate, the physical mass transfer coefficients 

were measured at different stirring blade positions; a gradual 
decrease in transfer coefficient from 10.69 cm/sec with 

the blades 1 mm below the surface, to 7.45 cm/sec at the 
surface and 1.998 cm/sec just slightly above the surface 
was observed. This phenomenon proves that the proper 
Maintenance of surface level is not trivial at all. The 
increase in Ky with the blades touching the solution is 

due to the increased mixing at the interface as well as the 
increased area caused by the ripple and vortex. 


It is known that CO, may react with water chemically 


2 


according to the following consecutive reversible reactions 


+ 
GOmt— Wie: HCO; 


CO. + H.O => H 


2 * 2 


A question might rise concerning whether the Ky (CO.) 


measured is entirely due to physical absorption or a 
combination of physical and chemical processes. The CO, 
system had been studied thoroughly by panckwerts©2) ; athe 
was found that these reactions were relatively slow com- 
pared with the physical absorption process. Technically, 
the system is classified in a ‘slow reaction regime' where 
the mass transfer coefficient Kk is the same as K). 

The low partial pressure (2 x Thay atm at most) and 
solubility of ozone create problems in analysis. Even when 


a micro burette is used, serious scattering is observed at 


low stirring rates because cf the low solubility. There is 
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no interference other than the possible decomposition of O,- 
However, Since the background pH of the distilled water is 
5, the effect of decomposition is deemed to be negligible 
(cf. Equations (1) and (2)). The assumption of zero gas 
phase resistance is justified by the value of the gas phase 


resistance measured by J. Rael eae 


4.2.3 Characterization of UV Intensity 


In order to study the UV effect on the O.,-Ferricyanide 


3 
reaction, the effective UV intensity, E, in the reactor for 
different screens had to be measured. A sensitive chemical 
actinometer (potassium ferrioxalate) which produced an 
eaSily analyzable product (Fea) upon the absorption of UV 
was used for UV measurement. The Peau produced was analyzed 
by the absorbance of the red Fe*t-1,10-phenanthroline 
complex. The details are given in Appendix 5. 

The concentrations of ee and the corresponding 
absorbances are shown in Table A3 in Appendix 5. A calibra- 
tion curve is also constructed by plotting the absorbance 
against the concentration of Rea and is given in Figure A6 
MmmeAppendix. 5... This,is,a straight line passing through the 
origin with a slope equal to e&, the product of molar 
absorptivity and the light path. From Figure A6, the term 


€2£,the product of molar absorptivity,and path length is found 
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2+ 
tion, the number of moles of Fe produced, n,, the number of 
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quanta absorbed by the ferrous complex per unit time, Ios 
the equivalent wattage, E, and the relative intensity are 
tabulated in Table 3. E was also measured directly by a 
Black Ray UV intensity meter. The results are given in 
Table 4. For this method, E is calculated by multiplying 
the reading on the meter by the area of the opening, which 
is measured to be 11.4 cmx 5 cm = 57 cm”, Comparing the 
results from these two methods, the values of E from meter 
measurement are lower owing to the fact that the meter 
cannot respond to wavelengths larger than 2700°A. The 
values of E by chemical method give the total incident 
energy with wavelengths smaller than 5400°A, whereas the 
readings from the UV meter give only the fraction of energy 
with wavelengths not greater than 2700°A. In spite of the 
difference in absorptivity, the relative intensities based 
on full UV (as 100%) are expected to be the same for both 
methods for a given light source; however, they differ 
substantially as revealed in Tables 3 and 4. The discre- 
pancy comes from the uncertainty in the reading of scale 
from the UV meter at full UV intensity, because it is just 
at thesedge of the scale. It is confirmed by comparing the 
ratios of E(fine) to E(coarse), and they are found to be 
2.6 and 2.4 for the chemical method and meter method 
respectively, so the error must originate from the measure- 


ment in the ‘full UV' case. 


The possible error sources for the chemical method 
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i) the undesired photo-decomposition of ferrioxalate 
after it is taken out from the reactor. The 
dim light used in performing experiment and the light 
source in Spectronic 20 may cause further reaction of 
ferrioxalate. 

iv) the assumption of no contribution to eee production 


for wavelengths greater than 5400°A. For more accur- 


ate work, the values of dp has to be known continuously 


for all wavelengths, or a monochromatic light source 
must be used. 

iii) the error in y due to the uncertainty in the energy 
distribution function. To.obtain a.better result 
requires the knowledge of continuous energy distribu- 
tion function or the use of monochromatic light. 

In spite of the uncertainty in the chemical method, 
the values of the relative intensities measured using this 
method are thought to be more reliable and are used 


for the subsequent experiments. 


4.3 Analytical Procedures 
4.3.1. Ozone Analysis 


Ozone was measured by a well-documented iodometric 
methoa ‘’7? . In this method, the gas was bubbled through a 


neutral 4% KI solution and ozone reacted to liberate iodine 


H,O0 + 31 WROTE o> OF eh al eae OH . 


Under neutral or alkaline conditions, the interference of 


0. is minimal. The neutral I, solution (e.g. in bubbler Bl) 
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was then acidified by 30 ml 1:10 H,SO, solution, and 
titrated immediately by 0.1N Na,S,0, solution which had 
been standardized against standard 0.1N I, solution 
(Fisher Scientific Co.).. Since the I. in the solution was 
subject to air oxidation after acidification, the titra- 
tion must be finished quickly. A magnetic stirrer was used 
to facilitate thorough mixing. Starch solution, as recom- 


mended by Taras et Ab eae, 


, waS not sensitive enough to 
give a noticeable change at the end point; the end point 
was indicated by the complete decolorization of the I, 
solution, which would usually be reached in 30 ml of the 
titrant added. It had a remarkable sensitivity of +1 drop. 

The ozone in the second bubbler, B2, was measured 
Similarly. The volume of titrant required was generally 
around 3 drops and seldom exceeded 10 drops. The blank 
was prepared by passing the same volume of air through the 
bubblers, and the blank solution was titrated for I,. The 
blank was approximately 0.05 ml, and was subtracted from 
the volumes measured above. After correcting for background, 
the sum of the titrant used should indicate the total 


amount of I and thus 03, in the bubblers. Detailed cal- 


oo 
culation and gas phase sampling are shown in Appendix 2. 


Gicndis 2 Analytical Method for Free Cyanide Determination 


The measurement of free cyanide is well documented in 


the Titec’ wt. Free cyanide can be isolated from the 


interfering matrix by simple distillation in an acid medium, 
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and then determined by titration, colorimetric measure- 


78 ' 
we) or cyanide selective electrode monitoring ‘/??, 


ment 
Titration is deemed to be the better method as far as con- 
venience and accuracy are concerned, thus it was used for 
this research project. The procedure for CN. titration 
has been described in full detail in reference (77). 


Standard 0.01 N AgNO. was prepared by dissolving 


3 
Oa Oo, OL AgNO. in a 2 liter volumetric flask, and the 
indicator solution by dissolving 0.02 g paradimethylamino- 
benzyl rhodanine in 100 ml acetone. The titration was 
Carried out at pH 11 or above, seven drops of indicator 
solution were added to the receiving flasks, the color of 
the solution was canary yellow and would change to salmon 
pink at the end point. The color change was not very sen- 
sitive, thus a white background helped to improve the sen- 
Sitivity significantly. 

The blank was prepared by diluting 50 ml of 2N 
NaOH solution to 250 ml and the same procedure as above 
was used to titrate the blank solution. The blank should 
be around 1 drop. 

This modified Liebig method, when used at a cyanide 
level above 1 ppm, has a coefficient of variation of 2% for 
distilled samples, and the sensitivity is approximately 0.1 
ppm of CN . However, the color change 1S NOtoatstlucteac 
this point. At 0.4 ppm, the coefficient of variation is 


four times that at a concentration level greater than 1 ppm. 
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The reaction between Ag* and CN’ can be shown as 
+ _ - 
hae 20N Ag (CN) . 


so the concentration of cyanide in solution is 


+ 
2x[Ag ] (volume of Ag™ used - blank) 


sample volume 4 Cr, 
4.3.3 Ferricyanide Analysis 
ou) ion exchange chromatography (81) , nor {82) and 


(77,83,84) are all considered to be 


cyanide distillation 
feasible for determining ferricyanide in solution. The 
first three methods require expensive instruments yet the 
accuracy is at most on the same order as that for distilla- 
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t10n The measurement of total cyanide content in 


ferricyanide by cyanide distillation is involved due to the 


#2 


stability and the inertness of the complex towards reactions. 


Ferricyanide can liberate HCN only under drastic conditions 
such as high acid concentration, reflux, presence of 
catalyst or strong complexing agent. A procedure of air 
circulation to remove HCN liberated from the boiling ferri- 
cyanide solution has been proposed by Gould, Afghan and 
Brooksbank ‘°>? who claimed that a sensitivity of 0.5 ug CN 
in 15 ml of sample and a recovery of 100% at CN concentra- 
tion above 1 ug/% could be obtained by titrating the 
liberated HCN with Ag’. However, the experiments were 
repeated using exactly the same set-up and procedure as 


recommended by the authors and the result as claimed could 
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not be reproduced. The percentages of recovery were spora- 
dic, and ranged from 40% to 80% even when air was bubbled 
through the solution for more than 2 hours. The unrecovered 
portion of CN was found to be left in the solution and it 
could be driven out by boiling. Even though the ferri- 
Cyanide had been broken down completely into HCN, the 
experiments with recirculating air failed because of the 
high solubility of HCN gas in cool water. A large percen- 
tage of HCN in the carrying air redissolved into the water 
condensate near the top of the distillation column and thus 
returned to the bulk solution. A combination of long bubbl- 
ing time and large air recirculating rate might drive the 
process to completion. 

It has been possible to develop a method that does 
give quantitative results by modifying the method of Gould 
et al. A schematic of the distillation apparatus used is 
shown in Figure 4. In this new method, the HCN formed by 
decomposition of ferricyanide was carried out of solution 
with steam, and dissolved in the condensate which was 
then mixed with 2N NaOH solution in the receiving flask. 

An appropriate volume of the condensate should be collected, 
too large a volume might cause an insensitive indicator 
color change at the end point while an insufficient volume 
might lead to a low recovery. The volume to be collected 
could either be determined by a CN probe to monitor the CN 
level, or just by trial and error. A volume of 150 to 200 


ml was recommended. 
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1) Reagents 


Sulfuric acid conc., 'Baker Analyzed' Reagent, BDH 
Sodium hydroxide, Fisher Scientific Cromer 

Silver nitrate, Analab 

Potassium cyanide, 'Baker Analyzed" Reagent, BDH 
Potassium ferricyanide, 'Baker Analyzed' Reagent, BDH 
Magnesium chloride-6H.0, Fisher Scientific’ Co. ; 
Mercuric chloride, 'Baker Analyzed' Reagent, BDH 
5-(P-Dimethylaminobenzylidene) rhodanine, 


Eastman Organic Chemicals. 


a) Apparatus 


Any common open distillation set-up could be used for 
complex cyanide distillation. Being a distillation device 
for HCN, it also served as a reactor in which the breaking 
down of the complex cyanide took place (Figure 4). The 
joints used in the apparatus were all of the size 24/40; a 
3 neck round bottom flask with 2 stoppers was used to 
contain the solution to be distilled. A 380 watt Canlab 
heating mantle acted as a heating source. A magnetic 
stirrer provided adequate agitation during the run. The 
flask was connected to a condenser through an elbow, and 
the condenser was placed at a slightly inclined angle so 
that the condensate together with the HCN could drain by 
gravity into the receiver containing 2N NaOH solution. 

The tip of the condensate drainage tube was dipped into the 


NaOH solution to prevent HCN from escaping into the air. 
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An enlarged section was placed between the condenser and 
the rubber tubing leading to the receiver, so that the 
solution in the receiver would not be sucked back into the 
condenser in case of sudden cooling of the condenser and 
the flask. 

The presence of grease in the condensate might interfere 
seriously with the cyanide titration. It was recommended 
that grease should only be used when needed, and a small 
amount of grease should be enough to stop leaking. AS a 
matter of fact, the joints were only greased when they were 
new, and there was no sign of leaking thereafter due to 
insufficient greasing. 

The tip of the condensate delivery tube was dipped 
into the 2NNaOH solution (approximately 50 ml) in the 
receiver. The cooling water was turned on and all the 
joints were tightened firmly to avoid any leakage during 
the run. One stopper was pulled out, and a known volume 
of ferricyanide solution was pipetted into the flask. 
Sufficient distilled water was added to make up the volume 
to approximately 600 ml. Then 20 ml 6.8% HgCl., and 10 ml 
513 MgC1,-+6H.O solutions were poured into the flask and 
followed by 60 ml conc. H,SO,- The stopper was put back 
into place hastily after introducing sulfuric acid into the 
solution. 

The stirrer and the heating mantle were switched on. 


During the heating up of the solution, bubbles were coming 
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out from the tip of the delivery tube due to the expansion 

of air inside the system. The solution boiled in 15 minutes 
and the characteristic yellow color of ferricyanide dis- 
appeared after boiling for 5 minutes, which indicated that 
the ferricyanide complex was completely broken down. An hour 


of reflux as recommended by Gould et ales 


was not required. 
After 120 ml of condensate were collected, the receiver was 
tilted so that the tip of the condensate delivery tube was 
just above the solution. Another 100 ml of condensate were 
collected. the purpose of this step was to wash out any HCN 
remaining on the wall of the condenser and the tube by the 
condensate. Then the receiver was removed and another 
receiver containing fresh NaOH solution was immediately put 
into place to collect further condensate to ensure that 
there was no cyanide left in the system. The distillation 
was stopped after 100 ml of condensate had been collected in 
the second receiver. 


The contents in the receivers were titrated for CN 


by standard silver nitrate solution as mentioned in 


is eae Analytical Method for Free Cyanide Determination. 


iii) Evaluation of Method 


ce the volume of ferricyanide sample added, ml 

vi the volume of titrant for the first receiver, ml 
V5 the volume of titrant for the second receiver, ml 
Vy the blank, ml 


[CN] the concentration of cyanide in the stock ferricyanide 
fe) 


solution, M 
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The reaction between Ag* and CN can be shown as 


+ -_ - 
Ag’ + 2CN) ==> Ag(CN), 


<— 


so the amount of CN. recovered is 


NV eer 


+ 
l ae 2V,,) [Ag ] 


The percentage of recovery R is 


_ amount of CN recovered ‘ 
R = amount of CN- added 500% 


‘a + 
2(V, + V, - 2V,) [Ag] 


i ee ee nn OPM OE) (92) 
V_ [CN] 
a O 


The results are tabulated in Table 5 and the amount of 
CN added (column 6, Table 5) is plotted against the amount 
measured (column 5, Table 5) in Figure 5. It is found that 
in the range of 3.155 mmole to 0.2984mmole CN, the aver- 
age percentage recoveryis 98.7%, which can be regarded as 
quantitative. 

At a cyanide level lower than 0.06231 mmole, only a 
portion of cyanide can be recovered as revealed by the down- 
ward deviation from the ideal line in Figure 5. The yield 
drops from 95.9% at 0.06231mmole of CN to 90.5% at 0.0249 
mmole, and would be expected to drop further at reduced 
cyanide contents. The loss of cyanide at low concentration 
might be due to the hydrolysis of cyanide catalyzed by 
eee or Mg 7. 

It can be concluded that the proposed open distillation 


method for complex cyanide determination is considered to 


be valid at cyanide contents higher than 0.07 m mole in 
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CN 

[cn] . Vea. 

M ml ml 

0.1246 25 150.92 
10 59.28 
5 30.04 
2 hl eye Ba 

-3 * 

1.246x10 50 5.997 
20 ar A ay 

0.1493 25 182.51 
10 H2o22 
2 14.64 


[Az ] = 0.01012 


Vy = 0.0 


* Titrated 


Table 5. 


ferricyanide by the open distillation method. 


by micro-burette 


CN 
measured 
m mole 


061 
0239 
-6100 
2443 


OOF Ww 


0.05976 
0.02256 


3.697 
1.468 
0.2987 


RZ 


Percentage of cyanide recovery at various levels of 


79 


ji 2 LV a 


HoOoI SH 


ar 
Pn 
¢ 
borite Sar 
iiom Ae 


Leos He .O 
vit : ¢ | 
@ OL. 
PAT ‘ +f 1) 


\ 4 cA Heh) 
20. ae 
+5 i oe 

Hy. 


16 ¥csVase7 Sblasayp Po — st 


‘notmniltiath gage 


= Saribd-e tole rad asanset 


Ws 


OES || Ee 

‘GE, RR ot 
AO, OE . 
e,d1 ¢ 


i aie 
“ me ; 


BS 
8. 7 _ 


rina’ ~ 


ea2 oar . ‘ 


Cyanide Added (m mole) 


0.01 0.03 0.06 011 0.3 06 1.0 3.0 6.0 
Cyanide Measured (m mole) 


Figure 5. Plot of cyanide content measured vs cyanide 
content added in the determination of ferri- 


cyanide at 22°C. 
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700 ml solution, beyond which an accurate calibration curve 


is required. 


4.4 Experimental Procedure - CN /0,/UV Reaction 


The experimental procedure is described fully in 
Appendix 2, General Operating Procedure. The inlet and 
outlet cyanide concentrations, the gas phase ozone concen- 
trations in the inlet and outlet streams, pH, and the gas 
and the liquid flow rates were measured. The cyanide 
feed solutions were prepared by dissolving 0.6-5 gm of KCN 
powder ('Baker Analyzed‘ Reagent, J.T. Baker Chemical Co.) 
and 10.00 gm NaOH (Fisher Scientific Co.) in 9.435 liters 
of distilled water. The pH was measured by a Fisher 
Accumet Digital pH meter. The experiments were run at 
108.6 rpm, 0 and 100% UV. The gas phase analysis is des- 


cribed in Appendix 2. Cyanide was measured by silver 


fverate ¢icpatione (ef244iGa2 Analytical Method. for Free 
Cyanide Determination). 


The experiments were carried out at a liquid teed 


rate of 0.765 ml/sec, a cyanide concentrations of 


3 3 3 3 


eye 1. B4x 10°, 3.701 x10  , 0.6113 ~ TOske 


7.694x 1072 M, and UV intensities of 0 and 100%. Then 


it was repeated at the feed rate of 1.094 ml/sec and five 


: = 
different cyanide concentrations from 1.140 x 10 to 


5 


7.880x 10 - M without UV. 
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4.5 Experimental Procedure - Ferricyanide/O,/UV Reaction 

The detailed operating procedure is given in Appendix 
2, General Procedure. The ferricyanide concentrations of 
the feed and the outlet stream, the gas phase ozone inlet 
and outlet partial pressures, the liquid flow rate, the 
temperature and the pH were measured. The ferricyanide feed 
solutions were prepared by dissolving 0.2 to 40 grams of 
potassium ferricyanide solid K3Fe(CN) (, "Baker Analyzed' 
Reagent, BDH, in 19.95 liters of distilled water. Sodium 
hydroxide or phosphate buffer was used to keep the pH of 
the solutions at 12 and 7 respectively according to the 


scheme shown below. 


pH Ferricyanide Chemical used 
ey concentrations below 500 ppm Fisher Scientific, 
19.95 2 of solution PHe/epuULce” 
concentrate, 100 ml. 
7 concentrations above 500 ppm buffer concentrate, 
19.95 & of solution 300 ml. 
12 all concentrations sodium hydroxide 
S2435..%. of solution TORG. 


The pH of the inlet and outlet solutions were measured by a 
Fisher Accumet Digital pH meter. The experiments were run 
at a stirring rate of 108.6 rpm and 4 different UV inten- 

sities ranging from 03 (UV lamp orf), to 100% (O.s19pwattcs)), 


with the ozonator setting at 100%. 
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CHAPTER 5 
DATA TREATMENT 


The reactor system has been modelled to relate the 


Co’ Ko and K to the measura- 


ble quantities such as cyanide depletion, temperature and 


unknown variables e.g. RT, C 


ozone concentration. The ozone usage, RT, can be obtained 
from a gas phase mass balance or it can be obtained by the 
simultaneous solution of the equations that relate RT to 

the amount of cyanide oxidized and the reaction rate cons- 
tant. These model equations depend on the reaction regime, 


mechanism and stoichiometry assumed. 


5.1 Summarization of the Equations to be Solved for the 


Rate Constant 

For the transition from slow to fast reaction regime 
(cf. 3.3.1 Transition from Slow to Fast Reaction Regime - 
Film Model) and for the reaction scheme outlined in Section 


3.2, the following equations apply: 
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where the number 1.8389 x LOG? is the diffusivity of ozone 


at 20°C calculated by Equation (90), and the asterisk * 
indicates that the values of these variables may change 
with models and experimental conditions. 

There has been difficulty in obtaining a reliable 
gas-phase mass balance for ozone because the inlet ozone 
concentration does not differ significantly from the outlet 
one. Therefore, an average ozone consumption ratio Sor 
which is obtained by averaging the values for all runs, is 
used to smooth the values of RT in Equation (94). By using 
Equation (94), RT has been smoothed and the idea of the gas 
phase material balance is still preserved. Zero gas phase 
resistance is also assumed in the above derivation. 

Equations (47), (48) or its equivalent, which des- 
cribe the transfer of ozone across the interface for a 
defined regime, are obtained directly by solving the 
differential equations; no experimental work is involved. 
Equation (55) is to relate ozone consumption to the amount 
of cyanide oxidized, which can be measured accurately. RT 
can also be calculated by the gas phase ozone mass balance, 
Mee sBouations (56) and (94);"so™that there™are*threewri 's 
from three independent sources. If Cro is equal to zero, 
Equations (47) and (48) can be reduced to a Single equation 
(Equation (57)). The unknowns RT, K, and K can be. solved 


from these three RT equations (Equations (57) and (55) and 


one of (34), (94)). Li Cro iS snot equal to zero, there 
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would be one more unknown and One more RT equation (Equa- 
tions (47) and (48) instead of (57)), the unknowns can 
still be solved readily. However, one of the RT equations 


is not required if the assumption of 


ai 1) kK (34) 


is made. It is a logical decision to drop the RT obtained 
from the gas phase mass balance (Equations (56) and 

(94)) because it is the most uncertain one. Equation (52) 
is derived from the well-mixed model, and provides the 
relationship between cyanate concentration and the rate 
constants Kar Ky. Equations (96) and (1) are from the 
Mreererure (cr. 2.1.1 Ozone |\Chemistry and 2.2.1 Solubility 


of Ozone in Aqueous Media), while (97), (98), (99) are from 


direct measurement. Equations (100), (101) give the value 
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of diffusivity (cf. Appendix 4 Detailed Physical Mass Trans- 


fer Measurement) and Equation (102), the value of physical 


mass transfer coefficient (cf. Appendix 4) at various tem- 


peratures. 
The variables Pye dh Cai’ Cao! Fue V8 Sel On), Ky 
and D. are measurable. Inserting the values of these 


A 
variables into the key equations (Equations (47), (48), 


(52)>,22(55)-andone of (34).,- (94) and. (56), Ke can be solved 


by Mueller's agitate ncaa (cf. Appendix 1). 
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Bice CN/0,/UV Reaction 
5.2.1 Fast Reaction Regime 

The theory developed in 3.3 Modelling of the CN /0,/0V 
Reaction can be greatly simplified by assuming zero ozone 
concentration in the bulk phase, which is valid for fast 


reaction. Equation (25) in 3.3.3 Fast Reaction Regime can 


be used 

Br PAS (58) 
where 

Ki = KeCn, eg) Sh Pa KC G5 + Ky < 


Kk. can then be related to RT after correcting for the gas 


phase mass transfer resistance (kg,cm/sec) which is equal 


85) 


to 2.15 cm/sec at 22°C as measured by Rowley | uSing the 


same reactor configuration and operating conditions 
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The process variables Pyar ep Cais Cao’ Fue V, S and pH 


are measurable, and D gees Ce can be generated by 


ytd pe 
Equations (95) to (102) listed in 5.1 Summarization of the 


Equations to be Solved for the Rate Constant. 


Four different cases are tried in this regime 
fi section 5.1). 
Case A: RT is obtained directly from the gas phase mass 


balance 
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Ky and K, can be solved simultaneously by Equations (103), 


foe), ao2) and” (55). 


Case B: RT is calculated from R multiplied by the average 


measured ozone consumption ratio So 


Rot] he oUR, =24.5 BONC Cae! (94a) 


(85) obtained a So ValGelor wl ese cees aise 
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Case C: Rowley 


Rl rik whe FL Ca; ~nh Sve (94b) 


Case D: The gas phase mass balance is ignored in this case, 


a constant ratio of Kp/K, is assumed 


Kote = heel | € es (34) 


Bocetions. 103), (34), (52), and, (55) are used to solve K, 


and Ky. The key equations are Equation (103), 
K_H 
rat MORK PR 
RT = a -ARA (103) 
22 /DimaK 
RT ui 
where 
Ki = K,*CNF a Ky a en oer K, °CNO ° 


Pouations (>), (>2) and one of ‘the Equations (56) to (34), 
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RT = RO case A (56) 
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Ky and K, can be solved by the computer program THES using 
Equations (103), (55), (52) and one of the equations from 
(56) to (34) in addition to the basic equations from (95) 
to (102). The input variables are surface area ARA, inlet 
and outlet cyanide concentrations CNI, CNF, rate of ozone 


consumption n Panter ozone partial pressure P pH, tem- 


A ye / 


perature T and liquid flow rate Fo: 
The whole computer program THES as listed in Appendix 7 
consists of a main program and two subroutines, MULR and FMU. 
MULR is the program for Mueller's method, which requires 
the inputs of initial guess XLI and XRI, maximum number of 
iterations IEND, convergency testing number EPS and the other 
experimental numbers to be fed into FMU, for instance, diffu- 
sivity DT, outlet cyanide concentration CNF, physical mass 
transfer coefficient KO, interfacial ozone concentration 
(saturated) CS, amount of cyanide reacted per unit time R, 
measured ozone consumption RO, initial concentration CNI, 
ozone decomposition rate constant XKD and an indicator II, 
which varies from 1 to 4 denoting cases A to D. Equations 


Prot 52.)), (55) and (56)' to (34) are stored in the sub- 


routine FMR. 
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The general algorithm for cases A to C is to 
use Equation (56) for case A, (94a) for case B and (94b) 
for case C to generate the value of RT; Equation (52),the 
value of CNO;and Equation (104), a combination of Equations 


(52) and (55), the value of Kye 


ys eh taaci 
B BME 1 


. BRE t K-°CNF 
K-°CNF-(=-- 1)-K, Kg 


(104) 


1 
CNIG— CONE 


Once a tried Ke value is fed from the program MULR into FMU, 


first K, and then CNO are defined. Another comparison func- 


tion F is introduced in FMR where 


F =RT(from Equation(103))- RT(from Equations (56) or (94a) or (94b)). 


The value of F is brought back to MULR for comparison and 
generation of a new K, value which is then fed into FMU until 
F is less than the set limit, EPS. 

The algorithm for case D is less involved. The 
values K. and CNO are obtained from Equations (34) and (52) 


B 


respectively. 


F=RT(from Equation (103)) - RT(from Equation G5:5))) 3 


These equations are stored in FMU and the same trial and 


error method is used to solve for Kes 


ou 


ot gS Transition from Fast to Instantaneous Reaction Regime 
The same THES program can be applied to this regime 
with modifications. The key equations are Equations (34), 
C22 to oo and (sc) (cr. 3.3.45 Transition trom Fast "to 
Instantaneous Reaction Regime). M and E. are functions of K 


is 


Stumcno we wiitcn,) sin turn, tsearrunction oL Ky. 


F = RT(Equation (68)) -RT(Equation (55)) 


Two additional subroutines FNU and NULR are required in this 
regime because Equation (68) is an implicit function of EN, 
the enhancement factor. The initial bounds of Ky are set in 
MULR, which then directs a tried K- value to FMU to define 
E; and M. Then, EN is solved by the other subroutines FNU 
and NULR based on the same principle as that of FMU and MULR. 


The initial guesses of K, and EN are critical in this program, 


Ae 
as an improper choice of initial conditions will render a nega- 
tive value to the term M(E, - EN)/(E,- 1). The square root of 
a negative number is not defined, and thus the value of EN 
does not converge. Several sets of initial K- and EN 

values have to be tried in order to get the results converged 
for the cyanide concentrations studied; a unique general set 


of initial conditions applicable to the whole cyanide concen- 


tration range covered in the experiment does not exist. 
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Wes: Ferricyanide/0,/UV Reaction 
The reaction scheme and the mass transfer model have 
been developed in Sections 3.2.3 and 3.4. As it has been 
shown in the 03/CN- reaction that the gas phase mass trans- 
fer resistance only accounts for 2% of the total resistance, 
the correction for the gas phase resistance is not necessary 
in this experiment because the 0,/UV/Ferricyanide reaction 
is substantially slower. The same THES program is used to 
solve for K except the following differences: 
it The ozone consumption cannot be obtained from the gas 
phase measurement, the outlet ozone concentration is 
essentially equal to the inlet concentration owing 
to the slow reaction, the insignificant ozone usage 
makes the gas phase ozone mass balance impossible. 
The ozone consumption is calculated from the liquid 


phase cyanide mass balance by Equation (55). 


vhs) Kp has a value of 3740 Meececua: whieh mseanboutel/. 5 aL 
of the K, measured in the 0,/CN- reaction (19083 
MERA.) . 


The key equations for the flat CNO profile are Equation 


(105) 
K,tanhyDT-K, 
RT = ———————_ -: ARA > Cas LOS) 
DT-K, 
where 
Ky = K-CNF + Ky + 3.25 x 3740 x CNO 


K-CNF + Ky +7 L2.05>CNO 


= 


7 4 iw 


and Equations (55) and (52). K can be solved by using these 
3 equations in addition to the general basic equations from 
CS )eeto a( bO2)) invSectiony5... .The.algorithm is identical. to 
that used in the free cyanide case. 

RT is expressed in another form if a quadratic CNO 


profile is assumed (cf. Section 3.4.2) 


ARA+DTvb, be be \ 

RT = —————_— {52+ (Cy. - 52) Seer TR ATT LL (106) 
Sinhvb, 6 4 4 
4.25 K°-CNF + K 

= D 
b, =—-_—___ 
4 DT 

Rae 5 FF -CNO 

b_= L 


4 DT -ARA- 6 


where Equation (106) is obtained by combining Equations (85) 
Promivojvem es nquation (105) for the flat CNO profile is 
replaced by Equation (106) for the quadratic CNO profile, 
and the computational procedure would be the same in both 
cases. 

Equation (89) is used for the penetration model (cf. 


section 3.4.3) 


exp (-K, 9) 
RT =ARA-C /DT-K erf/K_DT [1+ ~~] + —————_ (89) 
s i uf 2K.98 —— 
1 Y7K, 6 
K = K-CNF + 12155 CNO + Ky 
4DT 
@ = 1Ke 
x 
eri. x = = | exp (-x7) dx 
v0 0 


ae 


94 


In addition to the subroutines MULR and FMU, another sub- 
routine INTG is used to integrate the function exp (-x*), 
which is declared as FUNCTION FCT(x) in the computer pro- 


gram THES to provide the error function. 


CHAPTER 6 
RESULT AND DISCUSSION 


aa | CcN_/0,/UV Reaction 

The input variables and directly measurable data are 
tabulated in Tables 6 and 7; data for a flow rate of 
0.765 ml/sec are in Table 6, and a flow rate of 1.094 ml/sec 
in Table 7. CNI, CNF are the feed and outlet cyanide con- 
centrations respectively. Py is the ozone partial pressure, 
T the temperature and RO the amount of ozone consumed per 
unit time, as calculated by the gas phase ozone mass 
balance. Utilizing these numbers, the rate constant of the 
free cyanide/O, reaction, Ke, can be evaluated by the com- 
puter program THES listed in Appendix 7. A brief explana- 
tion of the program can be seen in Chapter 5 Data Treatment. 
Four different methods as mentioned in Chapter 5 are used 


to obtain K and the results are shown in Tables 8 and 9, 


f! 
in which A designates the case based on the ozone consump- 
tion RT directly from the gas phase ozone mass balance RO; 
B, RT from the cyanide consumption, R, multiplied by an 

ozone consumption ratio of 1.5; c, RT from R multiplied by 


iaenand D, Ko. = K,/5.1 instead of RT. All cases are 


B 
regarded to be in the fast reaction regime. 
The variables and parameters of interest for the case 


D e.g. the physical mass transfer coefficient Kor cm/Sec; 


the interfacial ozone concentration CS, M; the calculated 
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cyanate concentration CNO, M; the reaction and diffusion 
time TR, sec, TD, sec; the reaction diffusion time ratio, 
TRD; the calculated ozone consumption 03, mole/sec; the 
measured and calculated ozone consumption ratios So and Sii 
the rate of cyanide reacted R, mole/sec; the fractions of 


ozone consumed by decomposition and cyanate oxidation FRD 


and FRC are printed out in TablesA4 and A5 in Appendix 6, 


Detailed Analyses of the Ozonation Reactions, where 


SeRO 
So = ae (107) 
Jy03 
Se ela (108) 
Ky 
PRD = 3a oe (109) 
Ky rae ag he Kp CNO + Ky CNF 


Ch K,°CNO 
- CNF 


Ky ney ey ae K,°CNO aS Ke 


FRC = : (110) 


Comparing the values of Ke from the four different 
methods (A to D, Tables 8 and 9), Terr SerOoundmeLoe GELicwcie 
from ozone material balance (A) does not furnish any useful 
information due to the scattering. Owing to the large 
uncertainty in the term RO, the amount of ozone reacted, 
some K, values are negative. The assumption of constant 
ozone consumption ratio is questionable. As the cyanide 
concentration increases, ozone decomposition becomes a less 
significant factor in the whole process and the ozone con- 
sumption ratio would decrease. The values of K, in B and C 


by assigning constant ozone consumption ratios and then 
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forcing the ozone mass balance are in the same order of 


magnitude as that of D, except that K, varies from some 


negative numbers at low cyanide concentrations to a value 


it 


GLaha, 000 M “Gee at a cyanide concentration of 6.6 x 1072 


M. It seems that models based on constant ozone consump- 
tion ratios are not valid. D is more acceptable, because 
a constant K-/K, ratio is assumed. 

Despite the heavy scattering of K, values, several 
conclusions can still be drawn. First, Ke is independent 
of cyanide concentration which is revealed by the plot of 
K, vs cyanide concentration (Figure 6) where a horizontal 


line can be drawn through the points. The average value of 


Ky is 19086 ut eecie with the standard deviation of +18%. 


The power law rate expression (first order in ozone and 
Cyanide or cyanate) and the validity of the model are con- 


firmed. Secondly, in contrast to the statement that "the 


(20) 


unnecessary UV may slow down the reaction" neither a 


beneficial nor a detrimental effect of UV on the reaction 
rate is observed which is in agreement with the observation 
that UV irradiation does not appreciably promote the ozone 


decomposition in the gas phase. 


K. is also solved using the same computer program by 


f£ 
setting Kp = 0 andeq-=s8ls5 tnstead or L523 and 3225 respec= 


tively. Ke remains unchanged for cases B and C in which 


a constant ozone consumption ratio is assumed, while 


it decreases almost 25% from 20,000 mt ees to 16,000 
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a 3% reduction in K, is observed by changing q from 3.25 to 
1.5. The values of Ky and q are chosen to be 15.3 sec + 

and 3.25 not only because they have more theoretical 

Support but also they give an average calculated ozone 
consumption ratio of 1.5 (S5 in Table A4) which is in exact 
agreement with the average measured value (S, in Table A4 and 


0 

So in Table A5 are not used because the ozone inlet concen- 
trations have not been measured accurately in these experi- 
ments). The small reaction time to diffusion time ratio, 
VER Dae Shasta nay in Tables A4 and A5 in Appendix 6 indicates 
that the reaction is in a fast reaction regime. The values 
of Ky can also be solved by the equations developed in the 
fast reaction regime (cf. 3.3.3 Fast Reaction Regime - Film 
Model). 

If the reaction is found to be in the fast reaction 
regime using the reaction time -diffusion time ratio as the 


criterion, there is a good chance for the reaction to be in 


the transition from fast to instantaneous reaction regime; 


further characterization is required. The dimensionless 
parameter 
D, 8B 
Vtp/tp/ (1 + D.C ) 


is used as the criterion to differentiate the fast 
reaction regime from the transition from fast to 
instantaneous reaction regime. The reaction 

is said to be in an instantaneous reaction regime for a 


parameter value of 10 or larger, a fast reaction regime for 
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a value of 1/2 or less, and a transition regime for the 
values in between. The value of this parameter at various 
cyanide concentrations lies between 3.2 to 0.75, therefore 
the reaction would be in the transition from fast to ins- 
tantaneous reaction regime. Two sets of initial guesses are 
used; the first set is good for Runs 1-3 in Table 6, and 
the second for Runs 3-21. The value of Ke calculated by 
assuming a transition regime is higher than that of the 
fast reaction regime by about 5% as shown in Table 10. The 
reaction, being in a fast reaction regime is still applica- 
ble. The liguid mass transfer coefficient at a cyanide 


3 


concentration of 6.66x10~-M is about 0.06 cm/sec, 


which is small compared with the gas phase mass transfer 


(85) The gas phase resistance, 


coefficient (2.15 cm/sec 
which amounts to 2% of the total resistance, can be dropped 
without affecting the result. The effect of gas phase 
resistance need not be considered for the complex cyanide 
as the reaction is in the transition from slow to fast 
reaction regime. 

As shown in Tables A4 and A5, the physical mass trans- 


3 cm/sec, 


fer coefficient, Kos is on the order of 2x10 
which is significantly smaller than the gas phase mass 

transfer coefficient (2.15 cm/sec). The calculated cyanate 
concentration CNO is only 1% of that of cyanide. The frac- 


tion of ozone consumed by decomposition, FRD, ranges from 14% 


to 34%, being more significant at low cyanide concentrations 


\tuvit .it¥s0 0) SE sewed yey einai 
ox heh: tent ont ai ad biucr ote 


: a | 


a it. 
isitivl to ete of » sh Be acho ed's esyosn Ri 


; ae 
t~l BHOR 307% Goop as Jae peers ont in ye 


asa? tadpirl et omtoedt moss fanwat a OM. 
wote a6 #@ 2apiis oe anipes noisoset 

Bi €@Mivey sotiore. sen? 6 ak puted ack 2692. 

ts tnatorttsso tehegsys segm Btopdl “t -okd 


" 7 ce 
war ers 
st M~ OL xad.d 26 nolveisnad Ter 


eel aap ots djinm becaegaigs coat eke 
(eB : ute 
‘owen o1,9) donio obtsi 


? 


teen iasoi ats Yo PS LS | ore: x 

ine2ie 2aAT ALizes ett tlornine 3 ros 
vr 

Mtasyo xefqaoo eft 102 Bessh.eaoo ed Jou been vonase tabn 


rt 


= | » ? « ter > i _ bias 
3 1 Yt nostianeisy oft wh. ws cotsose ae 


“eisai 2e6a leanieviq ands .2A bas bA eofdal at a 
i r % M 

.398\M> OL *45 to zebxo oft no Bs 4.8 

. oe | 


y ix 
220m seadq esp siz ned? ta Sa lane 


S7eaays hbessivolso efi? . (O98 \arp 2h.S).ain ‘ £e>. 
i sie air 
262% ade cohheaing 2 Yvan 20 ino -2k Om 
a pee: 


StL moxk espras va 0 aoa 


oe 
"2 6 

i: 
— 


7 Ses 


Run UV 
1 NO 
2 NO 
3 FULL 
4 NO 
5 NO 
6 FULL 
7 FULL 
8 NO 
9 NO 
10 NO 
11 FULL 
12 FULL 
13 NO 
14 NO 
1 NO 
16 FULL 
17 FULL 
18 NO 
19 NO 
20 FULL 
as FULL 


Ke Rate constant for free cyanide calculated by the transition 


regime, 


K. Rate constant for free cyanide calculated by the fast reaction 


regime (Case D, Table 8) 


Table 10 Comparison of K, for the cy"*/0,/UV reaction calculated 


by assuming the transition from fast to instantaneous reaction 


1.840 


3.701 


Gol 


7.694 


regime and fast reaction regime. 
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where the bulk concentration of ozone is high. Since the 
decomposition is the major side reaction and strongly pH 
dependent, an accurate measurement of pH is warranted. On 
the other hand, only 13% to 7% of ozone is destroyed by the 
cyanate reaction (FRC). In contrast to the decomposition, 
it is more significant at higher cyanide concentration. 

The total consumption of ozone by these side reactions 
ranges from 27% to 41% as the cyanide concentration 


eM tory eGo4 x 107° M. 


ancreases from 1.14x 10 
For an ideal case, RO, the number of moles of ozone 
reacted per second obtained by gas phase material balance, 


is equal to RT, the amount of ozone reacted per second 


calculated by cyanide balance, i.e. 
Rie= R(K_ [CN] ay oltre K, [CNO] + K,)/(K, ICN] ) 


However, the uncertainty in RO makes the comparison impos- 
Sible. This can be illustrated by plotting Sy» the ratio 


Cen Ounwagainst S~,the ratio of RT to Riin Figure 7 


1 
(So, S. are known as the measured and the calculated ozone 


1 
consumption ratios). The points show an uninterpretable 
pattern scattering in a random fashion, which complies 
with the observation in column A, Tables 8, 9. The 
scattering is natural considering the subtraction of two 
close numbers. The mass flow rate of ozone is on the order 


Orme one x ee mole/min, whereas the amount of ozone reacted 


is at most 1x 107° mole/min, which represents a 1/22 or 58% 
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Calculated Ozone Consumption Ratio 


1.0 3.0 


Measured Ozone 
Consumption Ratio 


Figure 7. Plot of the measured ozone consumption ratio 
So vs the calculated ozone consumption ratio 
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G12 


ey 


li 


a 
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decrease in the exit ozone flow rate. As the relative 
error is concerned, it is almost impossible to obtain RO 
by taking the difference between the inlet and outlet ozone 


flow rates according to the error propagation theory 


ROt= Nas ~ Na (212) 


where ae hag are the inlet and outlet ozone mass flow rate 
respectively. By a simple mathematical manipulation, the 


following relationship can be derived 


2 Nak | 2 An 7 
ARO 
(GS) = (+ + (4h) (112) 
fn tT hai — fag 


where ARO denotes the uncertainty in RO caused by the mea- 


Surement errors An... and Anns: Assuming a 5% of measurement 


Al 
error, 
2 = 5 ihe wD 
(FE) es gx (2202 E08 A 2.2) SM ays, (113) 
10 
ARO _ a ° 
Ro 1.56 = 1563 ‘ (114) 


156% only represents the inherent uncertainty of the term 
RO; the actual error including determinate errors and 

the errors from the data other than in and Nag may be even 
larger than that. In spite of the scattering, the average 
Measured ozone consumption ratio, Sor agrees with the cal- 
culated one, Si- They both indicate that an average of 


1.5 mole of ozone reacts with 1 mole of cyanide. 
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The error in K, comes from the following possible 


sources: 


i) 


=i) 


Experimental errors: Besides the measurement errors, 
additional errors:'can be caused by the irreproducible 
surface area and the liquid hold-up volume in the 
reactor. The vortex and ripple on the surface due to 
vigorous stirring may increase the actual surface area 
to a value larger than the geometrical cross section. 
In the ozonation experiment, liquid level was main- 
tained by manually adjusting the liquid outlet valve, 
an insufficient control in liquid level may change the 
geometrical surface area because the cross section of the 
OVF column is not uniform. Furthermore, the ozone 
partial pressure fluctuates with the atmospheric 
pressure outside the building probably due to the 
change in oxygen content and the air compressor inlet 
pressure. For the best result, it is advisable to use 
a liquid level controller to maintain the proper 
liguid level and compressed oxygen as the oxygen 
source to generate ozone. 

Errors from the chemistry of the system: The order 

of the reaction is still controversial, a value of 
Tan) | te or phar 46) (for cyanide) has been 
reported by various workers. The stoichiometric 
TecLo Jeno (moles of ozone required to react with l 


mole of cyanate) is uncertain too because a compli- 


- a 
cated mixture of Nj, NO3, NHg is produced. In the 
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present work, an order of 1, a stoichiometric ratio 

of 3.25, a rate constant ratio (cyanide/cyanate) of 
5.1 as observed by Balyanskii are chosen. The consis- 
tency of Ke values at various cyanide concentrations 
validates the choice, but still, the other possibili- 
ties cannot be excluded. 

iii) The assumption of flat concentration profiles for 
cyanide and cyanate may introduce certain errors in 
the calculation. Cyanide concentration may be re- 
garded as constant since it is the major species. As 
for cyanate, the same argument does not hold owing to 
the relatively small concentration. It is obvious 
that the concentration of cyanate is high inside the 
film where the reaction takes place and low in the 
bulk phase. 

In conclusion, the cyanide/ozone/UV reaction is in 


the fast reaction regime (t,/t, =5%* TONae with the 


first order rate constant of 19086 ut secs. 

Because of large uncertainty in the ozone gas phase 
balance, Ky has to be calculated from the cyanide 
liquid phase balance correcting for the ozone decom- 
position and cyanate oxidation. The film and pene- 
tration models serve equally well for the above reac- 
tion. The rate constant of cyanate, Kar which is 
used in the next section in the ferricyanide/ozone/UV 


reaction is equal to 1/5.1 of Ker ie. 93740 MeeSena 


-l1 -l 
for a measured Ky value of 19086 M Sch, 
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6.2 Ferricyanide/0,/UV Reaction 

The rate constants of ozone/ferricyanide/UV reaction 
at different UV intensities were calculated using three 
different models, namely, the film model with both flat 
and quadratic CNO concentration profiles and the penetra- 
tion model with a flat CNO profile. All cases have been 
derived mathematically in 3.4 Modelling of the Ferricyanide/ 
0,/UV Reaction. A single program can be applied to all 
cases with minor modifications in the form of Ki: The 
penetration model is more complicated and it requires one 
more subroutine INTG to generate the error function. Since 
the decomposition of ozone is important at high pH, the pH 
of the solution has to be measured accurately and fed into 
the Stumm correlation incorporated in the program to get 
Kp: On the other hand, at pH 7, a value of 8.5~x 169 sei 
as predicted from Stumm's correlation is accurate enough 
to account for the ozone decomposition. The main programs 
for film model, flat and quadratic CNO profiles, and pene- 
tration model all at pH 12 are listed in Appendix 7. Com- 
paring the programs for free cyanide and ferricyanide, 
striking similarity can be found. 

The input variables and the directly measurable data 
e.g. the cyanide contents of the feed and outlet ferri- 
cyanide solutions, CNI (6 times the ferricyanide concentra- 
tion), and CNF, temperature T, and the ozone partial 


pressure P are tabulated in Tables 11 and 12. Data for 


a 
pH 7 and a flow rate of 0.4944 ml/sec are reported in 
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LIRECTLY MEASURABLE EXPERIMENTAL VARIABLES 


Table 11 


FOR THE FERRICYANIDE/03/UV REACTION AT A FLOW RATE OF 


No 
No 
Coarse 
Coarse 
Fine 
Fine 
Full 
Full 


No 
No 
Coarse 
Coarse 
Fine 
Fine 
Full 
Full 


No 
No 
Coarse 
Coarse 
Fine 
Full 
Full 


No 
No 
Coarse 
Coarse 
Fine 
Fine 
Full 
Full 


No 
No 
Coarse 
Coarse 
Fine 


0.4944 ML/SEC, pH 7 AND TEMPERATURE 21°C 


Temp 


(°C) 


C5? Ca a se Le ea ee Gee es Ne a ee I a. 785. 6) See, Cs ae. S&S Ons S.- G2 Sg «& 3 on ed ae a or a 
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Ozone 
Pressure (PA) 


PREP PNR PRP RP BBR PRP PRP RPP RP RPP BP EPO RP RP RP BP PRP RP RP RP Re 


(atm x 102) 


-917 
- 933 
~925 
s97Z 
#92) 
e921 
-826 
.826 


944 
96% 
997 
-965 
964 
.006 
-977 
-905 


843 
843 
984 
«985 
- B99 
023 
+952 


-890 
.888 
.810 
+852 
-858 
~858 
=967 
.928 


084 
-960 
-961 
-937 
948 


Cyanide 


Feed (CNI) 
(M x 10°) 


RPRRrHre OOOOCOOCoCoOoOOoO eooocooo°ono eooqoooo°cneo oo°o°0o°0oec]e 


1843 
-1843 
1843 
1843 
.1843 
1843 
-1843 
1843 


73329 
#3329 
aa5e9 
3929 
~3329 
Migs FAY) 
73529 
~ 3329 


-4489 
4489 
-4489 
4489 
4489 
4489 
-4489 


-9609 
- 9609 
-9609 
- 9609 
-9609 
-9609 
-9609 
-9609 


394 
394 
394 
~394 
394 


Concentration 
Product (CNF) 


(M 


FPRrFrFrF ooaoaooocooaooaoo aoaoooooo aqagooooooo oo0o0o000 O° 


x 103) 


- 1583 
.1614 
- 1508 
#529 
1475 
1456 
«1372 
-1482 
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2854 
a2 423 
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- 2609 


- 3960 
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-4002 
-4061 
43937 
3741 
- 3666 


sOlh/ 
~8545 
8452 
8544 
- 8390 
-8469 
- 8358 
.8378 


240 
- 240 
Be Wf 
720u 
20 


115 


ay 


a * 


7 
RO ALERY SAGE <a 3 ah SEAS ae 
90 StWR WT ACTA a . ‘ 
S7 rE SAUTAAGE 


ant ames ehscerrsl) «yi Saran) 
(iO): Joebert (Hj) hanw 


+ 4 7s FSS! mf) 
Ef) eee 
E64 «G ehef.0 
TS & ak GD 
i.e £ xa 0 


“>. a 


rte 2 

Eien f) 
. = ¢B.Q 
iiea.0 
O8ES.0 
@388,0 
Ret B.0 
ot 


116 


Table 11 (cont'd.) 


UV Temp Ozone Cyanide Concentration 
Pressure (PA) Feed (CNI) Product (CNF) 
- [aG) (atm x 107) (M x 10%) (M x 10%) 
Fine 23.6 1.934 1.394 1.245 
Full 25.0 1.942 1.394 L2G 
Full 25.0 1.907 12394 P2354 
No 23.4 1.936 Zee e 2.044 
No 2 ie 1.908 Deere 2.042 
Coarse Fe 1.945 Ziatew we Pl 09 be 
Fine PAS Neg! eta 2n2i2 2.080 
Fine 2a.4 1.842 A aeep ay 3 1.996 
Full 23.4 1.881 PE i Wells 
Full 204 1.881 Zacks 2.084 
Full 23.4 eo ae PENILE 2.093 
No 24.3 Teo 06 Sey: 52 E19 
No tA | 1.924 5.487 epg Ba RA 
Full tts Pa 1.929 5.487 aloo 
Full 24.2 De Gy De SOy Sooo 
Full 23.4 1.860 9.342 9.007 
No 2350 1.968 pA OA oo 
Full 23.8 1.943 LOSZE 9.805 
Full 23, 0 1.994 LOZ TOL 
No 26.1 1856 17290 Leo 
No Zo. E 1.890 17590 17346 
Full 26.1 1.874 17.90 Leet 
Full 25.0 1.927 Ghee iy fet fi) 
No 24.0 1931 36.74 B07 OL 
Full 24.0 1.867 36.74 36.10 
Full 24.0 1.905 36.74 35240 
Full 25 oe 1.908 30502 35.92 
Full ae, 2.000 36.82 SO OL 
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Table ll. Table 12 gives data for pH 12 and a flow rate of 
0.5525 ml/sec. The values of K calculated from the three 
different methods are shown in Tables 13, 14. kK, KQ, KP 
denote the rate constants based on the film model with 
flat CNO profile, film model with quadratic CNO profile and 
penetration model with flat CNO profile respectively. 
Comparing the values of K from the three different 
models, it is found that K is larger than KP by 8% or less, 
whereas KQ is slightly higher than K at low cyanide content 
and substantially lower than K at high cyanide content. 
The deviation at the high end may be as large as -100%. In 
order to show the trend more clearly, the rate constants 
from three different models, i.e. K, KP and KQ are plotted 
against cyanide content at various UV intensities on 
Figures 8, 9, 10. K can be correlated with cyanide content 
and UV intensity. In Figures 8, 9, K and KP at 100% UV are 
almost independent of cyanide content as represented by 2 
horizontal lines (1100 for K, 1000 for KP) in the graphs. 
The effect of UV on the rate constant is not obvious until 
the cyanide content is below 1.5*x ihn M, i.e. 40 ppm, then 
the rate constant decreases sharply with depletion of 
cyanide in the absence of UV. At a fixed cyanide content 
(below 40 ppm), the rate constant increases gradually and 
finally reaches a plateau value as the UV intensity is 
increased to 100% (200 watt). While KP agrees with K within 
8%, KQ shows a substantial discrepancy (Figure 10). The 


whole curve is skewed relative to that of KP and K, and 
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Run K, 10M 4s xg, 10% stoke, 102M As 2 
1 0.9238 1.271 0.8415 
2 0.3072 0.3898 0.2888 
3 D5 Git 5 Aah 2s 370. 
4 1.820 1 De 420 1.636 
5 3.805 4,802 3.401 
6 4.534 5.621 4.047 
7 10.29 11.68 9.191 
8 4.643 5.807 4.156 
9 2.193 2.735 1.960 
10 2,033 2.552 1.817 
11 4.967 5.589 4.427 
12 4.159 4.794 3.706 
13 7.4168 8.242 6.945 
14 7.406 7.906 6.621 
15 12.31 12,30 11.10 
16 12,53 12.50 11.31 
17 4.202 1633 3.752 
18 0.7821 10.41 0.7040 
19 3.130 3.580 2.790 
20 2,187 2.617 1.949 
21 4.683 5.088 4.185 
22 9.523 9.378 8.609 
23 10.48 10.19 9.502 
24 8.048 7.189 7.430 
25 8.557 7.580 7.917 
26 10.98 9.388 10,24 
27 9.029 7.932 8.367 
28 11.68 9,910 10.92 
29 10.25 8.855 9.547 
30 11.02 9.449 10.29 
31 11.01 9.434 10.28 


Table 13 Comparison of rate constants K,KQ,KP from three 
different models for the ferricyanide concentrations 
listed in Table ll. 
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Table 13 continued 


32 9.827 Seiit7 9.282 
33 TTT 8.768 10.20 
34 20.68 al St 19.97 
35 9.213 7.659 8.679 
36 8.543 wealeeA 8.027 
37 9.851 ie hy 9.302 
38 10.52 8.588 9.961 
39 ieee 9.662 11.47 
40 12.59 9.361 12.18 
41 13.08 9.666 12.66 
42 9,532 75352 9.150 
43 9.588 Teor be p! 
44 19.45 13.59 18.98 
45 6.363 5.169 6.027 
46 9,193 72108 8.811 
47 8.665 6.744 8.289 
48 bie 0 N25 TEES 
49 10.30 6.899 10247 
50 9,331 6.336 9.143 
51 9.378 6.366 9.191 
52 6.411 4.216 6.286 
53 7.854 5.010 7.748 
54 F227] 4.684 There.» 
55 6.621 4.330 6.518 
56 9.120 5.174 12720 
57 5.349 3.283 5.286 
58 7.093 aay 9.062 
59 Ts s02 6.417 12.99 
60 13.02 6.417 12.99 
61 4.814 opie 5.703 
62 10.80 5.450 10.76 
63 9.226 4.729 9.189 
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Run gy NA ROME Os Manes ade sinalt eee 

a 9.911 9.366 9.842 

Z 14.57 LZ O1 jWey sv. 

3 OL 9.910 DZ eae 

4 10.68 9.037 Pa ee 

5 25.75 18.35 25200 

6 9.387 7.685 11.62 

a) 14.57 11.18 17.81 

Table 14. Comparison of rate constants K, KQ, KP from three different 


models for the ferricyanide concentrations listed in Table 12. 


¥ 1 


qi - a a 
fi 
7 Di 
a 


we > sae 2 
pa to NOE 
- treqme 
woizye ner 
it off 
1% 
x02 


Lee 


3000 


S) 
° 
3S 


© 100% UV (0.319 watts) 
© 40.4% UV 


O 15.5% UV 
A 0% UV 


T= 22°C 


Rate Constant, K, M7! sec” 


ES ES ESES ES ESES ES> ES 
ay Sige dem por reky fo Go Ba Qu By 
oO fo) ro) fo) Oo fo) fo) o) fo) 
or Dips Oa One On alte Sep Se 
N Nx OW x SW Boo 3g leh $0 ak rete te 

Oo ~ © © xe) Oo - © 
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Figure 8. Plot of K from the film model, flat CNO profile vs 
cyanide content, CN, at 90° Os pH 7, 108 rpm and 
various UV intensities. 
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Figure 9. Plot of KP from the penetration model, flat CNO 
profile vs cyanide content, CN, at 22 C, pH 7, 
108 rpm and various UV intensities. 
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Figure 10. Plot of KQ from the film model, quadratic CNO 
profile vs cyanide content, CN, at 22 C, pH 7, 
108 rpm and various UV intensities. 
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KQ depends strongly on cyanide content. It implies that 
the model has been over-corrected by the assumption that 
the CNO profile is quadratic, and the rate constants 
obtained by the assumption of an uncorrected flat CNO pro- 
file are more reasonable. kK and KP at concentrations 
higher than 40 ppm (1.5x 107° M) and below 40 ppm with 100% 
UV irradiation are independent of cyanide concentration. 
The effect of UV on the rate constant can be seen more 
clearly by transforming Figure 8 into a ln kK vs $% UV plot 
at fixed cyanide concentrations (Figure 11). K increases 
with UV intensity as well as cyanide content. The curves 
converge to the horizontal line at high cyanide level 
(1.5x 10° M, or 40 ppm). 

Samilarly, Keand KQ at pH"12, 100% UV are plotted 
agaenstecyanide content at 22°C. In Figure 12, a horizon- 
taleinne ws observed/tor’ K (K=1100) and a severely bent 
curve for KQ, as represented by a dotted line in the graph. 
The graph reveals that the flat CNO profile is more satis- 
factory than the quadratic one. The value of K at pH 12 
agrees with that of pH 7. The agreement implies that the 
model and the assumptions made in the model are valid, 
therefore the value of K is independent of pH and ferri- 
cyanide concentration. 

The fact that the rate constant apparently decreases 
with decreasing cyanide content at low cyanide concentra- 
tions in the absence of UV imposes a problem for the treat- 


ment of ferricyanide. The rate of oxidation diminishes in 
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Figure 11. Plot of K vs UV intensity at 22°C Spa, 1 08erpm 
and various cyanide contents. 
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Figure 12. Plot of K, KQ vs cyanide content, CN, at 0 and 
100% UV, 22°C, pH 12 and 108 rpm. 
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the final stage where cyanide is almost depleted; the oxi- 
dation is so slow that in the absence of UV the desired 
discharge standard of 1.5 ppm cannot be achieved with a 
reasonable contact time and flow rate. To facilitate the 
oxidation reaction, UV has to be applied to the reaction 
System at low cyanide content (40 ppm) where the effect of 
UV is pronounced. For instance, at 1x 1077 M (2.6 ppm), the 
rate constant with 100% UV (K=1100) is higher than that 
With 40°4¢°UV\ by. a factor of 3, 15.5% UV -by 9)#and 0%) UV 
veel s 0,/UV is indeed a powerful combination to reduce 
the wasted cyanide content to the allowable discharged 
level. 

Most of the important parameters are shown in Tables 
A6, A7 in Appendix 6 for pH 7 and 12 respectively. At all 
UV intensities, K, is approximately 2.2 x 107? cm/sec and K 
ranges from 1.7 x LOM. 


L 


£0 36077 cm/sec, which means an 
enhancement factor of 7.6 to 13 irrespective of the pH of 
the solution. The values of EN are considerably lower than 
those for free cyanide, for which EN may exceed 20. The 
value of TRD (~0.5) indicates that the reaction is in the 
transition regime. The calculated cyanate concentration, 


6m to 9x ji Meka a Ncontent 


CNO, varies from 3.6 * TOn 
increases from 1.8 x Tlie LOuo so 10° °M, which is also less 
than that for free cyanide at the corresponding concentra- 
tion because of slower reaction. The trend can be observed 


more clearly on a log-log plot of CNO vs CN as given in 


Figure 13. A linear relationship exists and the lines for 
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free cyanide lie above those for ferricyanide. 

The percentage of ozone decomposed, FRD, shown in 
Tables A6é, A7 depends strongly on the PH of the solution. 
At pH 7, ozone decomposition is not significant and accounts 
POre0.2)L0)0.001% of the total ozone consumption as cyanide 


4M to 3.6x10°° M. 


concentration increases from 1.3 x10 
The decomposition of ozone can thus be ignored, so that 
accurate measurement of pH is not required; a Kp value of 
Osa Xx ia can be obtained from Stumm's correlation and was 
used throughout the program. However, decomposition 
becomes the major reaction as the pH is raised to 12, under 
this condition, decomposition claims 70~90% of the ozone 
consumed. Utilizing ozone for waste water treatment at 
such a high pH is uneconomical. FRD for ferricyanide and 
free cyanide at pH 7 and 12 are plotted against cyanide 
content in Figures 14 and 15. Comparing the FRD for free 
cyanide with that for ferricyanide the decomposition reac- 
tion is still significant (40%~10%). However, it is far 
less than that of ferricyanide because the faster cN /0, 
reaction tends to reduce the ozone concentration in the 
solution. For all cases, FRD decreases with increasing 
cyanide concentration which is in agreement with the model 
used. An increased cyanide level will reduce the ozone 
concentration in solution, and thus diminish the decom- 
position. 

| The fraction of ozone reacted with CNO, FRC, is 


complementary to FRD, large FRD is accompanied by small 
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Figure 14. Plot of fraction of ozone decomposed FRD vs CN for 
ferricyanide at pH 7 and 12, 100% UV and for free 
cyanide at pH 12, flow rates of 0.4944 and 0.765 
ml/sec. 
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Figure 15. Plot of fraction of ozone decomposed FRD vs CN 
for ferricyanide at pH 7 and 12, 100% UV and 
for free cyanide at pH 12, flow. rates of 0.4944 
and 0.765 ml/sec. 
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FRC and vice versa; and it accounts for the graph obtained 
by plotting FRC vs cyanide content (Figure 16). FRC is 
determined by the cyanate concentration as well as the 
ozone concentration. Even though CNO increases with CN, 
FRC declines because the decrease in ozone concentration 
can more than compensate for the increase in CNO. How- 
ever, the case of ferricyanide at pH 12 holds an opposite 
trend, it is probably due to the fact that at highnnoH, 
ozone has already decreased to such a small extent that 
itis not sensitive to the cyanide content, and then the 
Shape of the curve is completely governed by CNO. 
Conclusion 

Comparing the ferricyanide/O,/UV reaction at pH 7 and 
12, the dissimilarity is quite striking. At pH 12, ozone 
decomposition is dominant and consumes 90~70% of ozone, 
0,/CNO- and 0,/ferricyanide reactions only claim the remain- 
ing 3~10% and 7~20%. On the other hand, at pH 7, 03/CNO~ 
is as important as the O,/ferricyanide reaction (70~30% 
and 30~70%), and the decomposition is negligible. Despite 
the adverse predominant reactions and diversified condi- 
tions, after correcting for decomposition and CNO oxida- 
tion, the same K can be obtained at both pH values (K=1100 
at 22°C, PH 7 and 12). The excellent agreement as well as 
the applicability of the same model to the diversified con- 
ditions prove that the model proposed is indeed valid, and 


the argument in the previous chapters is substantiated. 
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6.3 Proposed Mechanism for the Ferricyanide/0,/UV Reaction 


The influence of UV and cyanide content upon the oxi- 
dation rate constant can be explained by a plausible 
mechanism involving the formation of an ozonide inter- 
mediate. It has been known in organic chemistry that ozone 
Can react with alkene to give an ozonide. The reaction of 
ferricyanide with ozone may also involve the following 


ozonide as an intermediate. 


6- 
C SF ll N—— ¢ 4 Fe(cN), 
i: x 
ys af 
Sa Sol 
ii rent 


Only a single form of ozonide would exist in that rearrange- 
ment by breaking the C=N bond is not feasible due to the 
double bond character. Furthermore, the C-O bond in the 
ozonide is strengthened by the resonance form II shown 
above, because of the 6+ charge developed at the carbon 
atom as a result of the stable bond between C and Fe and 
the electron withdrawing properties of Fe. The stability 
of this ozonide would account for the slow oxidation rate 
of ferricyanide. However, the decomposition of the ozonide 
could be catalyzed by another CN group, as it would provide 


a positive center to facilitate the electrophilic reaction 
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(NC) _Fe-C =N (NC) _Fe-C =N 
5 5 | 
SNM cabins 
O O 
\ VS ay if 
OG) 2a aife = N (NC) ,Fe + CNO 
OH 
Fe(CN) > + OO) + CNO a 
6 2 
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vs Fe (CN, (OH) ) — 


Schematically, the reactions can be written as 


k 


Ane Bs 


Ky 


k Ko 
B +B— > Ct+B 


where A denotes ozone, B cyanide and C cyanate. 


Rates can be expressed as 


Tr. a kiO,C, + k_1Cp 
a ee owcls Ebi de uel 
Bee Pon Cp ee 4 Gr 28 BAB C* 
Using steady state hypothesis, the term ths 
to be zero, then 
c* = k1C,°p 
Ea icmre tii Co 
2 
. Geb hyi we 
“<< EM jhavs Koeoct «kxG 


(5) 


(116) 


can be assumed 


CLL?) 


(118) 
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Several cases can be discussed: 


i) As cyanide concentration is high, k Cp >? k 


2 mo eo 


eee = kiC,C, (119) 


which is in agreement with the constant values of K determined 


pa GIN eee eos nto M (40 ppm) by the horizontal line in 


Figures 8, 9. 


ii) As cyanide concentration is low, k5C, << kia and 
k,k 
1 Z 
r = — — <= 
B mae CAB (120) 


iii) If UV stabilizes the ozonide intermediate and thus 
reduces ki or promotes the bimolecular ozonide reaction 
(ko), then k5C, >> ki 


The rate expression would have the same form as that in i), 


at 100% UV even at low cyanide level. 


and the rate is first order to both ozone and cyanide con- 
centrations. If the UV intensity is less than 100%, Ky is 


suppressed but still larger than koCar then 


Lae? 22 os 
rp = say Cr.Cp = KCC, (302) 
Kok 
ae K ue i 2 ( a (122) 
Sas B 


From this expression, K is expected to depend on Cy nel 
a linear fashion, and to increase with UV intensity. Indeed, 
it is what has been observed in Figures 8 and 9. The slopes 
of the curves on a log-log scale at cyanide concentrations 
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03 UV tees 
Loe. UV Bi al 
40.43 UV 0.95 


which are close to the value of 1 indicated by Equation (122). 

The reaction mechanism proposed gives a satisfactory 
explanation of the K values obtained, however, it should be 
noted that this is not the only mechanism which can account 
for the experimental K values. At this time, it is somewhat 
naive to claim that this mechanism is realistic, unless the 
actual experimental work iS carried out to positively iden- 
tify the intermediates. Further work has to be done to 


clarify the mechanism. 


6.4 Other Possible Mechanisms 
6.4.1 Free Radical Mechanism 


It has been reported that ozone dissociates upon the 


absorption of ce oe 
os 
Oo. + I —~ O+0 
3) a Z 
DS) 


The atomic oxygen produced in the dissociation process is 
an active species and can react with ferricyanide as repre- 


sented by the following reaction 


ee 
5 


-3 Ie = 
8S CN. ——> CNO + CN 


where cne> and cn,” denote ferricyanide and its product 


after one cyano group has been oxidized, and I, the UV 


intensity. Using steady state hypothesis 


‘aj 
‘At 


> 2THAOLEMTE 4 
: a) Shad 

3Vip beaoqna mere ee moisgnex It ‘ie 
teyvewr:d Bbanreagda saulay 4 ef3 io sata ri 


sada 


<= 


| 

mi 
- i pA ies 
1B! itw webomicom yina ait Sow at ata bat 


h» 
Ps —~ 


it ,gmMts eins PA .euehaev 5 Lasasmtiogxe 3 50 


Or iveiisnt 2) minsdoew ekdd peas ataie: oe evi 
i el te aoe 
oF cuo beatuse of Dhow lsdneatagxs teu 
ty wits «86 eegtet bowssad orl: 
al 
i rerio +a 
ae 7 
ia g ae sagt APD E BSS 
t j ) Ita Jira 


22950%q nottsiouzeip sts we 


4 25 ShIK6YS i358) vi2gew PoReey 


ope, eth ae ape 
Wish (are 


L329 


-~ 


Fy = ~k,I[0,]~k,[0] [0] -k,[0] [on 


x. 


] = 0 (123) 


Zonz3 = cr CONLCN ee (124) 


By solving these two equations, the rate of oxidation of 


ferricyanide can be expressed as 


5 
kk, [03] [CN, JI 


ro .30 = = -K'[0,] [CN] (125) 
CN =3 3 6 
6 k,[0,] +k, [CN, ] 
ely ede 
K' = 13 —- (126) 
k,(0,) + k,[cN,7] 
kyI 
lim K! = TT ETE (127) 
[CNZ*] += el 
ki kaI 
ee eels (128) 
[cn,71+0 k, [03] 


where K' is the lumped pseudo rate constant for the reaction 


assumed to be first order to 03 and Sita The effective in- 
tensity of UV in the reactor depends on the ferricyanide 

concentration because ferricyanide solution absorbs UV. The 
dependence of K' on [CN, 7] and I can be shown schematically 


by the following diagram. 
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6.4.2 Mechanism Involving Ferricyanide Dissociation 


JED, 


Ohno indicates that ferricyanide will dissociate 


into a free cyanide ion and an intermediate containing 5 


cyano groups. 


L 1 . z 
H.O+Fe(CN)> +I => Fe(CN).(H.0) 2 + CN 
2 6 : Bo 


2 


The CN produced is responsible for the oxidation reaction 


-_ k, - 
Hee > CNO P+) Onaae 


Similar to the approach in 6.4.1 Free Radical Mechanism, 


the rate expressions can be written as 


Boy) NS TE LCN, *)-k, [CN ] (CN, *]-k, [CN 1 [0, }=0 
(129) 
r 13 ='- k ,LICN, *] +k, (CN, 21 CNT] (130) 
CN 
6 
k,1 
: tae =k, [CNZ >} (0,){ 3, ____| (131) 
o° CNG ko {Co [CN 7] }4+k, [0, ] 
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pete 
k,{Co- [CNG7]} + k3105] 


where Co is the initial ferricyanide concentration, and is 
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k4I 
cae Rue: i = (134) 
[CN, i 2 ICN, ] + k,[0,] 
K' T 
=-3 
[owe] 


The shape of the K' vs [CNT 7] curve for the mechanism 
involving ferricyanide dissociation is similar to that of the 
free radical mechanism. The only difference is that K'in the 
former case decreases with increasing ozone concentration; 
in other words, the reaction order of ozone in the ozonation 


of ferricyanide is less than l. 


6.5 Conclusions 

Comparing the experimental K vs [Fe (CN) 77] plot with 
that from different possible mechanisms, the molecular re- 
action one through the formation of ozonide seems to fit the 
data best. However, the other two mechanisms cannot be ex- 
cluded. The dip in the experimental K vs [Fe (CN) 27] plot 
(Figures 8, 9) remains unexplained by the ozonide mechanism. 
Since both the free radical and ferricyanide mechanisms 
predict a downward trend for K at high ferricyanide concen- 


trations, this dip may be due to the presence of these two 
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mechanisms together with the ozonide one. The free radical 
and ferricyanide dissociation mechanisms will not be dominant 
at high cyanide concentrations, otherwise K would depend on 


the UV intensity and decreases with increasing ferricyanide 


concentration. 
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CHAPTER 7 
FUTURE WORK 


UV/ozone reaction systems have received wide atten- 
tion this year (1977) as revealed by the number of papers 
presented in the International Ozone Institute Symposium 
held in Toronto ‘86? | The present trend shows clearly that 
ozone is a potential chemical which can replace the conven- 
tional chlorine in the foreseeable future in the municipal 
and industrial water treatment. The main barrier to bar 
the further usage of ozone in water treatment is the cost 
which can be overcome gradually with the advent of techno- 
logy. However, in order to utilize ozone more effec- 
tively, the ozone reaction mechanism should be studied 
more rigorously. The ozone generating technology and ozone 
reaction mechanism are the two important areas that should 
be worked on. 

As for the continuation of this piece of work, the 
following studies are recommended: 

a) Most of the difficulties in the experiment arise from 
the uncertainty in the gas phase ozone balance due to 
the limitation of the surface area in the reactor. 
The transfer of ozone can be enhanced and thus the 
difference in ozone concentrations in the inlet and 
outlet streams be detected easily by a reactor of 
larger interfacial area. A bubble reactor capable 


of providing reproducible bubble size would be a 
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good choice. The interfacial area can be obtained by 
correlation or by carrying out a well-documented 
simple reaction, and the same surface area can be 
used for the ferricyanide/0,/UV reaction should the 
experiment be carried out under the same conditions 
used in the surface area calibration. RT, the ozone 
consumption can be measured directly with good accu- 
racy which allows the solution of K and K, simul- 
taneously. 

In the treatment of data, ozone decomposition is 
assumed to be in parallel with ferricyanide oxidation, 
therefore ozone decomposition will decrease the effec- 
tive ozone concentration in the solution and thus 
decrease the ferricyanide oxidation rate. According 
to the free radical mechanism, ozone decomposition 
would be in series with the ferricyanide oxidation 
reaction and thus increase the oxidation rate. The 
free radical mechanism can be tested by using the 


same computer program THES except that K is given by 


' -3 
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ky [0.] + k, (CN. ] 


where Ip is the incident UV intensity, e' is the molar 
absorptivity multiplied by the distance from the 
bottom quartz plate of the reactor to the interface. 


[0] can be assumed to be equal to the saturated 


concentration at the partial pressure of oxygen in 
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A sophisticated program can be written to include the 


free radical mechanism by assuming the following reac- 
tion scheme: 


0, + CN, ———-+ CNO molecular mechanism 
CNe i 
0, meio O ——> CNO free radical mechanism 
et 0. ozone decomposition 


The contribution of the ferricyanide dissociation 
mechanism to the ferricyanide/O,/UV reaction can also 
be tested by running the reaction at different ozone 
partial pressures. If the rate constants calculated 
by the program THES decrease with increasSing ozone 
partial pressure, the ferricyanide dissociation 
mechanism would be present (cf. Equation (132)). 

The rate constant and the stoichiometry of CNO /O, 
reaction should be elucidated. This can be done using 
the same set-up as mentioned previously except that 
KCNO solution is passed through the reactor instead 

of K,Fe (CN) _. 
The molecular mechanism proposed in Section 6.3.1 has 
to be tested by identifying the ozonide. 

The rate constant obtained in this work can be cross- 


checked with the stop-flow technique developed in 


chemistry. 
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£) The same technique can be applied to other cyano 
complexes such as nickel, cobalt, silver, cadmium 
Cyanide complexes. The rate constant might be able 
to be correlated with the instability constant of the 
complex. Then the rate constant of a cyanide complex 


can be predicted by the known instability constant. 


Recently, the arsenic and cyanide problem in the gold 
mines operated at Yellowknife became a hot issue and 
received criticism from the public. Ozone finds its imme- 
diate application to this problem. Arsenic can be oxidized 
by ozone to the Ag(V) state, which will form a calcium 
arsenate precipitate at high pH. The treatability as well 


as the rate constant should be explored. 
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NOMENCLATURE 


& Concentration 

CAO Bulk concentration of A 

ey Effective ozone concentration in solution 
CNF Bulk concentration of cyanide 

CNI Feed concentration of cyanide 

CNO Concentration of cyanate in solution 

CS Interfacial concentration of ozone 

D Diffusivity 

DT Diffusivity of ozone in solution 

E UV intensity, W 

EN, En Enhancement factor 

F Flow rate 

FRC Fraction of ozone reacted with cyanate 
FRCN Fraction of cyanide destroyed 

FRD Fraction of ozone decomposed 

H Henry's law constant 

I, UV intensity, quanta per second 

J Molar flux crossing a unit surface without 


chemical reaction 

J' Molar flux crossing a unit surface with 
chemical reaction 

K Rate constant for the O,/ferricyanide 
reaction using the film model with a flat 


CNO profile 
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Rate constant for the O3/ferricyanide reaction 
using the penetration model with a flat CNO 
profile 

Rate peneran for the 0,/ferricyanide reaction 
using the film model with a quadratic CNO 
profile 

Rate constant for the reaction between A and B 
Rate constant for the 0/CNO— reaction 

Rate constant for the ozone decomposition reaction 
Rate constant for the 0,/CN- reaction 

Gas phase mass transfer coefficient 

Liquid mass transfer coefficient 

Piysicalemass transier  coehiicient 

Overall liguid mass transfer coefficient 
Optical path length 

eee formed upon UV irradiation 

Ozone flow rate 

Partial pressure 

Stoichiometric ratio 

Amount of A or B reacted per unit time and unit 
volume 

Amount of cyanide reacted per second 

Amount of ozone consumed obtained by direct gas 
phase ozone mass balance 


Amount of ozone consumed by calculation 


Surface area 
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Measured ozone consumption ratio 
Calculated ozone consumption ratio 
Moles of ozone required to react with one mole 
of cyanate 

Temperature 

Time delay 

Reaction time 

Residence time 

Diffusion time 

Ratio of reaction time to diffusion time 
Reactor volume 


Distance from the surface 


Greek Letters 


6 


8 


Subscripts 
A 


B 


Thickness of the film 
Contact time 
Molar absorptivity 


UV frequency 


Gas component, i.e. ozone 

Liguid reactant, i.e. cyanide or ferricyanide 
Product of the reaction between A and B, i.e. 
cyanate 

Exit stream 


Gas phase 
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Gas phase 
Feed 

Liquid phase 
Bulk phase 


Surface 
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Appendix 1. Mueller's Method 


This method 5455) has been listed in IBM SSP scienti- 
fic package bearing the names of RTMI and DRTMI (double pre- 
cision). The whole subroutine is used in the computer pro- 
gram as MULR for free cyanide and ferricyanide with minor 
modifications and a change in the dummy arguments. RTMI or 
MULR are essentially the subroutines used to determine the 
root of the nonlinear equation f(x) =0 in the range of x 
from s-osup to *x 


Li 


tion scheme of successive bisection and inverse parabolic 


Ri (by input) by means of Mueller's itera- 


interpolation. The procedure assumes f(x, 5) -f£(x,,) = Atle, 


Starting with X= XG and Xp = Xpis the middle of the inter- 


xX. is computed as shown in the following diagram. 


val Xr0 R 


/ inves roti 
In case f£(x_)°f(x.) <0, th pate oe 
: m R / 


/ 


° #(x (2. step) 
x, and x, are interchanged yeh ke 


to ensure that 


(x) (4. step) p—— de 


£(x,,)°£(x,) a0. 2 


mo x, (1. step) 
an x, (2. step) 
(2. step) 


x) 


Mueller’s iterative method 
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in case_of rapid convergence, i.e. the root is neither too 


close to xX, nor x, as described by the equation 
Zo ie - - a 
(x1 (x) E(x, )] £(x,) [f(x,) E(x, )] md, AN OA), 


Xp is replaced by Xn and the bisection step is repeated. 

If, after a specified number of successive bisections, the 
above inequality still persists, the error message is sent 
out. In case of slow convergence, i.e. the root is near to 
either X, OF Xp, SO that the above inequality is not satis- 


fied, then the new point is generated by an inverse parabolic 


interpolation 
xou x i -2f +£ 
Ox = £(%)) EEToE GC) pre () TET ee Teeter — T 
m it aR m on ra 
(A3) 
and x =x, — Ax. Xx is sure to be situated between Xp and xn 


In the next iteration step, x becomes Xr and Xn becomes 
Xp aot E(x)+f£(x,) > 0 or x becomes Xp pipe £(x)-£(x,) < OfeetThe 
convergence is either linear or quadratic if the multipli- 
city of the root to be determined is equal to 1 or greater 
than 1 respectively, and if f(x) can be differentiated 
continuously at least twice in the range X15 and Xpis Each 


iteration step requires two evaluations of f(x). This iter- 


ative procedure is terminated if either the two conditions 


below are satisfied: 
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i) |x,-x, | < €-max(1, |x, |) check in 
bisection (A4) 
and Jf (x,) - £(x,) | Sol Nel ets loop 
ii) Pome re max (), > |x 1} check after 


inverse parabolic (A5) 


and |f(x)| < 100-ce interpolation 


Modified Mueller method 
The modified Mueller method is an extension of 
Mueller method which is used to solve two nonlinear equa- 


tions with 2 unknowns. Considering the following equations 
0d tt I BG (A6) 


y = f£,() (A7) 


with the solution of (XorYs)e then a function f(x) is 


defined as 
ax) £, (x) - £, (x) (A8) 
ro eee Tet 6 
Ss Ss 


so Mueller method can be used to obtain Xo by solving the 


single equation f(x). 
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Appendix 2. Detailed Operating Procedure 


General Operating Procedure 


The procedure outlined here is general for all experi- 
ments pertained to the text with only minor variations for 


individual experiment. 


1) Wet Test Meter Calibration: The wet test meter was 
calibrated by a device of known volume. Sufficient water 


was added into or drawn out from the meter until a cubic 
foot of displacement volume corresponded exactly to one 
revolution on the meter. A black adhesive tape was placed 
on the meter chamber to indicate the correct water level 

as it might fall gradually as a result of evaporation. 

2) Leakage Checking: Distilled water was pumped 

through the reactor with stirrer on, bearing cooling water 
On and air regulator valve opened. The reactor was tested 
against leakage by pressurizing to 1.5 atm with air, and 
then the pump, and valves T3, T2 and T4 were shut down,with 
Tl turned to the position 'VENT'. The pressure should hold 
for an infinite period of time if the reactor was perfectly 
sealed. If so, the needle valve would be adjusted as the 
next step, otherwise the bolts holding the reactor had to be 
tightened until the leaking was stopped. 

3) Needle Valve Adjustment: 350 ml and 125 ml of dis- 
tilled water were added to the first and second bubblers 
(Bl and B2) respectively; the water level was also marked 


by the black tape. Both Bl and B2 should be appropriately 


aos] 


-2Ib Yo Lm 


: | ‘ti my + 
ibbesest 9c 7 apa Baa - Pr fn 


semanas 


‘> Iie 12 (exeanep at »isri bestia suneoom ae ae) 
. 


moitsliev szonkm yine dsiw Seg, | aad ow bent st 3s 


<san 09pe | f e 
al | 


even Jeo OW o;fT 2MOLI6 awit i fed ‘q036M. 
mate) 
foittva .saulov awondt Fo ative’ 6 va | bac 


all rae 
° 
ina’ 1Ssd6e0 sd? mows Seo fwexh a ovat bebbe 
i ber a 


: jiguexe Debrovesrt00 smulov jemarne nS 


: a aN i 
fleonbs zposid A&A .tezem odd 20-80 PULO 


‘9 to tiveaa s es yilesbsxp Stes 

cit: aew wersw belilvara paiaoee> 2pBAs 
prrlooco PRiTsed a | xi232 it kw to20867 i 
1099893 eff .beaago syviav witente 4 

7 


¢ 
Mere 


tis dziw mse ¢.1 oF pnisixasesowq yd oon 


wie ssaw $f bas ST ET soviny 5 
’ *s 


vode Surgeitq eff ,'Tway* Rolsivog & 
s28q 200 TeyoewT oY Bh emis to Beg 


1 betesthe e€ biuew sviav oben > ss 


- + Jigen awe 
ei ies ee FE she 


160 


greased to avoid leakage before putting into place. The 
pumps were started again, T1T3 switched to the bubblers 'BUBB', 
TAagecO es EXIT". and the reactor pressure during sampling 
was recorded from M2. T3 was then turned back to the 
normal outlet poattione NEED") and the needle valve connec- 
ted to T3 was properly adjusted to give a pressure reading 
in M2 the same as the value recorded previously. The 
reactor pressure would then stay unchanged in sampling, 
such that the steady state would not be upset. 

4) Ozonator Settings: Tl and T2 were switched back to 
the normal position 'REAC' with ozonator on. The typical 


settings were: 


Air Regulator Pressure see E yeksiel 
Ozonator: pressure dial 6 psig 
SCFH dial ¥ 
$ ozone LOOST Oreos 


Under normal conditions the reading in the first manometer 
Ml would be around 27 cmHg. A reading far less than that 
indicated leakage. 

5) Liguid Level Adjustment: The liquid level was 
brought to the mark on the reactor by ralayener ate @ lili lrstep uti oa 
was closed had the level been low and turned to 'SAMP' 
position if it was too high. The inlet pump was delicately 
adjusted until the inlet and exit streams were exactly 
balanced; then the liquid level should remain unchanged 


for an appreciable length of time. 
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6) UV_Intensity: The reactor could be operated at one of 
the four UV intensities listed below in an increasing order: 
i) No UV: UV switch OFF 
ii) Coarse filter:. 
iii) Fine filter: 
ivieruli UV; UV switch ON without’ 111ter 
The reactor was allowed to run for 45 minutes to reach the 


steady state, then the reactor pressure P_ and temperature T 


R 
were measured. 

7) Outlet Gas Sampling: The bubblers were filled with 
distilled water to the marks, 16g and 6g of KI were weighed 
into Bl and B2 respectively, which were then tightened into 
place. As soon as Dl was opened and the outlet gas sampling 
valve T3 was turned to 'BUBB' position in which the outlet 
ozone waS sampled, the stop watch had to be started and the 
cap in Bl pressed down firmly by hand, otherwise the reactor 
pressure was just high enough to pop it up. After 1.5 
Minutes (0.1 cubic foot) of sampling time, the watch was 
stopped, and T3 turned back to the normal '‘REAC' position 
simultaneously. Dl was then closed and the volume of the 

gas sample collected was directly obtained from the wet test 
meter by taking the difference of the initial and final dial 


readings. Finally, the ozone trapped inside the bubblers 


was completely flushed out with air by opening D2. 
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8 ) Exit Stream Liquid Sampling: The liquid sampling had 


to be done in sequence right after the gas sampling because 
a fast drawing out of liquid would affect the liquid level 
and thus upset the steady state. In the sampling process, 
T4 was turned to 'SAMP'and a sufficient volume of liquid was 
collected in a flask depending on the concentration to be 
analyzed. The flask was weighed before and after the sampl- 
ing step so that the amount of sample collected was known. 
The sample was then analyzed according to the procedure listed 
in individual experiments. 

9) Inlet Gas Sampling: The procedure was essentially the 
same as that of 7) Outlet Gas Sampling except that the out- 
let gas sampling valve T2 was turned to the 'BUBB' position 
while T3 was still connected to the reactor 'REAC'. The 
inlet gas sampling had to be done after the liquid sampling. 
10) Restart: The reactor was washed several times by 
fresh feeding solution with stirrer on. The same procedure 
was repeated from 4) on after changing the variables to 
be studied, for instance, stirring rate, inlet concentration, 
inlet partial pressure, pH and so on. 

tay Shut-down: ALCer the ozonator was switched off, air 
should still pass through the ozonator for another 10 
minutes to flush out all the residual ozone which might 
cause extensive damage to the ozonator. The reactor should 
also be washed several times by distilled water. The solu- 
tion containing electrolytes might stain the quartz plate at 


the bottom of the reactor which could hardly be rinsed off; 
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then the reactor had to be soaked with oxalic acid solution 


until all the stains were gone. 


Gas 


Phase Data Treatment 


The following terms are defined 

the total volumetric gas flow rate at T and P atm, 
L/min. 

the volumetric gas flow rate registered by the wet 
test meter after ozone stripping, L/min. 

the volumetric gas flow rate from the reactor gas 
SitieLpyet.e. at T.and Pr, L/min. 

total gas molar flow rate, mmole/min. 

the outlet molar flow rate of ozone, mmole/min. 

the inlet molar flow rate of ozone, mmole/min. 

the partial pressure of ozone in the reactor, atm. 
the total reactor pressure, atm. 

the atmospheric pressure, atm. 

the gas sampling time, min. 

the blank for thiosulfate titration, ml. 

the volume of thiosulfate solution required to titrate 
ozone in the outlet/inlet samples, ml. 

the mole fraction of ozone in the exit stream before 
stripping. 


the density of water at T, g/cc. 
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F _ volume of gas registered in wet test MeLeL Eee OTe Lye 


W 


Inlet/Outlet Ozone Molar Flow Rate n 


Af 
n.. = 


Al 


Since Vy is 


Aiea 


Since two bubblers were used to sample ozone, 


be 2V,° 


ie ; 
,min 


Ai’—ar 


(Vi - 2V,,) x: [Na,S,0,] 


2x0 


(V, - 2V,) x [Na,S,0,] 


2xt 


equal to 0.05 ml, 


(Vv, -— 0.1) * [Na,S,0,] 


BS iw 


(V,- 0.1) x [Na,S.0,] 


2 


is 


Partial Pressure of Ozone, , 


The ideal gas law is assumed. 


a ay NyaRT 
F. =F +x) = 
G W W Patm 
Piim'c = nRT 4 
In the reactor: 
PAF, n,RT 
PF nRvT 


(A9) 


(A10) 


(All) 


(Al12) 


(A13) 


In the wet test meter, 


(A14) 


(A15) 


(Al16) 


(A17) 
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Solving the above equations simultaneously, 


pressed by the measurable quantities on the 


the equation 
P : -1 
Py = Gee ee Pe 
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(A18) 


Pa can be ex- 
right side of 
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Appendix 3. Mixing Pattern 


Procedure 
The following instruments were used: 
Fisher Accumet, Model 520, Digital pH/ion meter 
Fisher Scientific Co., Micro probe combination electrode 
13-639-92 
Sargent recorder. 


Bid 


M (pH 3) HCl solution was prepared and stored in 
a 20 liter plastic bottle. The end of the exit tubing was 
clamped in an upright position, so that air bubbles would not 
be trapped inside. The glass electrode was first standardized 
against a buffer solution of pH 4.00 (Coleman Certified Buffer) 
and then slipped into the exit tubing through the end opening 
until the tip of the probe was completely immersed into the 
effluent. A vessel was placed underneath the exit opening to 
receive the effluent. The reactor was cleaned thoroughly by 
pumping through distilled water with the stirrer and the 
cooling water on. The rpm was measured by counting the num- 
ber of turns per unit time . The flow rate was set to around 
70 ml/min by adjusting the displacement pumps. The inlet 
pump was then adjusted carefully to maintain the liquid sur- 
face sharply at the marked level which was slightly below the 
middle stirring blades. 

An accurate measurement of flow rate was made by a 
graduate cylinder and a stop watch. The pumps were stopped 


and the inlet end of the tubing was dipped into the 20 liter 
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bottle containing 10°? MuHCL solution. As soon as the pumps 


were started again, the pH meter, the recorder and the stop 
watch were turned on. The pH values were read from the pH 


meter at the time intervals shown below: 


Interval for taking reading time 
30 sec 0-24 24min 
10 sec 2.2-4 min 
15 sec 4-6 min 
20 sec 6-8 min 
30 sec 8-11 min 
1 min 11 min and on 


The experiment was terminated at 30 minutes, and the PH of 


M HCl solution was taken down. The same procedure 


the 10- 

was repeated at the stirring rates of 0, 59.5, 85, 102 rpm. 
The volume of the reactor was obtained from the weight 

of distilled water required to fill an acetone dried reactor 

to the marked level correcting for the total dead volume 

Vpi * Yoo’ The inlet dead volume was determined by taking 

off the inlet line at the reactor, displacing the water 


trapped in the inlet system by pumping air and weighing the 


water pumped out. The outlet dead volume Voo was determined 


in a similar fashion. 


Data Treatment 


Residence time can be defined mathematically as 


ay = | taF (t,) = | toF(tp)dt, (A20) 
0 0 
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where F(t,) is the residence time distribution function. 
Physically, it can be viewed as the fraction of the effluent 
stream that has a residence time less than tpi and dF(t,) 
or F' (t,)dt, is the volume fraction of the effluent stream 
that has a residence time between tp and t,tdt,- 

In the case of a CSTR, the following equation can be 


established: 


dCs .3 


Orctetsi< it Cte 
where 
hae Coan" fh 
7, ee D1 Do (A22) 
L F 
L 

se the outlet acid concentration 
C the inlet acid concentration (step size) 


Co the background acid concentration in distilled water 


Fi the liquid volumetric flow rate 
V the reactor volume 
th, the time delay due to the finite volume of the tubing 
Di the inlet dead volume 
Voo the outlet dead volume 
i = F (A23) 


The solution of the differential equation with the given 


initial condition is: 
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C-C t-t 

_a = F(t.) = 1 - exp-(—— =5 (A24) 
Cc (oy R 

f£ fe) R 
and 

OF ead OH a ehh 
ln = L 

Gach = ies 

R 


In the experiment, the outlet pH was monitored at certain 
time intervals. A typical pH vs t curve at the stirring 
speed of 85 rpm is shown in Figure Al. If the reactor is 
an ideal CSTR, the response should obey Equation (A24) and 
the inet) vs t plot (Figures A2, A3, A4, A5) should 
be a straight line with slope equal to ieee: or -F,/V and 


aivalue of zero at t < th: 
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Figure A3. Plot of log(C,-C)/(C,-C.) vs time at the stirring 
rate of 85 rpm. 
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Figure A4. Plot of log (C,-C)/ (C.-C) vs time at the stirring 
rate of 59.5 rpm. 
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Figure AS.’ Plot. of log(C,-C)/(C,-C) vs time at the stirring 


rate of O rpm. 
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Appendix 4. Detailed Physical Mass Transfer Measurement 


Theory 
Direct Measurement - 
0,/H,0 System 
Assuming a zero gas resistance, the material balance 


on A leads to 


ceo cAS cA) Eton 
F_C 
: 2 GA i 
le Ky = eae cc = Ky, (A25) 
S (5- - Cy) 


Indirect Measurement - CO,/H.0 System 
According to the penetration model, Ky (CO,) can be 


related to Ky (0) by Equation (A26) 


RetO24 TO 
Osne3 2 
where the subscript A stands for O3- If the contact time 


depends only on the hydrodynamic conditions of the reactor 


and not the nature of the species, then On should be equal 


to 840 and the above equation can be reduced to 
2 
RaCCO,) Yo 
0 2 B 
COM Da | ee. 
Oras A 


The value of D can be predicted by Wilke and Chang correla- 
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The effect of temperature on diffusivity is expressed as 
Dj(T5)u(T.) DJ(T))u(T,) 


ai “earaevernea werent oes) LOLS PEL h ws G ghey 5 
T, Ty 


The temperature effect on viscosity is given by Swindell's 


PREreiation 
: 2 
u (T) 1.3772(20-T) - 0.001053 (T-20) 
(esp ae rn ee 
Bem 20°C) T+105 ee 


20 LO00PC 
where T is in degrees centigrade. 


From the above equations, Ky (0) at T, can be related to 


Z 
Ky (CO,) at T by 
3 
Oe  ) ages 1 Vp ,CO ie 
alee ees CA see 
2 cali ae | RCE ND bA 
As 04 is concerned, Equation (30) is reduced to 
Keel © 2%, Te) par) T 3 
Soa e [(—e op] (A31) 
i.) es aaa eR SS 


Data from the Literature 


_D 
(64) 
. = a 
CO, : Yoo, 34 ml/mole 
0, EOULINGe point Of 0, 112 .C see 
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e. Vua = 48/1.358 = 35.4 ml/mole. 


|O 


The values of D can be obtained from Wilke and Chang 
correlation or from exper imental work. The values from 


various sources are tabulated. 


Compounds Tee Belem sacnerteed Dalowe donciaes 
Wilke and Chang experiment 
CO. 25 Wipe 
20 6 1.76 HOG) 
18 ake yt ial tees ia 
0, 25 1.984 1.47 
20 1.805 1a oes 


Wilke and Chang correlation works very well for CO however 


ae 


it is higher than the experimental value by 30% when it is 


applied to O Since 04 and CO, are of comparable size and 


aie 2 
polarity, such a large discrepancy is unexpected. In 
Matrozov's cain the pH of the solution was not mentioned, 


and the author claimed that there would be no ozone decomposi- 
tion due to short contact time. Such a conclusion is obviously 
premature. There is no doubt that brevity of contact means 


ratio; however, 


an enhanced K and thus an increased tp/t 


Of D 
several important aspects have been overlooked. First of 
all, the slow reaction assumption owing to brief contact time 


is not justified unless the criterion t, << tp, is satisfied, 


and the knowledge of pH or the ozone decomposition rate Ks 
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is required to calculate tp: Secondly even the surface reac— 


tion is minimized, yet the bulk reaction still goes on as 
the ozonated water travels from the absorption chamber to 
the ozone monitor. The decomposition tends to lower the 
ozone content in the solution, and thus accounts for the low 
Da experimental value. As a matter of fact, the Da from 
Wilke and Chang correlation is used to calculate Ky (03) by 


the indirect method (H0-CO system). 


2 
2 
O 


3? 
H, (T) = Hy (T,) of 2870(T-T,) = 67680 + 2870(T - 20) (A32) 


the choice is supported by Shambaugh et al. in their recent 


Cae) 


work in which the solubility of ozone was calculated 


based on the value of Hy quoted from Perry's Handbook, i.e. 
67680 atm/mole/cc. The temperature coefficient of 2870 


atm/mole/cc per°C is from Metrozov's work (cf 2.2.1 Solubi- 


lity of Ozone in Aqueous Media). 


co.? 
The value of Hao (T) is obtained from Perry's Handbook 
2 
by interpolation 
Hag (T) = 1420+ 44(T-20) atm/mole fraction (A33) 
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Procedure for Direct Measurement 

The ozone partial pressure and the ozone concentration 
in the liquid phase were measured (cf. Appendix 2 for 
detailed operating procedure). The stirring rates of 61, 
75, 94, 112, 131 rpm were tried in the absence of UV. In 
step 8) of Appendix 2 Exit Liquid Sampling, the flask to 
receive the liquid sample contained 100 ml of 5% KI solution, 
so that 0, in the liquid sample was fixed immediately by 
reduction. The liberated I, was analyzed by iodometric 
method as illustrated in 4.3.1 Ozone Analysis, except that 
a microburette and a 0.005 N Na,S,0, solution were used to 
titrate the sample. The blank should consist of 100 ml of 


5% KI solution with the addition of the same amount of 1:10 


H5SO,. 
Data Treatment - Direct Measurement 
R ing ea of 0.005NNa,S,0, x [Na,S,0,] x density 
ey 2x weight of sample collected 
FEC 
ae ae (A25) assuming zero gas phase resistance 
=r 
S = 84 cm? 
H,= 67680 + 2870(T-20) (A32) 


FL= 66 ml/min . 


Ky can be solved. The assumption of zero gas phase resistance 
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(85) 


is valid as indicated by J. Rowley's work that the gas 


phase mass transfer coefficient is 1.53 cm/sec at 22°C. 


Procedure for Indirect Measurement 
The experimental procedure was essentially the same as 
that used for 03-H.,O except that no gas sampling was re- 


quired, because the partial pressure of CO, was directly 


2 
given by the total reactor pressure corrected for the 
saturated water pressure at the reactor temperature. The 
detailed procedure in Appendix 2 could be followed starting 
from 6) without UV; 7), 9), 11) were skipped. In step 8), 
exit liquid sampling, the flask used to collect the liquid 
sample should contain 20 ml of standard 0.02N NaOH solution 


delivered by burette, so that the CO, in the sample was 


2 
fixed immediately by neutralization once it was out of the 
reactor. The flask was weighed before and after sampling to 
give the weight of sample collected. Only 20 ml of sample 
was required, and then the solution was titrated for CO.- 
The stirring rates of 112, 102, 92.7, 86, 60.2 rpm and the 


blade positions of slightly above (normal), touching and 


1 mm below the surface were tried. 


(67,68) 


co. Analytical Procedure 


The following symbols are used: 


[OH] Concentration of NaOH solution as standardized by HCl 


solution with phenolphthalein as indicator, M. 
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[HCO 3] Background HCO, concentration in NaOH solution due 
to CO. absorption, M. 

VirV5 Volume of HCl required to titrate the sample to the 
first (phenolphthalein) end point and the second 
(bromocresol green) end point respectively, ml. 

V Volume of NaOH added to the sample, ml. 

Vy Volume of HCl required to titrate 125 ml distilled 
water to the phenolphthalein end point, ml. 


Circ Concentration of CO. in the sample as calculated by 


Vi and V5 respectively, M. 


C avg Average of Ci and Cor M. 

le) Density of water at temperature T, g/cc. 

wt Weight of sample collected, g. 

[Coo]. Background CO. concentration in distilled water, M. 
It has been known that NaOH absorbs CO, in the Suge SE) 


2 


so the standard 0.02N NaOH solution should contain a small 


amount of HCO, besides OH . The addition of NaOH solution 


to the sample to fix the CO 


2 would cause a background HCO, 


error which has to be corrected for. In order to deter- 
Mine the amount of the background HCO,, an indicator, such 
as bromocresol green which may change color at higher pH is 
added after all the OH has been neutralized to the phenol- 
phthaleéin end point, and the solution is further titrated 

to the bromocresol green end point. The volume of acid used 


for the latter end point corresponds to the HCO, in the 


standard NaOH solution. 
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Phenolphthalein OH + HY — HO 


Bromocresol green HCO, Eat H,CO, 


Volume of NaOH solution 
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Volume of HCl used for phenolphthalein end pointx[HC]] 


_Volume of HCl used for bromocresol green end point x[HCl] 


Volume of NaOH used 


A back titration method is used for the determination 


in the sample, because the co. dissolved in water at 


1 atm of CO. pressure may escape into air once the pressure 


2 


is removed in the sampling step. NaOH can neutralize CO. and 


fhus) stabilizesit in the solutionain an ionic form HCO, 


Before titration, the sample solution contains background 


HCO 


and the excess OH . 


3° 


and the HCO, formed by the reaction 


co. +. Otat— HCO, 


gives the amount of excess OH in the solution while the 


bromocresol green gives the total HCO3- 


(v_[OH ] - V, [HC1])° 
Pee a raepadt (A34) 
rT wt 


(V,[HC1] - Vv, {HCO,],)° 


M A35 
2 wt ‘ ( 


C ces Pee 4 ‘ (A36) 


Finally, the CO, background in distilled water is 
V,, [0H ] 


[con]. mI yi eon (A37) 
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Where 125 ml is the volume of distilled water to be titrated. 


Data Treatment CO,-H.0 System 
The Partial Pressur 
ure of CO., Eco, 
Assuming the gas phase in the reactor is saturated 


with water vapor at the reactor temperature T, then it can 


be shown that 


P = Pp p ; (A38) 
CO, R H50 
P. was found to be higher than Been by a pressure equivalent 
Pome Let nich HO (70) so 
a Fa =3 
P, = Bern ae db dk ewe Vodal HO = Piam eee PaTM ets ed Ee) 
Pi She 0.0241 aceesOcc 
2 
— ett = —_— = 
ta: ccs = eee O02 oer Usa O me aya 0.0214 (atm). 
The Interfacial co. Concentration [CO] . 
-1_ ©%% 
Xe = 5 (A39) 
Co 


where X. is the interfacial co. concentration expressed in 


mole fraction, Hao is the Henry's law constant for co. at 
2 


temperature T, atm/mole fraction. 


Hao = 1420+ 44(T-20) atm/mole fraction (A33) 
2 
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The Physical Mass Transfer Coefficient of cO., K,(CO,) 


Dhseeah (Cavge 7 102) 4) Fy, 
0 2 ([C0,] 5 Cayg? S 
5 (Cy ~ [co.] ) FL 
: (A41) 
10000 
Cc S) 
18H avg 
CO. 
> 10 BAS 
CO. 
The Physical Mass Transfer Coefficient of COzr Ky) (03) 
Ky (0,) = 0.988 Ky (CO) from Equation (A30) 
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Appendix 5. Detailed UV Characterization Procedure 


Through the years photcchemists have determined the 
quantum yield, $, of products of many photochemical systems 
as well as the intensity of the light beam by using the 
thermopile galvanometer combination ee Complex and tedious 
procedures are required in the calibration and application 
of the thermopile-galvanometer system, thus the trend among 
photochemists in recent years toward the use of chemical 
actinometers for light intensity determination is understan- 
dable. Besides convenience, results obtained with a well 
chosen, sensitive chemical actinometer are more reproducible 
and reliable than those obtained with the thermopile method. 

The intensity of a light beam passing through the 
reactor can be measured chemically by a suitable reactant A, 
which absorbs light of the wavelength of interest and for 
which the quantum yield, dar of some product B is known 
accurately for specific experimental conditions. Consider- 
ing the following photochemical reaction: 


iN NV B 


the quantum yield of B, dpe is defined as the number of 
moles of B produced upon the absorption of a single photon 
at a specific wavelength. Then reactant A can be placed in 
the photochemical reactor to be studied with the exact 
optical train etc. which is to be used in the photochemical 


reaction i.e. 0,/ferricyanide/UV reaction. The temperature, 
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pressure and the concentration of the reactant Aare adjusted 
to the appropriate range, and an exposure of A to the light 


beam is made for t sec. The following symbols are used for 


subsequent derivations. 


I the intensity of the light transmitted 

I the intensity of the light beam incident on the 
peaccor, 

nN, the amount of product B after exposure time t sec 

dp the quantum yield of B 

ue the exposure time, sec 

[A] the concentration of A 

€ the molar absorptivity of A 

x the light path cm. 


According to Beer's law, the fraction of the incident 


light absorbed by A is 


2-2 =1- 10 °!412 (A42) 
Oo 
“3 (A43) 
L =r 2068 quanta/sec! \s A43 
fe) pat (1 - 10 e[A] 2, 


For a reactant A with high molar absorptivity, where practi- 
cally all the incident light is absorbed, the equation can be 


reduced to 


a3) 


iy ee 


Since all terms on the right hand side are known, Ty can be 


I quanta/sec . (A44) 


calculated readily. 
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Potassium Ferrioxalate 
By far the best solution-phase chemical actinometer 
for photochemical reasearch today is the potassium ferri- 


C7 42s It 


oxalate system developed by Parker and Hatchard 
is very sensitive over a wide range of wavelengths and yet 


Simple to use. When sulfuric acid solutions of K Fe(C, 


3 4)3 
are irradiated in the range from 2500 to 5770 A, simul- 
taneous oxidation of oxalate and reduction of ferric to 
ferrous state occur. The quantum yields of Feau formation 
and the fraction of light absorption by ferrioxalate have 


(74) and are shown in 


been accurately measured by Parker 
Table Al. In the UV region at wavelengths up to 3660 A, 
almost 100% of the incident light is absorbed with an 
average Peo quantum yield of 1.24. The ferrous ion thus 
formed is made to be highly absorbing and easily analyzable 
by the formation of the red-colored 1,10-phenanthroline - 
ee? complex. 

The concentration of ee in the sample is read from 
the calibration by the net absorbance after correcting for 
the background i.e. the absorbance of the sample without 

2+ 


UV irradiation. The number of Fe produced in the reactor 


can be calculated by 


concentration of Rests LOO: M (A45) 
n ae a x SW er eee Oe [ 
B 20 
VY = 0.499 & . 
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According to Equation (A44), the number of quanta per sec is 
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Table Al. Quantum Yields of Fe** in the K,Fe(C,O,)3 Chemical Acti- 
nometer at 22° 


Fraction of Light Absorbed 


Wavelength, A [K3Fe(C,0,),], (J = 15 mm) pare 
5770 
5730] 0.15 0.118 0.013 
5460 es 0.061 0.15 
5090 - O15 0.132 0.86. 
4800 0.15 0.578 0.94 
4680 0.15 | 0.850 0.93 
4360 0.15 0.997 1.01 
0.006 0.615 1.11 
4050 0.006 0.962 1.14 
3660 0.006 1.00 . 1.21 
1.26 
0.15 1.00 1.15 
1.20 
3340 0.006 1.00 1.23 
3130 0.006 1.00 1.24 
2970 
a 0.006 1.00 1.24 
2537 0.006 1.00 1.25 


Table A2. Energy Distribution in Low- and Medium- 
Pressure Mercury Arcs 
Relative Energy 


Low-Pressure Medium-Pressure 
Wavelength, A Mercury Arc* Mercury Arc? 
13,673 Ae 15.3 
11,287 ser 12.6 
10,140 as 40.6 
5770-5790 10.14 76.5 
5461 0.88 93.0 
4358 1.00 77.5 
4045-4078 0.39 42.2 
3650-3663 0.54 100.0 
3341 0.03 9.3 
3126-3132 . 0.60 49.9 
3022-3028 0.06 23.9 
2967 0.20 16.6 
2894 0.04 6.0 
2804 : 0.02 9.3 
2753 0.03 2.7 
2700 tee 4.0 
2652-2655 0.05 15.3 
2571 ne 6.0 
2537 100.00 16.6° 
2482 0.01 8.6 
2400 Ley 7.3 
2380 aha 8.6 
2360 a 6.0 
2320 Sak 8.0 
2224 sa 14.0 


* Hanovia Lamp Division, Engelhard industries, Newark, 
N.J., SC-2537 lamp. 

» Hanovia’s Type A, 673 A, 550 w lamp. 

® Reversed radiation. 
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n. x 6.023 x 1029 


%, = 1.24 up to 3660A. 


The energy associated with the quanta flux I is 


Ix hx x 6.023x 107° 
E = 7 watt (A46) 
10 


where h is the Planck constant, and equal to 66256 % 1007" 


erg sec and v is the weighed average frequency of the medium 
pressure mercury lamp. Being a continuous light source, the 
medium pressure mercury lamp emits a spectrum covering a 
broad range of wavelengths. The average light frequency 

can be obtained only if the energy distribution function is 


known. 


) 

= 1 

pytatelwac 5 (A47) 
) 
1 


where oi is the relative energy at the frequency v, (g =100 

at Ay = 3650 ~ 3663). Only the wavelengths up to 5461 A 
needs to be considered; beyond this the quantum field is so 
low that it can be ignored. The average frequency is found 
to be 3310 A from the energy distribution function published 
by Hanovia Company ©! >”) Table A2. The power of the UV 


light incident on the reactor can then be written as 
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np x 6.6256 x 10° *7 x 6.023 x 1073 x 3x 101° 
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Procedure 


Synthesis of K 3Fe(C.0 ) 


ee 
Reagents Ferric chloride, Fisher Scientific Co., 
Potassium oxalate monohydrate, Baker Analyzed 


Reagent,J.T. Baker Chemical Co. 


The procedure was carried out under reduced light condition. 
600 ml of 1.5 M K5C,0, andi200: miniof .l.5 M FeCl, were mixed 
with rigorous stirring at 80°C, and the remaining undis- 
solved solids were filtered while the solution was hot. 

The filtrate was seated in the ice bath to be crystallized. 
The green crude crystals of ferrioxalate had to be recrys- 
tallized once by dissolving in 200 ml of water and repeat- 
ing the same procedure above. The crystals were dried ina 
current of warm air (45°C) by an air blower. The solid 

could be stored in the dark for long periods of time without 
change; however, a noticeable color change of the pure 


green crystals to a yellowish brown occurred rapidly, 


accompanied by crystal decomposition, when the crystals 


were exposed to direct UV light. 
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Preparation of Standard Solutions 


Reagents: Ferrous ammonium sulfate, ACS Reagent, 
The Nichols Chemical Co., Ltd. 
1,10-phenanthroline monohydrate, Certified ACS 
Fisher Scientific Company 
Sodium acetate, Certified ACS, Fisher Scientific Co. 
Sulfuric acid, Baker. 

The following solutions were prepared: 

OFO015 M Ferrioxalate solution: 29.5 g of KFe(C,0,) ,° 3H,0 
in a 4 liter brown bottle filled with 0.1N H,SO, 
solution. 

0.2% 1,10-phenanthroline solution: 0.4g in 200 ml distilled 
water. 

10 °MFe** standard solution: 0.426g of ferrous ammonium 
sulfate in a1 liter volumetric flask containing 
O2LN §H.S0O"> solution. 


2 4 


Buffer solution: 600 gm of Na0,CCH.,, and 360 ml 1NH,SO, 


diluted to l liter. 


Standard Fe-t solution for calibration: A burette was used 


to deliver 0, 2, 3, 4, 5 ml of standard 107° M 


rec solution to five 100 ml volumetric flasks. 
Hen ope em) Of 0.1NH.SO, solution were 


added to the above flasks respectively to make 
up the volumes to 5 ml. After introducing 2 ml 
of phenanthroline solution and 5 ml of buffer 


solution to each flask, the solutions were 
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diluted to volume by distilled water and 
allowed to stand for 60 minutes before 


measurement. 


UV Irradiation 

The reactor was sheltered by black paper and the 
inlet tubing wrapped by black tape. The reactor was first 
rinsed by distilled water, then filled up to the mark by 
0.015M ferrioxalate solution with T4 closed (referring 
figure 2) and all the fluorescent lights in the room except 
the one in the hood off. The stirrer and the cooling water 
were turned on and the exposure time was recorded as soon as 
the UV lamp was on. The exposure time would be 90 sec for 
100% UV, 100 sec for ‘fine filter', 150 sec for ‘coarse 
filter'and 150 sec for no UV. Then, the UV lamp and the hood 
light were both switched off and the solution in the reactor 
was collected in a beaker from which 20 ml of solution was 
pipetted into a 100 ml volumetric flask containing 2 ml of 
0.2% phenanthroline and 5 ml of buffer solutions. The 
flask was then stored in a dark place after filling up 


with 0.1NH,SO, solution. 


Spectrophotometric Measurement 


The absorbance was measured by a Bausch and Lomb 
Spectrometer Spectronic 20. The procedure had been des- 
cribed in reference 76. A blank solution containing no 


ferrous ion was used to adjust the 100% transmittance. The 
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samples had to he mee aras as fast as possible once they 
were taken out from the dark storage place, under a reduced 
fluorescence light condition to avoid further photo- 
decomposition of ferrioxalate. The calibration curve is 
shown in Figure A6, and the absorbances of the standard 


solutions are tabulated in Table A3. 


Intensity Determined by UV Intensity Meter 


Intensity can also be measured by a Blak-Ray Ultra- 
violet Intensity Meter, J-225 (for wavelengths from 2300A 
to 2700 A), Ultra-Violet Products, Inc., San Gabriel, 
Calif., U.S.A. The UV assembly was removed from the mount- 
ing pole, the UV intensities were taken by the meter at the 
center and the edge of the opening under various conditions 
i.e. 0~ 100% UV. The dimension of the opening was required 


inethe calculation. 
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2+ 
Run Fe added [re**] Transmittance Absorbance 
nt 10M Y, 

1 a yd! ee 11.0 0.9586 

2 4 4.348 Pao tee 0.6383 

3 6) cPyAcn A Dita 0.5560 

4 2 Pay Ee 42.0 0.3768 

< 0 0 100.0 0 


Stock Fe? solution concentration 1.087x10-M 


Table A3. Absorbances of standard ae solutions for the calibration 


of Spectronic 20 spectrometer at none 
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Figure A6. Calibration curve by plotting absorbance 
against the ferrous ion concentration :at 22°C. 
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APPENDIX 6 


Detailed Analyses of the Ozonation Reactions 


The symbols used in Tables A4-A8: 


KO 


KB 


cs 


CNO 


FRCN 


Physical mass transfer coefficient, eT 

Rate constant for the cyanate oxidation reaction, M 
Interfacial ozone concentration, M 

Calculated cyanate concentration, M 

Measured ozone consumption ratio 


Calculated ozone consumption ratio 


Reaction time, s 


Diffusion time, s 


Ratio of reaction time to diffusion time 


: -1 
Cyanide reacted per unit time, mole(s —) 


1 


; ; : -1 
Calculated ozone consumption per unit time, mole(s ) 


: } -1 
Measured ozone consumption per unit time, mole(s ) 


Fraction of ozone decomposed 


Fraction of ozone reacted with cyanate 


Fraction of 


cyanide destroyed 
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Listing of Program THES 


(i) Free Cyanide, Film Model, Flat CNO Profile, pH 7 and 12 
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Si bea T No ,efDtord OW 201% .tebot ait 5 


| 
| Se P X TANG?) 1 3eme eres 19.35 210 
yay FOR THES eva 5.19 037 
*INCS(CARN,1443 PRINTER) 
REAL KNAT,KN 
DIMENSTON GP(60),GT(60),G6CNI(60) »GCNE(60 GCNO(6 GKR a) 
NIMENS TON ERB eNEOT CKAI CEO) CENTS AMAR TT RAT MeO tent cueieas 
NIMFENSTOAN GTR(60) GTD(60) »GTRD( 60) sGRN (60) »GRC (60) GCS (60) 462 (60) 
DIMENSTON GRO (60) 9MM( 60) ¥GKO(60) »GINT (60) yGKL (60) 
DIMENSION GK4(60),GKR4(60) 
DIMENSTON GKL2(60)9GKL3(60) ¥GKL4(60) ,GEN3(60) »GEN4 (60) 
DIMENSION GN3(40) 
DIMENSION SO0(69),S1(60) 
Pe tAh sles) /1 20*0./ 
DATA GN32/60%*0./ 
NATA GKL 2 sGKL3,9GKL45GEN3, GENG/300*0./ 
DATA GK49GKR4/120%0./ 
DATA GKR39GTR sGTDyGTRD»GRNyGRC GCS yGRyGRN,GKN,GINT»GKL/720%0. / 
DATA GP GT sGONT sGCNF,GCNN,GKRL yGKR2,GKR3yGENyGKI1 9 GK2/660%0, / 
NATA MM/60%0/ 
KK=0 
I JK=1 
FL=./65 
GO TN 1 
997 FL=1.094 
KK=0 
ITJK=TJK+1 
] RFAN(542)MyP yXNT9XN1 9 T,CNI gCNF 
? PORMAT{(I1;,6F10.5) 
TF (M-5)3,9,9 
3) KK =K K+1 
pict 
XKN=15.3 
Bla 4e25 
ARA=84, 
RN=(XNI-XNN)*.LO01L/ 60. 
R=(CNI—-CNEF )%*.O01XEL 
S$ =7 0302586%*(1.32772%(20-T)—-2.001053%( T-20 ) *( T-20))/(T4+105.) 
W=1L.002%FXP(S ) 
DP=(14+273.2)*1-805E-05/ ( 293. 2*1)) 
KNAT=.184*SORT(NT/1 «8 389F-05) 
KN=KNAT/ARA 
H=67680 o+2870.2%*(T-20. ) 
CS =P/H*1000. 
CALL MULR(XsXKReCNF NT KOGCS 9R pROGCNT 9 XKN GIT) 
meek UK) =X 
GKR1(KK )=XKB 
GKL (KK )=2.15°%*RN/(2.15*ARA*XCS/10N00 --RN) 
GEN(KK)=GKL(KK)/KQO 
TT=IT+1 
CALL MULR (UX sXKRgCNF gDT,KONGCS 9R gRO GENT, XKDN GIT) 
GK2(KK)=xX 
GKR2 (KK )=XKB 
GKL2(KK )=1.498*2.15*R/(2-.15*ARA*®CS/1000 --R¥1 2.498) 
GINT (KK )=GKL2 (KK) /KA 
TJT=II+l 
CALL MULR(XyXKRyCNF DT yKO9CS 92 gROGCNT y XKN GIT) 
GK3A(KK)=X 
GKR3(KK)=XKR 
GKL3 (KK) =160 200% 261 5%*R/ (261 5*ARA*CS/1000 .-R*1 0200) 
GEN3 (KK )=GKL3(KK)/KA 
Ti=lI+1l 
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CALL MULR(XsXKByCNF yDTyKOy¥CS gRyROQCNT 9 XKN GIT) 
GK4(KK)=xX 

GKR4(KK )=XKR 
CNON=1./(1./(CNI-CNF )+X KR/X/CNF) 

TR=1 6/(X*CNF+XKN+STXCNOXXKB ) 

TN=DT /KN/KN 

TRN=TR/TN 
RC=CNO*XKB*ST/ ( XC NFFX KNF STXKCNO*KXKR ) 
RN=XKN/(X*CNFE+XKDN4+S TXCNOXXKR ) 
GKL4(KK)=SORT (NDT ( X*CNF+ X KB S TCNO+ XKD) ) 
GFEN4 (KK )=GKL4(KK)/KN 

GN3(KK)=R%*( X*CNF+XKN+S Te XKBECNT) /X/CNE 
GCNT (KK )=CNI 

GCNFE(KK)=CNE 

GP(KK)=P 

CGCNN (KK )=CNNA 

GKN(KK)=KN 

See kKKiy=CS 

GTR(KK)=TR 

Crate ak l=1) 

GTRN(KK)=TRD 

GRN(KK)=RN 

GRO(KK)=RC 

GR(KK)=R 

GRN(KK)=RN 

SO(K K)=GRN (KK) /GR (KK ) 

S1(K K )=(X%*CNEF+S THX KR*ECNO+ XKN) /X/CNF 
MM(KK)=M 

Sire LOL 

ILINF=0 

imesh tsi + KK 

TLINFHTLINE4+1 

Pen E11) 1O0;01,10 

IF (ILINF-42)12,11,12 

WRITF(6513) 
PORMAT('1!.///// 9 OX oI RUN! 6 3X%e*tUV' 4 5X_ "CNT! 5 9X_'CNFI,9X, 


pee nee oes 81! 5X, tRO!) - 


WRITEt(6.111) 

FORMAT(® *,18X%o!M'g11X og *Mtgl1iX, "ATM! ,8X,'DEG ClolX-'*MOLE/SS') 
WRITE(6 520) 

FORMAT (/) 

IL INF=ILINF+1 

TFC ILINF-42)75,11,7 

IP=I-1 

ITF (GCNI(IP)-GCNI(1I))697,6 

IF (MM(T)-1)14,14,15 

WRITE (6524) T »GONT (I) sGCNF(1I),G6P(1) »GT(T),GRO(T) 
FORMAT (*# 845X129 3Xo9IND' y5X 9 3( F110 4591X) 9 F5225,E11 04) 
iste Lie 3) 

WRITE(6529)I1e¢GCNI(1I)sGCNF(T) »GP(T),GT(T),GRO(T) 
meen! foe. 1 203Xe "FULL saxo 3Cbll st, 1X) ,F5.2,EL Fo) 
CONTINUE 

TioO=0 

11 =9 

ITEST=90 

NO 60 JT=1,KK 
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ARM 


110=1 1041 
IF (1-1)32, 35,3? 

3? IF (1-8) 33548 y 33 

33 IJIP=I-1 
IF (110-11 )374 35535 

37 ~~ ~=IF(GCNI(IP)-GCNI(T))38, 39,38 

38 WRITF(6,34) 

34 FORMAT(/) 
ITEST=ITFST#H1 
GN TN 39 

35  WRITE(6,36) 

Be FORMAT( NL /S////47Xy "RUN NOe y4X_ Al gL 2Xy Bl gL 2Xq IC 44 12Xy ID?) 
110=0 

39 WRITE(6440)TyGK1(1) »GK2(T) 9 GK3( 1) »GK4(T) 

ME GRMAT (8s 855X591 255X%5o 4K! 4 3X, 30E116452X) 561104) 
WRITE(6541 )GKB1(1),GKR2(T) 9 GKR3(T),GKR4(T) 

4] FORMAT(* %,12X,"KB*, 2X ,3(E11.4,2X) ,Fll. 4) 
WRITE (6945 )GKL(I) GKL2(1)yGKL3( 1) »GKL4(1T) 

o> BeemAT(' "51 2X_'KL!, 2X 5 3(E11.4,2X) »Fll. 4) 
WRITE(6547 )GEN(T) GINT(1),GEN3( 1) ,GEN4(T) 

a7 FORMAT (' *,12X,'EN', 2X 53(E11.492X) »F11.4) 
IF (TTFST—3)515 52452 

5? TF (ITFST-7)5 3554554 


51 i ate, 
GN TN 60 
me «61d =8 
GO TN 60 
el 1 =7 
GN TN 60 
42 WRITF(6,55) 
55 FORMAT (* ',///,6Xy'TARLE THE COMPARTSON OF RATE CONSTANTS AND M 


*ASS TRANSFER',/,10X9*COFEFFICTENTS ORTATNED RY 4 DIFFERENT WAYS', 
P/>10X-s'A RT FROM GAS PHASE MASS BALANCE! ,/10X,'B RT=1.5%R', 
Pen eee oe eh —1.2%*R',/.,10%,'N KR=K/5.1, RT FROM R#) 
eer 10 35 
60) CANTINUE 
J2=-1 
41 S2=J527+1] 
T1=14+J2°% 
1 2=2+J2% 
132=3+J2% 
14=44+J52% 
15 =54+J2%5 
TF (JS 2-J2/2%2)65562965 
6? WRITF(6572) Tl 912913914915 
72. = FORMAT (194 /S/// 14 7X 9 RUN NOG 2X9 4( 1241 2X) 9 12) 
G0 1TO~79 
65 WRITE(6975 )T1L,T291391 4915 
75 FORMAT (' ty /S////47Xo!RUN NO’ '92X94( 129 12X) 972) 
79 WRITE(6581)GKN(T1),GKO(12) »GKN( 13) ,GKO(14) »GKO(T5) 
81 FORMAT (! ty / 4g 7X y*KO' y TXe4( EL] 049 3X) yEL1 4) 
WRETE(6 +83)GCS(11),GCS(12)+GCS(13),GCS(14) GCS (15) 
83 FORMAT (# §.6X%y*CS' se TX 4(Fll oe 49 3X) oFll o4) - 
WRITE (6984)GCNO(T1) GCNO(T2) sGCNO(T3)9GCNO(T4) »GONO(T5 ) 
84 FORMAT (t 1,6Xyt CNN, OX 9 4( E11 0 4,3X) Elle 4) 
WRITE (6987 )GTR(T1)9GTR(12) eGTR( 13) 9GTR(14) 9GTROT5) 
87 FORMAT (' 'y6X%e4 TR! 9 7X9 4( FLL 649 3X) FLL 4) 
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WRITE (6 488 JGTD(T1) sGTN(12) yGTD( 13) ,GTD( 14) gGTD(T5) 
88 FORMAT (* "56X9'TN'47X94(E11.453X),F11.4) 
WRITE (6289 )GTRD(I1) sGTRN(12) -GTRD(13),GTRD (14) »GTRD(IS) 
89 FORMAT (! "s6X9"TRD" 4 6X 94( E11 .453X) 5Fll. 4) 
WRITE (6491 )G03(T1),G603(12),603( 13) ,603(14) ,G03(15) 
91 Pena rUPN YS 6Ks"03%", 6X 54(E11.4,3X) fll. 4) 
WRITE(63995)S0(11)+S0(12),S0(13),S0(14),S0 (15) 
95° FORMAT(' ',6X,!'SO ", 6X 94(E11.4,3X) Fl. 4) 
Bee (Gy 997)S0(T1),S1(12),51(13),S91(14).S1 (15) 
Meer ORMAT(* '56X>,'S1 ',6X ,4(F11.4,53X) F114) 
WRITE(6 990 )GR(T1) »,GR(1T7?),GR(13),GR(T4) »GR(I5) 
90 FORMAT (* *,6X5'R',8X,4(E11.4— 3X) »Fl1.4) 
WRITF(6392)GRD(T1),GRD(12),GRD(13),GRN(T4) ,GRN(TS5) 
9? BRMAT (* "'56X>'FRD!',6X .4(E11.453X) ,EL1. 4) 
WRITE(6593)GRC(1I1)9GRC(1T2)yGRC( 13) ,GRC(14) »GRC(T5) 
92 Seem arT Att 2, 6X5 ' FRC! 56% 540611 24,3X),F11.4) 
TF (J2—KK/5 )61 461,999 
999 TF (TJK-1)997,997,998 
99R STAMP 
END 


FEATURES SUPPORTED 
ANE WARD TNTFEGFRS 


TOCS 
CORE RFOUTREMENTS FOR THES 


COMMON 0 TNSKFEL COMMON 0 
VARTARLES 4100 PRNGRAM 2018 


ENN) NF COMPTLATTOAN 


THES 
DMP FUNCTIAN COMPLETEN 
J// FOR FMI) fama 7 15.22.33 
SURROUTINE FMU( TOL» XK Re CNF a9 DT KO CS gReFeCNI ¢ XKD gROGTT ) 
Seat i. KO 
ST=3 625 
ARA=R84, 
Pea ti—2)1,2-3 
] RT=RN 
SO¢7TH410 
2 RT=1.498%R 
Brie tO. 10 
3 TF(ITI-3)5,5,6 
5 RT=1.2°«R 


10 XKBR =1 O0/(ST/(TOL*CNF*(RT/R-1 2) -XKN)—1. /TOL/CNF )/ (CNI-CNF ) 
CNO=1e/(16/(CNI-CNF )+XKR/TOL/CNF) 
6010; 7 

6 XKR=TOL/5.1 
CNN=1./(1./(CNI-CNF ) +X KB/ TOL/CNF ) : 
RT=R*( TOLECNE+STEXKRECNOFXKD)/TOL/ON 

7 F=2 15*RT/(2015*ARA/1LONO « *CS=RT)-SORT (DT *( TOL *CNF+ST*XK BC NOFX KD ) ) 
RETURN 
END 
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FEATURES SUPPORTED 
ANF WORD INTEGERS 


CORE RFOUTREMENTS FOR FMI 
COMMON 0 INSKFL COMMNN 0 
VARTARLFS 16 PRNGRAM 216 


END NF COMPTL ATION 
EMI! , 
AMP FLINC TIAN COMPLFETFEN 
J// FOR MIILR Per tiet a. he 5624603 
PeeerNe MULR(X, XKBeoCNF,DT»KO,CS >Re ROVONT, XKD,I1) 
REAL KN 
XEP=Looo. 
X¥RT=10N0NN0N. 
FPS=.01 
TEND =32000 
TFR=O 
4 ee a 
FRaXRT 
X=XL 
FPL eX 
CALL FMU(TOL,XKR GCN FyNT KO 9CS Ry FyCNTy XKNgROGIT ) 
TF (F-0.)1,14 91 
1 FL=F 
X=XR 
TOL=xX 
CALL EMU TAL sXKRyCNFgNTeKOQCS gR go FoCNTy XKN gRAG TT) 
TFE(F—-90)291469? 


2 FR=F 
; PEA GTGN(1.»FEL)+SIGN(1.,FR)) 25,3925 
2 T= 
TALFHLON.*FPS 
4 T=14+1] 


NN 12 K=1l-TEND 
X=.5%* (XL +XR ) 
TAL=xX 
CALL FMUCTOL yXKRGCNFANTZKOGCS sR,eF CNT, XKD RN, TT) 
TF (F-0.)551655 
5 TF(STGN(1 ey F)FSIGN( Ley FR)) 79690 
6 TOL=XL 
XL=XR 
eA EAL 
TOL=FL 
FL=FR 
FR=TAL 
7 THUSFHFL 
A=F*T OL 
A=A+A 
TF (A-FR*(FR-FL))89999 
R IF (I-IEND)1791799 
9 XR=X 
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FR=F 

TOL=FPS 

A=ARS (XR) 

ITF (A-1.)11511,10 
10 TOL=TOL*A 
11 meee XR—-XL)—-TOAL) 12512513 
12 IF (ABS(FR-FL)-TOLF)14,14,13 
13 CONTINUE 


TFR=1 
14 TF(ARS(FR)-ARS(FL))16,16,15 
we X=XL 

F=FL 


146 RETURN 
D7 A=FR—-F 
DX=(X—XL)*FLH( Le +F(A-TOL)/(A*(FR-FL)))/TOL 
SM=X 
FM=F 
X=XL-NDX 
TOL=X 
Peuk EMUI TNL sXKBeCNENTsKOeCS »ReFeCNI, XKD,RO,I1) 
IF (F)18,1465,18 
18 SERPS 
A=ARS(X) 
Pe({A-1.)20.20519 
19 TOAL=TOL*A 
20 Me CARS UNX)—-TNL) 21, 21,22 
7 PE CARS(F)-TOLF)16,16,22 
2? TF (STGMN(1.9F)+SIGN(1.9FL) ) 24923524 
a3 XR=X 
FR=F 
GN TN 4 
24 XL=X 
FL=F 
X¥R=XM 
FR=FM 
GA TN 4 
25 TFR=? 
RETURN 
FND 
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COMMON 0 TNSKFEL COMMON 0 
VARTARLES 26 PROGRAM 580 


END NF COMPILATION 

MUIR 

DMP FUNCTION COMPLETED 

// XFO THES maemo 78 15.24.04 
*CCEND 
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(ii) Ferricyanide, Penetration Model, Flat CNO Profile, 


pH 7. 
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*INCS (CARD, 1443 PRINTER) 
REAL KO,KOAT 
EXTERNAL FCT 
DIMENSTON GCNI(&80O) 
DIMENSTON GKP(&80) 
DATA GCNI/&80*0,/ 
DNATA GKP/80*0./ 
KK=0 
READ(552)M,yP,T,CNI,CNF,PH 
FORMAT(1I1,5E10.5) 
IF (M-5)2,9,9 
2 KK =K K+ 1 
ARA=84, 
Bea sho 75 
XKB=3470. 
QH=1.2E-14/10.%% ( —-PH) 
XKD=700. 0H 
R=(CNI—-CNF)*.OOL*FL 
§$=2.2302585%*(1.3272%*(20-T) -—. 001053%(T-20) *( T-20))/7(T+105) 
WeaLoog2rEXP(S) 
Pen 12773522) *1.805E-05/( 293..2%1)) 
ROA n=, 184*SORT(NT/1.8389E-5) 
KO=KNAT/ARA 
H=6 TORN +28706.*%(T-20.) 
ESsP/H*1000. 
Ye b=100000. 
XbIl=1. 
EP S$=.0000 1 
TEND=1000 
NMEMU=2 
CALL MULR (Xs F,XLI,XRIg FPS TEND, TERS CNE SDT 9KNyCSyR CNT, XK NEM! J 
CKP(KK)=X 
GENT (KK) =CNI 
CUrLo 1 
9 ITLINF=0 
PO 409 bale KK 
Pint he=1) INE+] 
IF(ILINE-1) 31,32, 31 
21 TE(TLINE-44) 23, 32, 33 
37 WRITFE(6,35) 
25 SOPMAT (I 1', /////+10X%s 'RUN NO.) 10X, 1 KE) 
37 WRTTE(6,38)1,GKP(T) 
3R FORMAT (' stk pl ceo elle) 
Teihe=TLINE+ 1 
TE(ILINF-44) 40,4 34, 40 
24 ReELMPHTLINE-1 
an TO 40 
aS [P=L—-L 
FEO ObP —GONT(1))372 39927 
39 WRITE(6536)1,GKP(1T) 
36 FORMAT (! beOX%s b2H15XsE11.4) 
40 CONTINUE 
STOP 
END 
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INCS 


CORE RFOQUTIREMENTS FOR PEN 
COMMON 0 INSKEL COMMON a) 
VARTARLES ZR4 PROGRAM 434 


END OF COMPILATION 
PEN 
DMP FUNCTION COMPLETEN 
// FOR MULR EIMSBeMT THTOVSTOESR 
SUBROUTINE MULR (XSF sXLT_ XRT ¥FEPSy TEND, TER sCNE (DT gKNGCS 9R gCNT 9 XK gh) 
FEAR PKG 
hER=O0 
rt ee ae 
XR=XRI 
=X i 
TOL=xX 
(et ts FMU( TOL, CNE,ZAT SKN, CS PeFe CNT, XKD,N) 
befr—-0.)1,16,1 


1 FL=F 
X=XR 
fOe= xX 


CAEL FMU( TOL, CNF,NT, KN CS_P 9F 9 CMT». XKDy, NM) 
ProGr-O0 «) 25 1652 


2 PReF 
PeUStTGN(1].,.FL)+SIGN(1.,FR) 125,3,25 
2 t=0 
TOLFH1LNNV EPS 
a pie] 


bieers. K=1L, PEND 
haope tal +XR) 
EOL =x 
CALL FMU( TOL,» CNF DT, KM, CS9R Fy CNT 9 XKN,N) 
Petr. )55 16.5 
5 IF (STGN(1e,F)+STGN( Leg FR) IF oAe7 
6 TOL=XL 
XL=XR 
eee etl. 
peices 
-L=FR 
FR=TOL 
7 TOL=F-FL 
A=FTOL 
A=A+ A 
TF (A-FR*(FR-FL))8599,9 
8 ITF (T—-TENND) 17,1799 
9 yR=yX 
FR=F 
TOL=FPS 
A=ABS(XR) 
TF (A-1.)11,11,10 
10 TOL=TOL*A 
a. TE (ABS(XR-X¥L)-TOL)12512,13 
2 TF (ABS (FR-FL) -TOLF)14,14,1? 


aaT3 yauny a nT 
f,.T£,@1 TY 942 Be Ve a “AN 
PAST 2945 TAK eS IX, HK) SAI Tie ame 


wy, i ep tee AD NAST 9°43...) Lie 
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rae 


13 CONTINUE 


TER=1 
14 Pee oor n ABS (FL )16,16,15 
oe. X=XL 

FHFL 


16 RETURN 
17 ASFR-F 
DX=(X—XL)*FL*(1.+F*(A-TOL) /(A*(FR-FL)))/TOL 
XM=x 
FM=F 
X=XL-NX 
TOL=xX 
ap fel FMU ( TAs CNFSDT, Kf» CSek 9 Fas CNT, XKDsy Ni) 
PeCr ats, 16418 
18 HOUSER S 
A=ABS(X) 
PLAST.) 20520,19 
ee TOLETOL*%*A 
20 Peer oti x) —TOLS2i~e 21s 22 
21 Pepa (PIL) 16s L622 
i? PePSTGN(T.5F)+SIGN(1.5FL) )24,23,24 
eiP xR=X 
f=P=—f- 
BO rTON4 
24 X|_ =X 
FL=F 
aa ai 
FR=FM 
ar TH) 4 
ae TER=2 
RETURN 
END 


EFATURES SUPPORTED 
OME WARD INTEGERS 


CNRE REQUIREMENTS FOR MULR 
COMMON ®- INSKEL COMMON 0 
VARTARLES 26 PROGRAM 916 


END OF COMPILATION 

MULR 

DMP FUNCTION COMPLETEN 

// FOR FMU MenSeh oT? LOeS Ts 55 
SUBROUTINE EMU( TOL, CNFE DT» KN,CS, Re Fe CNT, XK yN) 
REAL KO 
FEXTFRNAL FCT 
ARA=84. 
XKR=3470 
CNO=1./(1e/(CNI-CNF )+XKB/TOL/CNF) 
RTH( Let (XKN+ 36 25% CNN*XKB) / TALS CONF) *R 
Wl L=XK BE CNOX 36 25+ TOL*XCNF+ XKN 
W QaW 1*4,%NT/3.1416/K0/K0 
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XU=SORT(W2) 

IP S=1.0F-10 

NIM=100 

CALL INTG(Oe,9XU,TPSsNIM,FCT,FER ) 
TE=ARA*XCS*SORT(W1L*ENDT) (FER (Let. 5/W2)tFFXP(—W2)/SORT(3B.1416%W2) ) 
TE=TF*.OOL1 

P= Te—RT 

RETURN 

END 


FEATURES SUPPORTED 
ONE WORD INTEGERS 


CORE REQUIREMENTS FOR FMU 
~ COMMON QO INSKEL COMMON 0 
VARTABLES 34 PROGRAM 22 


Seow COMPILATION 
FRU 
NMP FUNCTION COMPLETED 
een (15 SEP 77 19.38.24 
SUBROUTINE INTG(XL+XtleFPSeNDIM, FCT, Y) 
DIMENSTON AUX( 100) 
DNATA AUX/100*%0./ 
Pee )=e Set Ch Xid+ FET (XU) ) 
H=XU=-XL 
TF(NNIM=-1)8,8,1 
x PRG 2, 10,2 
2 HH=H 
E=EPS/ABS(H) 
BELT =0'. 
P=]. 
JJ=1 
DIRE Vs f=2,N0 EM 
Miz ALXs( 11s) 
DELTIS—E La 
H[)=HH 
bH=.5%HH 
P= .52P 
X=XL+HH 
SM=0. 
NO 3 hePh rere 
ShM=SM+FCT(X) 
Ss] axe D 
AUX (1 pee SKAUX CI -1)+P*SM 
Q=1. 
hal sh 
NO 4 J=1,JI 
Li=I -J 
Q=0+0 
Q=0+0 
4 AUX (TL) =AUX CLIF1) 4 (AUX (TT+1 ) -AUX(TT))/(9-1.) 
NELTZ=ABS(Y—-AUX(1)) 
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MAL Oily eS 
CP IDEW T2Z—E) 105.105, 6 
Pte be Writ) ts Ll Ls LL 
JJ=JII+ JJ 
LER=Z2 
9 Y=Hx AUX (1) 
RETURN 
10 IER=0 
GO. TH 9 
V1 TER=] 
Y=aHouy 
RETURN 
END 
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PEATURES SUPPORTED 
ONE WORD TNTEGERS 


CORF RFOQUIREMENTS FOR INTe 
CAMMON QO INSKFEL COMMON 0 
mee TAR LES 2320 PROGRAM evra 


END AF COMPILATIAN 

INTG 

eee FUNCTION COMPLETED 

pee FOR FCT Poyocrodgd) Lo .3R. 42 
FUNCTION FCT(XF) 
POT =Tei Zk SR*E XP ( —XF=xXF ) 
RETURN 
END 


FEATURES SLIPPNRTEN 
ANE WORD INTEGERS 


mee REGQUTREMENTS FOR FCT 
COMMON QO INSKFL COMMaAn 0 
VARTARLES 4 PRNGRAM 26 


FND NE CAMPTLATION 
FCT 

DMP FUNCTION COMPLETEN 

// XFO PEN 15 SEP 77 19.38.54 
*CCEND 
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(iii) Ferricyanide, Film Model, Flat and Quadratic Profiles, 


pH in 


>a 


t esi Part 


oftesbau) baa 


ae a 
¥7> JCB M Pp X TANG Miser hth be 51) oe 2 
yy er OR THES bOGwece 1! 16.50.24 
*ITNCS(CARD,1443 PRINTER) 
REAL KO,KL,yKOAT 
DIMENSTON GCNI(80) »GCNF(80) ,GP( 80) ,GT(80) »GK (80) ,GKO(80) ,GNT( 80) 
DIMENSTON GKO(80) »,GKL(80),CEN( 80) ,GCS(80),GCNN(80) »GTR(80) ,GTD(80) 
DIMENSTON GTRD( 80) »,GR(80),GRT(80),GFRN(80) ,GFRC(80) »,GFRCN(80) 
NIMENSTON MM(80) 
DATA GONT »,GCNF4GP,G6TGKy GKO »GNT,GKN, GKL »GEN/800*06./ 
DNATA GCSyGCNO,GTR»GTDNsGTRI »GR_yGRT»GFRD »GFERC »yGFRCN/800*0./ 
DATA MM/80%*0/ 
KK=0 
READ(5,2)MyP_9ToCNI ,CNEF 
? FORMAT (I 1,4F10.5) 
| py Oi cee are BB 
3 KK=K K+ 1 
AR A=84, 
FL=.4944 
XKR =3 ZenTi Os 
XKN=8.5F-04 
Re (Ch T—-CNF)*.001*FL 
$=2.302585*(1.3272%( 20-T) —.001053*(T-20) *( T-20) )/ (T4195) 
W=1.002*EXP(S) 
mrt 036 2) 1 BUDE -O5/(29.7. 251) 
Fp. Lh4sSORT(DT/1.5389E-0 5) 
KN=KNAT/ARA 
H=6 7680.42870.*(T-20.) 
GS=P/H* 1000. 
XRT=100000. 


— 


XLI=l. 

Pro se G00 1 

TEND=1000 

NEMUS 1 

CALL MULR(X Fy XLI9XRI 9 FPS, TENDg IER GCNF (DT 9 KN 9CS4R CNT y XKD yNFMU) 
GKO(KK) =X 

NFMU=2 


CALL MULR(X%,FoXLI_XRI 9 FP SoT ENDg TER 9 CNF QNT gKMyCSyRyCNT y XK 1D NFM ) 
CNO=1./(1./(CNI-CNF)+XKB/X /CNF ) 
RT=(1let( XKN+XKB*CNO* 3.25) /X¥ /CNF)*R 
KL=RT/ARA/CS* 1000.4 

EN=KL/KO 

GK (KK) =x 

RR=R*1000./CNI/FL 
TR=le/(X*CNF+ XKN4 36 25% CNM XKB ) 
THEDT /KO7 KO 

TRO=STRATLIO 

RCHCNOXXKB* 36 25/ (X*CNF+ XKN 26 25% CNA XKB ) 
RN=XKKN/ (X*CNE+ XK D+ 3-6 25% CNO *XKB ) 
GP(KK)=P 

GT(KK)=T 

GCNI(KK)=CNI 

GCNF (KK )=CNF 

MM(KK)=M 

GR(KK)=R 

GRT(KK)=RT 

GKL( KK) =KL 

GKO(KK)=KO 

GDT(KK)=DT 

GEN(KK)=EN 
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GCS(KK)=CS 
GCNO(KK)=CNO 
GFRCN(KK)=RR 
GTR(KK)=TR 
GTDN(KK)=TD 
GTRD(KK)=TRD 
GFRC(KK)=RC 
GFRD(KK)=RN 

GQ TO 1 

TL INF=0 

DO 30 T=1, KK 

IL INE=IL INF+1 
BRACTeLNE 7 1))) 10.5), V1) 10 
IF (ILINE-42)12,11,12 
WRITE(6,13) 

FORMAT (91%, /////46Xo RUIN No hy 4X oF UV gy BX yg CNT! GL OXy HCNE 4G LOX, 


ae MME EO Xe UT ats) 


WRITE (6,20) 

FORMAT (/) 

TLINF=IL INF+1 

DEMWIBLNE —42: 7 9th Light 

bP = 1b) 

IF (GONT (IP )-GCNT(1))65756 

TF (MM(T)-1)14,14,15 

LEQMM-CL) =2:))16% 16,5 17 

TRUM CL —3) 1B, LA pl9 

WRITE(6524)1,GCNT(1),GCNF(1T),G6P(1),G6T(T) 

FORMAT (* "45 XyT2y8X_*NON' 4 6X ,F1164,2X%y Fl o4y2Xy Fl 1 0452 X%yF 502) 
Gir Fhe 30 

WRITE(6,26)1,GCNI(1I) »GCNF(T ),GP(I),GT(T) 

FORMAT( 1 ,5X eT 2¢8Xe*COARSE',2X,FLL eo 4e2Xe Elle 4e92X5FlL1 049 2X yp F502) 
Gatti 730 

WRITF(6,28)1T,GCNI(1I),GCNF(1T),GP(1),G6T(T) 

FORMAT(" %,5XyT298XgtFINE', 4Xyg Fl Lo 4y 2X gFlL lL e4y2XgFll oho 2XyF 502) 
GO TO 30 

WRITFE(6,29)1,GCONI(1T) sGCNF(T ),GP(1),G6T(T) 

FORMAT (! 5X eT 2e 8X eo tFULL® p4Xe FL 1 0492X pFLL oe 2X gE) 049 2X 9F522) 
CONTINUE 

TLINE=0 

NO 40 T=1,KK 

ILINF=JLINF+1 

PreesINE—1) 31,3231 

PROT LINE—44) 33-32-33 

WRITE (6,35) 

FORMAT(' 1's /////y10X9* RUN NOR ty LOXe*K! 20K yt KO!) 

WRITE (6,5 38)1,GK(I),GKO(T) 

FORMAT(* %4/,10X912915X-Fl16499X,F 11.4) 

IL JNE=IL INE+1 

TLINE=IL INE-1 

GO TO 40 

R= E=<1 

IF(GCNI(1IP)-GCNI(1))3 7939 27 

WRITE (6, 36)1yGK(1) GKO(T) 

FORMAT(* §,9X_912915X9F116499X9F 11.4) 

CONTINUE 

J2=-1 
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41 J2=J2+1 
Tl=1+J2%5 
T2=2+J52%5 
13=3+ J2*5 
T4=4+ J2%5 
Peab+ J2*5 
IF (J2-J2/ 2% 2) 455424545 
42 WRITE(6,52)11,172,13,14,15 
52 FORMAT (9 1"%5/S////47TXg*RUN NO O49 2X41 12, 12X) 912) 
SUeTO 59 
45 WRITE(6,55)11,12,13,14,15 
55 RPURMATC® 8 o/////47Xe* RUN Nile 2K s 402th weit eu 
59 WRITE(6,61)GKL(11),GKL(12), GKL( 13) »GKL(14) ,GKL(I5) 
61 BORMAN CISL edXo * KL, 7Xe4( Fl1.453X),611.4) 
WRITE (6,62)GKO(1I1)+GKN(1I 2), GKN( 13) ,GKO(I4) »,GKO(I5) 
6? FORMAT(!? ",6X9"*KO',7X,4(Fl1645 3X), Fl1-4) 
WRITE(6,63)GEN(I1),GEN(12), GEN( 13) >GEN(14) »,GEN(I5) 
62 FORMAT(* *,6X,'EN',7X,4(F110.45 3X) ,F1I1.4) 
WRITE (6,64)GDT(11),GNT(12), GOT( 13) sGNT(14) ,GNT (15) 
64 BORMATCP# II OXs' DT, 7X,4(F11.453X)eE11.4) 
WRITE(6,65)GCS(I1),GCS(12), GCS(13),GCS(14) ,GCS(15) 
65 FORMAT(* ',6X%,'CS',7X,4(F11.453X) ,F11.4) 
WRITE(6,66)GCNO(I1),GCNN( 12 ),GCNO( 13), GCONO(1T4) ,GCNN(T5) 
66 FORMAT(* *,6X_5'CNO!,6X%,4(F11.24,3X),E11.4) 
WRITE(6,67)GTR(I1),GTR(I2), GTR( 13) eGTR(T4) ,GTR(I5) 
67 PUMA t(’ §'66X5'TRYEgTXe4(F11.453X%),5F1L.4) 
VRITF(6,68)GTD(I1),GTN(12), GTN(13),G6TN(14) ,GTD(15) 
6B FORMAT(*® *,6X,'TD*',7X,4(F1124,3X),E11.4) 
WRITE(6e69)GTRN(I1) sGTRD(I2 ) »GTRD( 13), GTRD(T4) ,GTRD(IS) 
69 FORMAT(* ',6X%,'TRD! 46X_o4(F11.24,3X),F11.4) 
WRITE(6,71)GR(1I1),GR(1E2),GR (13) ,62( 14) sGR( 15) 
aL FORMAT(* #,6X9!R!',8X,4(F11 453X%),F11.4) 
WRITE(6,72)GRT(I1),GRT(12), GRT(13)+GRT (14) ,GRT (15) 
Pe FORMAT(* ",6Xy,'*RT" yg 7X40 Fl 1049 3X) F114) 
WRITE(6,73)GFRN(1I1),GFRN(I2 )»GFRD( 13), GFRD(T4) ,GFRD (15) 
<3 FORMAT(! #!,6X,'FRD! ,6Xe4(E11.4,3X),Fl1-4) 
WRITE (6,74) GFRC(1I1),GFRC(I2 ),GFRC( 13), GFRC(T4) ,GFRC(15) 
To& FORMAT(! 1,6X%e'FRC!',6Xe4(F112-4,3X),Fl1L.4) 
WRITF(6,75)GFERCN (I11),GFRCN( 12) ,GFRCN( 13) ,GFRON(14) ,GFRCN(T5) 
75; FORMAT(' *,6X%_,"'FRCN!',5X_4( Fll.493X)9F11 24) 
IF (U2-KK/5) 41,41,999 
999 STOP 
END 


FEATURES SUPPORTED 
ONE WORD TNTEGERS 
IOCS 


CORE REQUIREMENTS FOR THES 
COMMON @) INSKEL COMMON 0 
VARTARLES 3368 PROGRAM 16° 
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DMP FUNCTION COMPLETED 
// FOR MULR POpoe eit? 16.51.53 
SUBROUTINE MULR(UXsF yXLIyXRT,EPSy TEND, TFRyCNF DT oK09CS yRyCNI gy XKN NN) 
REAL KO 
IER=0 
XL=XLI 
XR=XRI 
X=XL 
TOL=X 
CALL FMU( TOL, CNF DT, KN, CSR 9 Fe CNIy XKDy N) 
PROR—-O.) 1; 16,1 - 
7 FG=F 
X=XR 
TOL=X 
CALL FMU(TOL, CNF DT, KN, CS, Ry Fy CNI 9 XKD GN) 
FR R00.) 2516,2 
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Z FR=F 
PSICSTGM(1L.sFLI+FSIGN( 1. 5FR)) 2553,25 
2 1=0 
TOLFH=100.%FPS 
4 T=[]+1 


DO 13 K=1,ITEND 
X=.5%(XL+XR) 
TOL=X 
CALL FMU(TOAL, CNF DT KM, CS 5 Ry Fe CNI» XKNQN) 
leit =-0.)5-16,5 
. ITF (STGN( Leg F)+STIGN( LeeFR)) 19657 
6 TOL=XL 
XL=XR 
XR=TOL 
TOLS=FL 
FL=FR 
ERe rE 
qT TOL=F-FL 
A=F*# TOL 
A=A+A 
TE (A-FR*(FR-FL))89999 
8 IF (I-ITEND) 17,1739 
9 XR=X 
PRA 
TeG=EeS 
A=ABS(XR) 
FE CAs.) 11,11, 10 
10 TOL=TOL*A 
aoe IF (ABS(XR—-XL) -TOL) 12,12,13 
ed TE(ARS(FR-FL) -TOLF )14,14,13 
14 CONTINUE 


TER=1 
14 IF (ABS(FR)-ABS(FL)) 16916915 
15, XeXb 

F=FL 


16 RETURN 
17 A=FR—-F 
DX=(X—XL ) *FL*( 1+ *(A-TOL)/ (A*¥(FR-FL)))/TOL 
¥M=xX 
FM=F 
K=XL=~-DX 
TOL=X 
CALL EMU(TOL, CNF yg DT9KO_CS_ Re Fy CNT gy XKD9N) 
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IF(F)18,16,18 
18 TOL=EPS 
A=ABS(X) 
IF(A-1.)20,20,19 
19 TOL=TOL**A 
20 ERGABS (0X) -TOL P21 21,22 
21 IF (ABS(F) -TOLF) 16,16, 22 
2? Pee PON AL eet S1GN( 16 5FL.)) 24.23.24 
no XR=X 
FR=F 
GO TO 4 
24 XL=X 
FL=F 
XR=XM 
FR=FM 
GO TO 4 
25 TER=2 
RETURN 
END 


FEATURES SUPPORTEN 
ONE WORD INTEGERS 


CORE REQUIREMENTS FOR MIUILR 
COMMON @) INSKFEL COMMON 0 
VARTARLES 26 PROGRAM 576 


END OF COMPILATION 

MULR 

NDMP FUNCTION COMPLETED 

// FAR FMU LOSE Pr doe 20) 
SUBROUTINE FMU( TOL, CNFsDT»K Oe CS¢Re Fe CNI yXKD GN) 
REAL KO 
REAL 3KL 
XKB=3470. 
ARA=84,. 
FL=.4944 
GCNO=H=1./(1-/7( CNI —-CNF )+XKB/TNL/CNF) 
RT=( Let (XKN+ 36. 25*CNO*XXKB) / TOL/CNF)*R 
KL=RT/ARA/CS* 1000.6 
IF (N-1)1,1,2 

1 UU 1=3 .25%*FL*CNO*XKO/ARA/DT/( 46 25* TONL*CNF+ XKD ) 

UWU2Z=SORT((4.25* TOL CNFF XKD NT) 
UWUZ=(EXP (UI2/KO) -EXP(-UU2/K0))/2. 
FH=UU2*UUL/SCS/UU 3+ (CS-UUL) *tNU2/CS/TANH(UU2/K0) -KL 


GO TO 3 

2 UL=SORT(DT*( TOL*CNFE+ X KM 36 25% CNN%* XKB) ) 
F=U1/T ANH(U1L/KO) -KL 

a RETURN 


END 
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FEATURES SUPPORTED 
ONE WORD INTEGERS 


CORE REQUIREMENTS FOR FMU 
COMMON O INSKEL COMMON 8) 
VARTABLES 26 PROGRAM 220 


END OF COMPILATION 

FMU 
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